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PAPERS 


STATISTICAL ANALYSIS HYDROLOGY 


SYNOPSIS 


The basic principles the duration curve are reviewed this paper 
such manner outline the various ways which properly used 
hydrologic analyses. Some the details and differentiations discussed are 
small significance compared inaccuracies existing the greater part the 
hydrologic data available today. However, the proper use the duration 
curve essential fundamental understanding statistical analysis and, 
since proper use more difficult than improper use, the differentiations are 
deemed relevant from engineering standpoint. 


CONCEPTION THE CURVE 


its application flood-control work, the duration curve may looked 
upon refinement the practice basing the design the largest flood 
record. The difficulty this old practice lies not only the uncertainty 
that the maximum occurrence representative for the period record, but 
also the absence criteria for determining the variation magnitude with 
period occurrence. Both these deficiencies are remedied the proper 
use the duration curve. 

The theory the duration curve relatively simple. The occurrences 
the future (those that are for the future) are considered constitute 
series, arranged the order descending magnitude. there are in- 
finite number occurrences store, then this series may represented 
continuous curve such that shown Furthermore, considered 
that there cycle trend magnitude with respect time, may 
assumed that occurrences which happen within given period have each been 
selected random from the store occurrences. Lastly, since the number 
occurrences store considered infinite, the selection finite number 
will have effect the distribution magnitudes within the 
remaining series occurrences. 


comments are invited for immediate publication; insure publication the last dis- 
cussion should submitted February 1943. 


Asst. Engr., Engrs., War Dept., Los Angeles, Calif. 
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Since they have been chosen random from the store occurrences, the 
occurrences record (those used the analysis) may used esti- 
mate, the laws probability, the curve represented Fig. With this 
curve defined, forecasts concerning future occurrences may made. 


CONSTRUCTION THE DURATION CURVE 


The curve shown Fig. may represent all the occurrences ever 
happen may represent all the yearly maximum occurrences ever 


Relative Number Occurrences 


Magnitude Occurrence 


CuRVE 


happen, any one many other series. This curve called probability 
distribution curve and similar nature the normal probability curve 
Gauss and Laplace. Equal areas under the curve represent equal numbers 
occurrences. The duration curve the integral curve the probability 
curve. Thus, magnitude occurrence plotted against cumulative area 
under the probability curve, the duration curve obtained. Fig. shows the 


Relative Number Occurrences Smaller 
Than Given Magnitude 


Magnitude Occurrence 


Curve (UNCONVENTIONAL) 


duration curve that corresponds the probability curve shown Fig. 
The maximum value the duration curve ordinate equal the total area 
under the probability curve. this value called 100%, then the ordinate 
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seale Fig. will range from zero 100%, and this scale called the per- 
centage occurrence scale. Convention interchanges the ordinate and ab- 
scissa scales and reverses the order the percentage occurrence scale such 
that the maximum ordinate corresponds occurrence. The duration 
curve Fig. shown conventionally Fig. 


Percentage Occurrence 


Fic. (CONVENTIONAL) 


constructing the duration curve from the occurrences record, 
assumed that each occurrence representative for the period—that is, the 
largest occurrence 10-yr period considered the median the largest 
occurrences all the 10-yr periods ever happen. Likewise, the second 
largest occurrence 10-yr record considered the median the second 
largest occurrences from all the 10-yr periods ever happen, and forth. 
Since the total number greater than any median equal the 
total number smaller than the median, expected that half the occur- 
rences will exceed their medians and that half will smaller than their medians. 
each occurrence plotted the position its median, smooth curve 
drawn through the points will the best available representation the true 
duration curve. 

Fig. all the occurrences excess the 10% occurrence value 
constitute 10% the total number occurrences. one occurrence chosen 
random, the probability that will exceed the 10% occurrence value 
10%. Similarly, the probability that occurrence chosen random will 
occurrences are grouped random sets ten, then the probability that the 
maximum occurrence set chosen random will exceed the median the 
maximum occurrences all the sets 50% (since there are many above 
median there are below). Therefore, may stated that the probability 
that none the ten occurrences chosen random will exceed the median 
the maximum occurrences all the sets 50%. 
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Let equal the per cent occurrence the median. Then p), 
will the probability that occurrence chosen random will not exceed the 
median. The probability that ten consecutive occurrences chosen random 
will smaller than the median therefore and known equal 50%, 
0.50. From this, the value 0.933, and equals 0.067. The percentage 
occurrence corresponding the median the maximum occurrences from 
all the sets therefore 6.7%. 

similar manner, possible determine the percentage occurrence 
that defines the median all occurrences consisting the second maximum 
occurrence from each set. this case, the probability that least nine 
ten occurrences chosen random will fall below given value 50%. The 
probability that none will exceed the value The probability that ex- 
actly one will exceed the value Then the probability that least 
two ten will exceed the value one (certainty) minus the two foregoing 
probabilities and known equal 0.50. Thus, 


and, therefore, 


TABLE solution successive approximations 


For TEN yields value 0.836 for and 0.164 for 
OccURRENCES The percentage occurrence corre- 
sponding the median the second 
largest occurrences from the sets there- 
fore for the third largest, 
16.4 and others, may computed similarly. 
25.8 
35.5 The “plotting position” occur- 
54.8 rence its estimated percentage oc- 
64.5 
74.2 currence. Plotting positions for any set 
ten arbitrarily chosen occurrences are 


given Table 


INTERPRETATION THE DURATION CURVE 


plotting the duration curve, the assumption was made that the recorded 
occurrences were representative any set events that would occur pe- 
riod equal the period the record considered. This not true, however. 
The largest event may twice great would normally occur the given 
period, and the second largest may smaller than the normal second largest 
occurrence. Throughout the entire series occurrences, however, these dis- 
crepancies should counterbalance each other large extent. The longer the 
record, the more accurate the balancing out. Since the series plotted 
not accurately representative, will not plot smooth curve. the points 
that appear high were lowered somewhat and the points that appear low were 
raised such manner that they all lie smooth curve, the points then 
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would much better representation normal series occurrences for 
period equal the period record. plotting the curve, the largest ten 
occurrences was assumed equal the median the maximum occurrences 
all ten sets. the occurrence plotted did not lie the smooth curve, the 
assumption was error, and the point the smooth curve having the same 
percentage occurrence was the new estimate the median. 

Suppose that the occurrences for the future arbitrarily divided 
into sets ten occurrences each. one these sets chosen random, 
there 50% chance that the maximum occurrence will exceed the median 
the maximum occurrences all sets. This obvious since there are many 
occurrences choose from excess the median there are smaller than the 
median. Thus, the median the maximum occurrences all the sets the 
smallest value that probably will not exceeded single set chosen ran- 
dom. This value for ten-occurrence set read the smooth curve 6.7%. 

Since the smooth curve obtained the best possible representation the 
true duration curve, and since the curve independent the number occur- 
rences used its construction, normal series comprising any number oc- 
currences can obtained from the curve. Thus, normal series one hun- 
dred occurrences can obtained merely picking off values the respective 
points which each one hundred occurrences would plotted (their 
plotting positions). 

often has been necessary determine the frequency with which given 
value exceeded. the duration curve determined that the value 
question has 10% occurrence, then, the average, 10% the occurrences 
will exceed that value. Thus, fifty occurrences happen ten years, five 
can expected exceed the value. 


CURVE 


The necessity extrapolating duration curves short-period records 
make long-period determinations has arisen almost every hydrologic in- 
vestigation which duration curves were used. Unless the curve can de- 
scribed mathematically, the extrapolation necessarily based upon few 
occurrences the end the curve, and, therefore, liable large error. 
the form the curve can determined beforehand, however, every point can 
play its part determining the extrapolation. 

Many attempts have been made describe the duration curve mathe- 
matically. many equations suggested, only one has rigid mathematical 
derivation. This the Gauss-Laplace normal distribution curve. The funda- 
mental assumption upon which this law has been formulated that the arith- 
metic mean the variate equal the median the variate. The justifica- 
tion this assumption for many cases has been proved unquestionably 
innumerable observed distribution curves. However, the attempts that have 
been made fit rainfall runoff data the curve have resulted definite 
discrepancies. The result that many individuals have taken upon them- 
selves find more suitable distribution curve use rainfall and runoff 
analysis. Unfortunately, the ordinary treatment the problem has not been 
attempt discover the discrepancies the Gaussian law and thereby cor- 
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rect the fundamental assumptions, but rather adapt the law imposing 
limits using mathematical variations the equation. Most these not 
adhere the simple assumptions that form the basis the law Gauss. They 
not have rigid derivation, therefore, and will fit observed data only 
accident. 

Outstanding discussions the problem here involved appeared the 
Transactions the Society 1927? and Particularly enlightening 
the discussion Arne Fisher.‘ Referring Thiele, Mr. Fisher states: 


mentions that, the distribution the variate, actually 
observed, found to.be skew, deviate from the normal probability 
curve, may considered ‘fictitious’ observation the ‘pseudo- 
normal’ variate, h(x), which latter function supposed follow the 
Gaussian error law; that is, h(x) distributed according known 
normal error law exponential form, possible find the frequency 
distribution the inverse function, 


Again Mr. Fisher states: 


investigations led him the obvious fact that the geo- 
metric mean the observations yields the most probable value, the arith- 
metic mean the logarithms the observations must yield the most 
probable value the logarithms the observations, that the log- 
arithms the observations are assumed normally distributed, 
law error can found for the observations themselves logarithmic 
transformation the variate. Such law must necessarily lead the 
geometric mean the observations the most probable value.” 


logarithmic transformation the variate often has been used with ex- 
ceptional success. However, certain happen when short-period records 
are analyzed, many variations from the ideal curve were observed. was 
not apparent that there was theoretical justification for using the logarithmic 
transformation, was considered that other curves could used which would 
have good basis. However, substitute law seldom has been formulated 
that did not receive more adverse criticism than praise. Evidence supporting 
refuting the applicability the normal distribution curve rainfall and 
runoff analysis weak. indicate the applicability inapplicability the 
normal distribution curve rainfall and runoff occurrences would necessitate 
obtaining good many 1,000-yr records. 

When series rainfall occurrences consisting every occurrence record 
plotted, the existence the annual cycle may, and probably does, impart 
degree the series, since cannot stated that each occur- 
rence independent every other. the series selected consists the 
maximum occurrence each year, however, the series not affected the 
annual cycle. Such series, having the advantage normal distribution, 
often can used instead the series consisting every occur- 
rence recorded. necessary construct duration curve all occur- 


rences, the duration curve yearly maximum occurrences may used 
“Straight Line Plotting Skew Frequency Data,” Goodrich, Transactions, Am. Soc. 


Asymmetric Probability Function,” Slade, Jr., ibid., Vol. 101 (1936), 35. 
Ibid., 72. 
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extrapolate the curve extrapolation appears difficult, since one curve 
function the other. 

The writer has constructed duration curves for six rainfall records from 
southern California stations varying length from forty-nine sixty-two 
years. Four duration curves 24-hr rainfall for these stations are shown 
Fig. The data were plotted logarithmic-probability paper (on whicha 
duration curve will straight line the logarithm the variate follows the 
normal probability curve). Some good results have been obtained, and some 
sizable deviations exist. This, however, constitutes neither proof the 
applicability the law nor condemnation the law. The records are far 
too short serve basis for testing the law itself. They do, however, 
illustrate the typical duration curve and show what type deviations from the 
normal may expected. 


THE DURATION CURVE 


The writer has attempted analysis the accuracy obtained using 
duration curve. After choosing six long-period rainfall records stations 
the vicinity Los Angeles, Calif., duration curves yearly maximum, daily 
rainfall totals were prepared for the entire period each station. addition 
these, the periods were broken down into 10-yr, 15-yr, 20-yr, 25-yr, 30-yr, and 
40-yr periods, and duration curves for each period for each station were con- 
structed. The average period record for the stations was about sixty years. 
Each station, therefore, had approximately six 10-yr records, four 15-yr records, 
three 20-yr records, two 30-yr records, and one 40-yr record. Duration curves 
were drawn straight lines logarithmic-probability paper. 

The first step the analysis was determine for 10-yr period whether 
the 10-yr normal maximum occurrence computed the duration curves had 
greater correlation than did the actually recorded 10-yr maximum occurrences. 
The values obtained from the duration curves for given station were averaged, 
and the percentage deviation each value from the mean was computed. 
The thirty-three deviations thus obtained from the six stations were used 
compute the probable error single observation. similar procedure, 
the probable error single recorded 10-yr maximum occurrence was com- 
puted. The probable errors thus obtained were 12.6% for duration-curve 
value and 17.5% for recorded value. 

The probable errors for values the longer periods could not computed 
satisfactorily the same manner since there were not enough occurrences for 
each station establish good mean value. However, probable errors for 
these periods were computed assuming that the correct value for given 
station that indicated the duration curve for the entire record 
for that station. Deviations were computed percentages this long-term 
value, and probable errors were then determined. These are shown Table 
2(a). will noted that both the probable error and the largest deviation 
experienced the computations are consistently higher for recorded values 
than for computed values. Furthermore, the period record lengthens, 
the probable error computed value greatly reduced, whereas that 
recorded value reduced only small amount. apparent, then, that 
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one the advantages design the duration curve rather than the 
maximum occurrence record that greater accuracy obtained. 

provide better indication the variation accuracy with the period 
record, the probable errors single 100-yr value, single 50-yr value, 


(a) BEAR LAKE DAM 


ao ~ 


> 


Ww 


24-Hour Rainfall, Inches 


- 


Percentage Occurrence 


and single value were computed for the same periods record 
Table 2(a) and the same manner. The results these computations are 
shown Table 

Table applies only curves having approximately the same slope the 
curves Fig. the slope the curve greater, the probable errors will 
greater. good index the slope the curve the ratio the value 
the 99% value the curve. This ratio may used analysis relate 
probable error the slope (degree dispersion) the duration curve. 


MAGNITUDES FOR DESIGN 


the case ordinary flood-protection structure, the period for which 
the structure afford protection chosen economic analysis. this 
the case, the structure should designed control the normal sequence 


MOULTON RANCH 
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floods expected within the given period. the economic side 
secondary importance, may found advisable design the structure for 
abnormally severe sequence floods that have certain degree prob- 
ability occurring during the life the structure. The extreme case that 


24-Hour Rainfall, inches 


S 


0.6 


Percentage Occurrence 


Los ANGELES, CALIF. 


which the structure designed for the maximum possible flood which can 
take place within the lifetime the structure. 

loosely used term this regard occurrence. 
This term not definable unless one states the degree probability intended. 
the ordinary sense the word, means greater than 
chance happening greater than 50% probable. However, some engineers 
think ‘‘maximum probable” much greater than this. One attempts 
conceive something possible, but having probability happening. This 
impossible since zero probability synonymous with impossibility. 
would seem advisable view this misunderstanding eliminate the 
phrase, from the hydrology vocabulary. 

deriving the values shown Table the writer assumed that the value 
sought was that which has 50% chance being exceeded within given 
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period. However, the value sought, for instance, that which has 25% 
chance being exceeded one hundred years, then, instead using possibly 
value having 50% chance being exceeded in, say, two hundred years 
(as often erroneously done), value can derived easily that will determine 
flood having 25% chance being exceeded one hundred years. This is: 
100(1 0.287. Thus, not only inadvisable, but also 
unnecessary, dealing with structure designed for 100-yr protection, 
consider any frequency other than the 100-yr frequency. 


TABLE AND OBSERVED 
PROBABLE OCCURRENCES 


CoMPUTED AND RECORDED (%) 


Years Probable Largest 
record Normal Error (%) Deviation (%) 100-yr 
occur- No. 
rence curves 
Duration Re- Duration Re- 
corded curve corded 


protect dam against possible failure, would ideal, feasible, 
design spillway for the maximum possible flood flows. However, found 
that this quantity, estimated conservatively, huge prohibit design. 
his paper entitled ‘‘Maximum Probable Floods Pennsylvania 
Am. Soc. E., states: term ‘maximum probable flood’ 
used describe not the largest flood possible, but flood large that the 
chance its being exceeded greater than the hazards normal all 
man’s This definition has been criticized being too vague. 
However, the basic idea seems good one and, the probability de- 
fined can put down percentage probability such 10% (one chance 
ten), then the flood under consideration can determined mathematically, 
although large probable error may exist such remote point the dura- 
tion curve. This large probable error may cause some hesitate use the 
quantity, but there type analysis available which the error can 
Therefore, may stated that so-called ‘‘maximum probable” 
flood can determined the use the duration curve the degree prob- 
ability occurrence within given period time stipulated. 


his paper entitled Alden Am. Soc. 
E., states: Having only limited number items from which plot con- 


Transactions, Am. Soc. E., Vol. 106 (1941), 1453. 


50-yr l-yr 
; ; 20 20 15 11.2 15.6 37 42 14.4 13.0 4.5 
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tinuous curve, the assumption made that this continuous curve would pass 
through the center each block the actual apparent that this 
assumption slightly error, since the values thus obtained for Table would 
15, 25, instead 6.7, 16.4, 25.8, ete. number items increases, 
however, the error this assumption rapidly becomes negligible, and curves 
plotted the approximate method will very close those plotted the 
correct plotting positions. Since computation plotting positions for various 
numbers occurrences requires considerable time, has been impossible 
include table plotting positions. For the same reason, may found 
impracticable for individual compute plotting positions for certain 
studies. this case, good approximation may made computing the 
plotting position the largest occurrence and estimating the plotting positions 
the remainder the occurrences, using for guide the approximate plotting 
positions referred Mr. Foster and expressed the equation: 


will noted that the true plotting positions are always closer 50% than 
are the approximate ones, and that the difference between the true and ap- 
proximate values becomes smaller 50% approached. 

using the duration curve for flood-control design, the design flood 
usually defined that which will exceeded once (on the average) within 
the lifetime the structure. Since com- 
mon sense dictates that flood which TABLE 


probably will exceeded ordinarily OccURRENCE VALUES 
poor design flood, appears advisable 

Life structure Percentage occur- 
adopt the design flood definition the (years) rence (design flood) 
minimum flood which will probably not 
exceeded within the life the struc- 
ture. Table gives percentage oc- 
currence values (one which was deter- 
mined the foregoing analysis) design 


floods for various periods for which struc- 
tures may designed. 

has been common practice, order facilitate analysis, plot only 
the upper portion duration curve for which great number occurrences 
are available. From engineering viewpoint, this procedure generally 
justified because approximately correct and eliminates great deal work 
that most cases unnecessary. However, the construction curve from 
partial series has two disadvantages—plotting positions are approximate, 
and partial series cannot have normal distribution. These disadvantages 
may some cases (in which more accurate work justified design considera- 
tions and available data) outweigh the advantage expediency. 

Probably one the most important weaknesses the duration curve re- 
sults from its use mechanical process, the accuracy which often over- 
estimated and sometimes underestimated. One therefore hesitates use the 
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curve for record short duration. However, addition the recorded 
occurrences, possible use all types rational comparisons and analogies 
assist constructing the curve itself. this can done properly, very 
satisfactory design can obtained. Any method used hydrologic design 
must flexible and adaptable varying quality data. Hydrology has not 
advanced far enough but that judgment must still play major part any 
analysis. 


Conclusions drawn from the preceding analyses and discussions are 
follows: 


(a) order represent best distribution curve using finite number 
occurrences, the occurrences should not plotted with equal spacing; 

Following from conclusion (a), the plotting series occurrences 
above base value complete series not mathematically rigorous; 

(c) The use the Gaussian normal distribution curve has way 
been proved inapplicable properly chosen rainfall runoff series; 

(d) The duration curve can used determine the occurrence that prob- 
ably will not exceeded within given period well that occurrence which 
will exceeded (on the average) once within given period; and 

(e) The ‘‘maximum flood, adequately defined, can deter- 
mined use the duration curve. 
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PERMEABILITY MUD MOUNTAIN DAM 
CORE MATERIAL 


WALTER,? AND HOWARD HARSTAD,* JUNIORS, 
Soc. 


During the early investigation the material contemplated for use the 
impervious core the Mud Mountain Dam (in northwestern state), 
was proved that the moisture content the material when compacted had 


much greater influence the permeability than did the density the 


material, within the range densities which would compacted dur- 
ing construction. Since the voids the soil act conduits for the water, 
postulated that the water content the soil the time compaction 
governed the arrangement the voids. high water content, the voids 
were apparently all uniformly small size and the permeability was low. 
alow water content, but the same numerical void ratio density, the voids 
were apparently wide range sizes and the permeability was high. Thus, 
wide narrow range void sizes may likened wide narrow 
range particle sizes. 

That the wide range narrow range void sizes may imparted the 
soil will, with absolutely variation density being necessary, proved 
the extreme range permeability among groups samples compacted 
very nearly identical densities but different moisture contents. 

This paper reports the results investigation inaugurated determine 
the permeability-voids ratio relationships the proposed borrow-pit material 
basis for the distribution the material the various zones the Mud 
Mountain Dam. pronounced were the effects slight moisture differences 
the permeability the material that moisture control during construction 
comments are invited for immediate publication; insure publication the last 


discussion should submitted February 1943. 


Associate Geologist, Chf. Seattle Dist., Soils Laboratory, Engr. Dept., Seattle, Wash. 
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became paramount importance, not much for the usual density control 
but primarily for control the permeability the core material. 
Throughout the paper reference made water and 
order accommodate those who may not conversant with 
this terminology, the following definitions are given: 
Optimum water water content which the maximum 
Void ratio the volume intergranular space the 
volume solid particles given soil mass without regard the pro- 
portions liquid, air, gas which may occupy the 
This paper will discuss turn: (1) The material tested, (2) and 
equipment, (3) test results, (4) interpretations, and (5) further investigation. 


MATERIALS TESTED 


Test results this paper apply only the materials used the impervious 
core the Mud Mountain Dam. The soil mixture stream gravel and 
glacial till. The average mechanical analysis the material after removing 
all boulders larger than in. shown Fig. Curves and represent 


Sieve Sizes 
0.75" 0.25"4 No. 


Curve Till 
Admixture (Inches) 


Grain Size, Millimeters 


ANALYSIS 


20% and 40% till admixture, respectively. The curves reveal extreme 
range particle sizes from colloids in. 

The colloids are found primarily the stream gravel and occur naturally 
coating around the particles the deposit. They have been identified 
the National Bureau Standards ray and petrographic analysis 
montmorillonite, alteration product volcanic glass and ash. Because 
its extreme flatness, montmorillonite has been described 
The extreme sensitivity the mineral excess moisture results 


Mechanics Manual Engineering Practice No. 22, Am. Soc. E., 1941. 
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marked loss stability the gravel when slightly above the optimum 
water content. 

For all permeability tests, samples containing only those particles passing 
3-in. were used. Mechanical analyses this material are shown 
Fig. curves and The range void ratios which the material was 
tested was limited those expected during ordinary construction the fill. 


EQUIPMENT AND METHODS 


13.54 in. diameter and in. high ft) were 
used throughout the testing program for compaction and permeability. The 
permeameters were arranged with screens top and bottom permit the free 
passage water and from the sample, and they were operated either 
variable-head constant-head upflow permeameters, shown diagrammati- 
Except for very pervious materials, all tests were run with the 
variable-head hookup. 


0.5 1.0 1.5 2.0 


Scale Feet 


> 
c = 
2 
2 
> 


Screens, Top and Bottom 


Overflow Trough 


Constant Head Tank 


Water Supply 


Valves 


No. Mesh and 


Compaction.—The material was compacted the compaction cylinder 
permeameter three layers with seventy-five 12-in. free-fall blows per layer 
25-lb hammer with 3-in. diameter face. This gave compaction which 
was equivalent that obtained when the material was compacted the 
embankment. After each layer was compacted, the surface was completely 
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break any surface skin that might tend the 
permeability. 

water was distilled and preheated above room temperature 
that there was falling temperature gradient through the sample during 
the entire run. This assisted the removal trapped air the sample and 
hence the prevention air-locking. 

apparatus for de-airing water was assembled similar that described 
elsewhere was impossible obtain enough water run 
all the tests and, since preheating the water seemed effective this 
material, the de-airing was discontinued entirely. 

Test sample was compacted the cylinder the desired 
moisture content and loaded with ton per ft. Water was started flowing 
through the sample with constant-head hookup shown When 
the sample had become saturated, rather when water appeared the upper 
surface the sample, the valves were turned for variable head, and duplicate 
permeability tests were run. For the long-time tests, the valves were turned 
back for constant head, and each thereafter, for period from one 
eight weeks, permeability determinations were made. 


AND RESULTS 


The specified moisture limits the embankment were within above and 
below the optimum moisture. allow for factors uncertainty field 
control and determine any dangerous difficult conditions that might 
approached the upper and lower moisture limits, laboratory tests were run 
the material ranging from below the optimum water content 
above. 

Tests were run the materials four groups, follows: Gravel with 
(1) 20% till admixture, (2) 25% till admixture, (3) approximately 40% till 
admixture from test fill, and (4) core material with approximately 35% till 
admixture, sampled during construction the dam. 

Group 20% Till cubic yard the core material containing 
20% till and 80% gravel was thoroughly mixed, and successive 
quartering was divided into eighteen samples approximately 
each. Three the samples were used determine the optimum water 
content the material. the basis this determination, three the 
samples were placed the permeameter below the optimum, three 
below, three the optimum, three above, and three above 
the optimum. Fig. shows the results compaction tests plus the unit 
weights obtained when compacting the material the permeameters the 
various water contents. 

each sample was placed the permeameter, was saturated and the 
permeability test run. All samples were placed nearly possible the 
same manner except that the water content was varied. There was, course, 
some unavoidable variation the void ratio. 


Experimental Investigation Protective Filters,” Bertram, publieations from the 
School Engineering, Harvard Univ., Cambridge, Soil Mechanics Series No. January, 


~ 


. 


September, 1942 MUD MOUNTAIN DAM 1093 


When permeability plotted against void ratio, semilogarithmic paper, 
straight-line relationship usually apparent for most soils. group 
however, there none the usual straight-line correlation between the 
permeability and the void ratio the entire series, shown Fig. 4(a). 
However, there good correlation between each set three points obtained 
from those samples placed the same water content. Fig. shows 
definite correlation between the permeability and the water content which 
the material was placed. 


Unit Weights Obtained When 
Compacting Samples Into 
Permeameters Various 

Moisture Contents 


Unit Dry Weight, in Lb per Cu Ft 


Percentage Moisture Portion 


Group 25% Till similar series tests was run for 25% 
till admixture with the same permeability-moisture correlation obtained 
group These tests are plotted, permeability versus water content, 
Fig. ratios are plotted Fig. and show none 
the usual straight-line relationships. 

Group 40% Till third series tests was run using samples 
selected random from test fill constructed core material with approxi- 
mately 40% till admixture. The permeability-moisture correlation Fig. 5(c) 
better than might expected when considered that this material varied 
considerably till content and mechanical analysis. Permeability-void 
ratios are plotted Fig. 4(c) with the same complete absence correlation. 

Group Embankment shows the results observations 
256 the samples taken for permeability control the core the dam during 
the early part construction, and directly from the borrow pit during excava- 
tion operations. The till admixture varied from 25% 40%. tests were 
run the water content which the samples were received the laboratory. 
Here again, the permeability-moisture correlation was surprisingly favorable 
when considered that the mechanical analyses the samples varied con- 
siderably. 

will noted that, both groups and there limited range 
water contents. Had the range been wider, there little doubt that the curves 
could drawn shown the broken-line extension the curve Fig. 
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Time samples that are placed the permeameters water 
content below the optimum, the permeability high. When water run 
through the sample for days, the permeability may drop certain cases 
the reduction generally occurs during the first few minutes after water 


0.0001 0.001 0.01 0.1 1.0 10. 
Permeability, Centimeters per Secx 


PERMEAMETER 


introduced into the sample. None the individual samples have shown the 
extreme variations found samples placed above and below the optimum. 
samples placed above the optimum, the permeability remains low in- 
definitely, except for occasional slight increase that may attributed 
removal trapped air, which increases the effective voids. 


INTERPRETATIONS 


When test results are inconsistent with the usual approximation 
permeability-void ratio relationship, the conclusion inevitable (1) that there 
are large experimental errors the tests (2) that there are factors other 
than the numerical void ratio that may influence permeability. 

After eliminating all experimental errors far possible without reducing 
the anomalies, becomes evident that the void ratio expressed number 
does not define the nature the voids. Numerically, the void ratio measure 
only the total volume voids. This total void volume may distributed 
infinite variety sizes and shapes voids, depending the distribution 
soil particles. 

granular soil uniform particle size, the permeability-void ratio 
correlation fairly consistent. This explained postulating that the size 
and distribution the voids are more less constant. well-graded soil 
with small amount clay included, consistent permeability-void ratio 
correlation can obtained only when, and if, the particle and void distribution 
can controlled. Control may exercised under two conditions: 
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(1) any given sample, the void ratio may reduced progressively 
consolidating the sample without appreciably disturbing the particle and void 
distribution. permeability tests are run concurrently with the consolidation 
and the permeability-void ratio relationship plotted semilogarithmic 
paper, the plot will approximate straight line. 

(2) series samples the soil, the particle and void 
distribution appears controlled exercising control the water content 
which the sample compacted into the permeameter. constant water 
content but with varying void ratio, the results permeability tests when 
plotted semilogarithmic paper again approximate straight-line perme- 
ability-void ratio relationship. 


The optimum water content sample Mud Mountain core 
material, with maximum particle size in. determined compaction 
investigation, 11.5% the dry weight the entire sample. the same 
sample, the water content the fractions passing screen 16.5%. 
Now, the optimum small sample using only those fractions passing 
screen determined, found 15.0%. Thus, when the 1-cu-ft sample 
containing the 3-in. maximum particle size its optimum water content, 
small aggregations those fractions passing }-in. screen are 1.5% above 
their optimum. With the fractions passing }-in. screen slightly above their 
optimum, the clay, and particularly colloidal fractions, are very plastic 
condition and are relatively free uniformly distributed throughout the 
soil system during the compaction process. The montmorillonite clay becomes 
very unstable the presence water and lends itself readily such uniform 
distribution. Any increase the water content which the sample com- 
pacted above this stage reflected general increase all voids the soil, 
and the permeability-void ratio relationship then approaches straight line 
(see Fig. 4), reflecting the numerical void ratio. 

this same sample were compacted into cylinder water content 
1.5% below the optimum, the aggregations finer particles will their 
optimum, will well compacted, and, therefore, will have relatively low void 
ratio. Surface tension and cohesion acting within the clay aggregations will 
contributing factors maintaining the aggregations their maximum 
density, thereby enhancing nonuniform particle and void distribution 
throughout the soil mass. The entire sample, being below its optimum, will 
have relatively high void ratio because cannot compacted densely 
when the optimum. Any decrease the moisture content which the 
sample compacted below this stage reflected further decrease the 
density and increase certain the voids the soil, and the permeability- 
void ratio relationship approaches straight line (see Fig. 4), again reflecting 
the numerical void ratio. 

The numerical void ratio appears govern the permeability the upper 
and lower limits water content which the materials were tested, but, 
relatively narrow range water contents near the optimum, enormous 
changes permeability take place with change void ratio necessary. 
This range water content marks the transition from condition more 
less uniform distribution the fines the wet side condition concen- 
tration fines into aggregations held together cohesion and surface tension 
the dry side. interesting note that this range moisture either 
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side the optimum almost the exact range originally specified for field 
control Mud Mountain Dam. 


APPLICATION 


determining the permeability certain pervious foundation materials, 
sometimes quite impossible obtain undisturbed samples test the 
material place, which case common practice determine the perme- 
ability remolded samples over range void ratios. this done con- 
junction with single consolidation test, has been suggested that the 
permeability-void ratio relationship should approach straight line. How- 
ever, series tests run samples the material compacted various 
water contents, the results will indicate the degree confidence that may 
placed the predictions seepage through the foundation materials. wide 
variations permeability are obtained that cannot attributed void- 
ratio change, some other means determining the permeability the founda- 
tion materials indicated. 

The early practice puddling cores earth dams accomplished the same 
end attained compacting this material some moisture content slightly 
above the optimum. Puddling results uniform soil particle distribution, 
uniform void distribution, and hence low permeability for that par- 
ticular soil. The fact that, the Mud Mountain core material, there only 
1.5% clay, and puddled core there may nearly 100% clay, 
does not alter the fact that high water content assists uniform particle 
distribution. 

INVESTIGATIONS 

The investigation has been limited thus far single material that 
rather distinctive its particle-size distribution. this material, the effects 
moisture are extremely pronounced that they completely mask, over 
certain range, all effects numerical void-ratio change. Investigations 
materials from four other sources, with narrower range particle sizes and 
montmorillonite clay, are being conducted. indicated that the results 
will not reach the almost impossible extremes the results obtained using 
Mud Mountain borrow-pit material, but that herein may lie the answer 
certain minor variations from the expected straight-line permeability versus 
void-ratio relationships. 
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OPEN CHANNEL FLOW HIGH VELOCITIES 


One the problems the engineering design chutes steep spillway 
channels the proper freeboard allowance made for the entrainment 
air high velocities. Calculations the channel dimensions also must 
modified with respect the application normal frictional resistance which 
altered the result air entrainment. all cases where water flowing 
steep chute allowed accelerate over sufficient distance, air mixed 
with the stream and the phenomenon results. The same 
admixture air with water can observed waterfalls natural stream 
channels where the velocity sufficiently high draw air into the water 
its agitation. 

Very few complete observations have been made for the measurement the 
volume air which will mixed with the water under these conditions. 
Model tests are futile sufficient velocity cannot attained create the 

phenomenon. gather information the behavior water flowing 
steep channels and measure the bulking the water, some observations 
actual spillway chutes were made 1936 and 1937. The Pacific Gas and 
Electric Company granted permission for measurements made three 
their chutes and cooperated maintaining uniform discharges during the 
period the tests. 

Incidental the measurement the bulking the water, was possible 
also observe the performance the stream while flowed around vertical 
and horizontal curves the spillway channels. The data gathered during 
these tests are presented this paper, together with theoretical analyses 
leading the development formula for computation air bulking and 

general conclusions its application. 


PROCEDURE 


Ordinary field methods cannot used measuring velocities ranging from 
per sec. The initial test was made Hat Creek No. chute 
western state. The spillway channel was laid out for observation July 11, 
1936, and measurements were taken the following day. obtain more 
complete data this chute second set observations was made April 17, 
1937. 


Hydr. Engr., East Bay Municipal Utility Dist., Oakland, Calif. 
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The second test was made Rapid Flume July and 21, 1936. The 
third test was made the South Canal chute another river project. The 
gages were located and levels run April 24, test was made 
May 

The general field procedure was first place gages along the chute channel 
for the purposes measuring the velocity and cross-sectional area. The 
distances were measured along the actual slope the chute rather than using 
the horizontal distances, stations being marked every scribing 
paint mark and the station number the side the flume. For the concrete 
chutes, gages were stenciled both inside surfaces the channel right 
angles the bottom slope, using prepared templet, the paint being quick- 
drying black lacquer applied with long-handled paint brush. Painted gages 
were not used the wooden channel Rapid Flume, but measuring point 
was established top the flume and distances were measured down the 
water surface. Levels were then run each station and the elevation 
the top the chute wall recorded. the same time the channel dimensions 
were checked. 

The procedure for the actual observations, following this preliminary work, 
consisted several steps: First, the discharge the power canals above the 
chutes was gaged before and after the test. Four discharge measurements were 
usually made from gaging bridge, using both the Price and Hoff current 
meters, two measurements being taken the beginning the test and two 
the end. Observations were also taken frequent intervals gage 
installed the power canal note any fluctuations the water level. Then, 
after water had been turned into the chute and the flow had become steady, 
readings were taken all gages along the channel. When the water travel- 
ing high velocity with entrained air, particles water leave the main 
body the stream with the result that the water level choppy and oscillates 
and down. Therefore, observations were taken the maximum and 
minimum levels each station over sufficient period time note the 
extreme range. mean these readings was then used the average water 
level. cross waves developed the chute, due curvature entrance 
conditions, additional readings were taken across the channel measuring 
down the water surface from board laid across the top the chute. 

Finally, velocity measurements were taken between various stations along 
the channel. Prior making the measurements coloring matter was mixed, 
using either flourescein potassium permanganate. Commercial flourescein 
was used, with lime first added the mixing water insure more complete 
solution the crystals. About one teaspoonful flourescein pint 
water was used for each 240 per sec discharge. For the potassium 
permanganate about ten teaspoonfuls pint water were used for each 
240 per sec discharge. The colored solutions were poured into half-pint 
pasteboard containers prior insertion the stream channel. Before the 
test was begun, the observers were placed their positions along the flume. 
Men with the dye solutions were stationed the head the flume the 
start the course timed. second observer with the revolver was 
located initial station, and the timers with stop watches were located 
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the lower station. Ata given signal, the container with the dye was dropped 
into the head the channel and the upstream observer noted the passing 
the color his station firing the revolver into the air. The timers the 

lower end the stretch started their watches this signal and stopped them 
when the center the slug color passed the lower station. Four six 
observers with stop watches were used time each color cloud, and the average 
their observations was used determining the mean time interval. 

Water Surface 


Lost Entrance 


Water Surface 368 +75 

Water Surface 115 
per Sec 


Bottom Channel 


2900 
> 
Hat Creek chute the Hat Creek No. power plant located 
short distance upstream from the penstock leading the power house. This 
steep straight concrete channel having vertical curves and dropping 189 air 
vertically distance 409.5 measured along the channel. The waste 
water from the flume was discharged through side-channel spillway weirs 


located both sides the canal. After dropping over the weir lip the dis- vel 
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charge was gathered into four collecting channels which unite the head 
the chute. 

The concrete the chute was rather rough, with pronounced form marks. 
The tie wires used constructing the forms had not been clipped off from the 
concrete and these added the roughness the surface. The channel had 
uniform width 5.75 with side-walls high except near the entrance and 
outlet where the height was increased. 


Due two the collecting channels entering the chute inlet right 
angle to-the direction flow, cross wave was developed creating rather 
turbulent entrance conditions and also increasing the percentage air the 
water this point. However, the cross wave and portion the entrained 
air had disappeared about The profile the chute shown Fig. 
The first 275 the chute consisted tangents connected short-radius 
vertical curves. Immediately below Station 2+75 there was short convex 
vertical curve leading steep pitch having inclination The 
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lower part the curve below Station 3+60 was concave vertical curve 
forming the outlet bucket. 

The flow the chute accelerated Station 2+75. this section the 
channel running waves formed along the side-wall and gave the appearance 
higher water levels than the center the channel. However, check made 
measuring the depth water the center portion the stream showed 
that the maximum height reached the water was the same that observed 
the side-walls. The upper surface the stream was very turbulent, with 
balls water that would dart from the main body the water surface 
high velocities drop back again short distance into the general mass 
aerated water. The entire surface the flow was composed these small 
waves rollers which ranged from 0.4 0.6 height above the main 
body the stream. 

Station 2+75 the curvature was too sharp followed the trajec- 
tory the flow, with the result that the upper layers water left the bottom 
the channel. the maximum flows some the spray overtopped the side- 
walls this steep section the chute. The cavitation caused some erosion 
the concrete. The flow returned the bottom the chute between Stations 
3+25 and 3+75, depending the volume flow. The velocity the water 
was appreciably reduced this section the channel, where the stream was 
spread through larger cross-sectional area. result, the water left the 
bucket the lower end the chute with lower velocity than that observed 
Station 2+75. The jet from the bucket discharged into the stream channel 
Hat Creek above the normal stream water level, the energy being dissipated 
eddies and turbulence the pool formed the channel the discharge. 
The general layout the chute and the flow conditions during the test are 
shown Figs. and view the flow conditions with discharge 
395 per sec presented Fig.2. this station the velocity accelerates 
from per sec persec. 2+70 water apparently slightly 
higher along the sides the chute than the center due wave action caused 
retardation along the side-walls. The velocity this point per sec. 


FLUME 


This chute located just above the penstock leading the high-head 
plant power house. The flume used discharge water from the Upper 
Standard Canal down the slope into natural stream channel leading 
Lake Tabeaud, where used afterward through the low-head 
wooden chute having horizontal curve to-the right near the upper end 
concave vertical curve near the terminal end the channel. The total vertical 
drop 145 distance 667 measured along the chute. The discharge 
water from the Upper Standard Canal into the flume controlled 
head gate. Water passes from the canal through culvert pipe and drops 
into the tapered inlet the flume, creating considerable initial turbulence 
the water. addition this tapered section the chute channel built 
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horizontal curve the right Station steep section the chute, 
starting Station (Fig. 5), gradually merges into long concave 
curve variable radius extending from Station 2+75 Below 
Station 5+25 the flume built flatter grade down the baffle box used 
check the velocity. The water discharges from the baffle box into the 
natural drainage channel described previously. 


The channel constructed wood with longitudinal battens, and below the 
taper section has average width 4.6 ft. The average height the side- 
walls 2.4 ft. Cross waves develop the upper section the flume, due 
the horizontal curve. These disappear the velocity accelerated the 
steep straight section. The air entrained the inlet the flume due 
entrance conditions partly lost and the entrained air again drawn into the 
water, the insufflation commencing along the side-walls the flume until 
the entire cross section below Station 1+50 has become aerated. Where 
turbulent entrance conditions not exist the entrainment air commences 
along the side-wall, described 1936 Lane,? Am. Soc. 

The velocity Rapid Flume accelerated maximum per sec 
Station 3+75 but gradually diminished the grade flattened until was re- 
duced velocity per sec the water passed into the baffle box. 
The profile the flume shown Fig. 
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CHUTE 


This chute located the end the South Fork canal which picks 
water below power house. planned ultimately use this water 
through Power Plant No. but the present time the flow only used for 
irrigation. The water discharges down the chute into natural channel. 
This channel 3,535.25 long measured along the slope and has total ver- 
tical drop 140 ft. The discharge from the end the South Canal con- 
trolled radial gate set across the channel. Water can also discharge 
side channel spillway which drops the water back into the main chute 
channel within distance 100 below the radial gate. Water dropping 
into the chute from the side channel augmented the percentage entrained 
air the head the chute. water had entered this route, much 
greater distance channel would traveled the flow before the entire 
cross section was aerated. 

The channel the chute was constructed concrete, the surface being 
rather rough and pitted, probably due sand being carried down the chute 
high velocities. the bottom particularly, the cement and sand had been 
worn away, leaving the aggregate fragments exposed. From the head the 
channel Station 23+00 the width the cross section was ft. From 
Station 23+00 Station 23+50 there was tapered section and from the 
latter station the end the chute the channel had width ft. The 
side-walls the 10-ft section varied height from 3.5 4.5 whereas 
the 7-ft section the height the side-walls varied from 3.5 5.25 ft. 

There are thirteen horizontal curves radii varying from 100 300 ft, 
the sharpest having internal angle 98° and radius 142 ft. The floor 
the channel these curves not superelevated the outer edge. 
opportunity was afforded here observe the flow water around curves. 
The total angular deflection the entire chute was 422° stated, the 
entire length channel was 3,535.25 ft, which 1,328.8 was located 
curves. There are numerous changes grade, including three re-entrant 
grades inserted decelerate the water the toe steep sections. Except 
for steep pitch the inlet, the steepest bottom slope 0.100 for about 300 ft. 
The other principal grades are 700 0.055, 600 0.030, 200 0.022, 
600 0.011, 715 0.025, and 170 three short lengths average 
slope 0.080. 

total time min sec was required for flow 334 per sec 
pass from the head the chute the outlet, the average velocity being about 
per sec. Due the many changes grade, the velocity various sec- 
tions the channel varied from per sec. The profile the channel 
not shown, due its great length. 

the channel bottom the curves was not superelevated, the water piled 
the outside each curve and produced standing cross waves the tan- 
gent section below the curve (see Fig. 7), addition the choppy roller and 
surface waves elsewhere noted normal chute flow. These cross waves lead 
error computing the wetted area some sections only side gage 
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readings were taken. Such errors are not sufficient magnitude affect 
the results for the entire chute. the 10-ft width channel the supereleva- 
tion the water surface the outer side the curve did not form uniform 
slope across the entire cross section, due larger percentage the flow 
following the outer wall. However, the 7-ft channel widths the slope 
the water surface was more less uniform that the difference level between 
the gage readings the outer and inner walls the curves indicated the 
slope the water surface across the stream section. the curve 142-ft 
radius, the 7-ft section, constructed 10% grade, the flow was very 
turbulent (see Fig. 8). The discharge this section tended proceed 


2600 


1+25 


Bottom 


Elevation, Feet 


series chords with result that water crowded certain points the outside 
the curve, permitting spray discharge over the top the wall. Between 
these points the water surface was well below the top the wall. 


NOTATION 


The letter symbols this paper are defined where they first appear and are 
assembled for convenience reference the list symbols appended the 
Symposium. 


THEORY FLOW CHUTES 


The objective, desired gained the tests, was the determination 
formula for the calculation air entrainment projected spillway channels 
large capacity. practical application, the method procedure 
should satisfy three requirements: (1) should conform nearly practicable 
established methods hydraulic calculations; (2) should sufficient 
simplicity readily applied; and (3) should yield results that are 
reasonable conformity with the measured tests air entrainment. 


Fic. VELOcITY, LONGITUDINAL WATER SURFACE 
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Meager data available indicated that with small concrete chutes, having 
velocities ranging from per sec, the entrained air varied from 15% 
35% the cross Measurements masonry lined chute 
Austria gave air entrainment ranging from 60% The use such 
large percentages for air bulking design would make large spillway chutes 
for the discharge river floods very costly. 

was indicated that the velocity, the distance traveled, and the roughness 
the channel had some bearing the air entrainment. Personal observation 
had led the conclusion that air bulking was reduced with larger cross sections 
and greater volumes discharge. 
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developing the theory flow chutes the following assumptions were 


made: 
The value the roughness coefficient the Manning Kutter for- 
mula constant for the particular type material which the chute con- 
structed; 
The air and above the water caused additional loss energy, the 
the specific gravity compensating for the added area; 
computation, computed hydraulic radius can used 
with smaller value than normal; and 


Channels from Surge Oberbauret Dr. Ing. Karl Innerebner, World Power 
Conference, 1924, Vol. II, 481. 
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The velocity head computed from the mean velocity can used 
without substantial error. 


Eq. the discharge cubic feet per second; the observed 
average velocity section, feet per second; and the computed wetted 

The true velocity head greater than that computed from the mean 
velocity, since the true kinetic energy the summation the product the 


velocity head and the weight water for elemental areas, the cross section. 
Then 


and 


channel; and acceleration due gravity. The value Eq. 
represents the average energy head per pound water flowing per second; 
always greater than unity and ranges from 1.02 1.20. Since the velocity 
distribution the chutes tested this investigation was not obtained, was 
assumed that the corrections were small and could neglected. This assump- 
tion has the effect slightly increasing the energy loss above the actual values. 
the design large spillway chutes, the inclusion reeommended and 
value ranging from 1.07 1.10 customary. 

When air entrained water, the basic formula must modified 


which the ratio water mixture air and water the relative 
density the mixture ratio air). Assuming that the flow 
channel has reached uniform velocity and neglecting internal friction the 
water, the relation between the force represented the weight the water 
and the frictional resistance (f) the channel walls for cross-sectional area 


which the length measured along the channel; the weight 
water 62.5 lb); the slope the channel degrees; and the 
wetted perimeter for mixture air and water. 

When air entrained the water, the normal assumption given most 
hydraulic textbooks the relation between the friction and the velocity 
for uniform flow must modified the entrained air reduces the density 
the mixture and the friction 


Hydro- and Aero-Mechanics,” Prandtl and Tietjens, Engineering Societies 
Monographs, McGraw-Hill Book Co., Inc., 1934, 71. 
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The observed hydraulic radius computed from the observed area and 
wetted perimeter the mixture air and water must used. Then, the 
normal Chézy formula will apply. 
One difficulty applying this formula that the ratio air the water 
must known advance. This ratio not constant but varies with 
the velocity, size channel, and roughness the side-walls. Hence, 
find the retardation factors for nonuniform flow with the data available 
from the chute tests, the values Kutter’s could most easily computed 
from the slope the energy gradient using the Manning formula 


2.2082 


desired, the corresponding values the Kutter formula then can 
determined the relation between Manning’s and Kutter’s for various 
values 

the nominal value. The relation between the values for the Manning 


formula would 


therefore, were also calculated. Some writers have recommended the use 
the observed area with smaller velocity, This procedure 
would necessitate the use larger value, which likewise has been com- 
puted from the available data. This method has the advantage utilizing 
the actual area, but there are several objections the use fictitious velocity 
less than the actual which will pointed out later. The most logical procedure 
would appear the calculation the actual velocity from the calculated 
depth without air and the determination the true depth air plus water 
from the percentage entrained air. The relation between the entrained 
and the velocity, hydraulic radius, and channel roughness from 
experimental data. 

Using Eq. the computed values are about the same 
channels similar material with velocities below critical. using the 
observed velocity the corresponding slope the energy and 
the computed hydraulic radius smaller value can computed desig- 
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nated here the other hand, using the observed hydraulic radius 
the computed velocity and the corresponding slope the energy gradient 


1.486 


The increment velocity flow progresses down the chute may ex- 
pressed function the time. However, this leads equation involving 
hyperbolic functions that are awkward 

simpler solution can derived expressing the increment velocity 
function the distance traveled. Such equation will developed 
applying rectangular section uniform width, Since many chutes 
are this type, this formula will have wide application. 


Applying Bernoulli’s theorem Fig. 


being the velocity head coefficient from Eq. 2b. 
Taking the derivative with respect 


rectangular channel, substituting the computed depth for the 
value the discharge per foot width becomes y-; but 


Flow Water Channels Under Steep Durand, Transactions, A.S.M.E., 
1940, (presented Semi-Annual Meeting A.S.M.E., San Francisco, Calif., July 10-15, 1939). 
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Substituting Eq. Eq. 11, 


but and substituting the value derived from Eq. 9a, 


The negative sign used Eq. since decreases with increase 


Equating Eq. and Eq. and solving for 


2.2082 


letting 


Then 


This the basic equation for flow steep rectangular chutes derived 
from the simplified relation rather than the more general equation 


the latter used, three variables are introduced making the solution 
more difficult. Then, 


solution cannot made unless f(V). Then fi(V) dV. 


| 
a V q 
(15) 
? V 
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Substitution Eq. Eq. gives 


The use Eq. would difficult. 
Returning the basic Eq. 17, the solution for using constant value 

determine exact value for chute computations, several assump- 
tions may may assumed that the slope the energy gradient 
may separated into two components representing frictional resistance and 
air resistance using the computed hydraulic radius This solution has 
fictitious character since the air resistance based the area occupied pre- 
sumably the water and not the air. was, therefore, discarded. 

Another solution could made assuming that constant for given 
value relation steep chutes between the slope and the entrained air 
varying with the hydraulic radius was determined Ehrenberger 
Then Eq. 19a becomes 


this ease, the actual radius should used Eq. Then 
similar steps used deriving Eq. 16b, 


(23) 


This solution best approximation since not constant for given 
value but varies with the velocity. many steep chutes the velocity 
accelerates throughout the length the channel. 

The most satisfactory solution was obtained using the Froude number. 
The air entrained the water would appear bear some relation the 
Froude number, air entrainment essentially related the kinetic energy 
the stream. The proper application the dimensionless parameter 
hydraulic phenomena has previously been pointed out® Boris Bakhmeteff, 
Am. E., and Arthur Matzke, Assoc. Am. Soc. 

the present case, the frictional resistance, addition the kinetic 
energy, has definite bearing the air entrainment. not incon- 
gruous use the radius the longitudinal parameter. 

After some experimenting with the data, was found that the best relation 
with the Froude number existed using the ratio the two ratios, air 


Die Wasserwirtschaft, Vol. 23, 1930, pp. 573 and 595. 


Hydraulic Jump Terms Dynamic Similarity,” Boris Bakhmeteff and Arthur 
Matzke, Transactions, Am. Soc. Vol. 101 (1936), 635. 
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Experimental data were plotted, will shown subsequently the paper, 


ro 


using ordinates and abscissa, and the following general rela- 


tion was found exist: 


which coefficient and constant for air entrainment. With 
much the data was found that but with smooth entrance condi- 
tions has finite value. transposition can determined, assuming 


Dividing Eq. Eq. 


Substituting the value from Eq. and replacing its 


value from Eq. 


246 and constitute basic formulas for the solution problems 
air entrainment spillway chutes. The use values obtained from 
Eq. 26, Eq. would appear afford the most rational theoretical solution 
for flow rectangular chutes. 

the design chutes having trapezoidal section, fact any shape 
other than rectangular, and transition sections any form which the 
width varies, Eq. cannot applied. Similar formulas can developed, 
but they are too cumbersome for practical use. easier solution obtained 
trial-and-error method. 

Starting with Bernoulli’s theorem from Eq. 10, the loss energy between 
two adjacent stations can expressed using subscripts and indicate 
upstream and downstream sections, respectively. Equating this the slope 
the energy gradient the Manning formula, the following equation results: 


1 p y2 U 
or, substituting for its value from Eq. 16b, 
i. p V B oF 
2/3 
72 
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using Eq. 27, value first assumed and the corresponding value 
the hydraulic radii the two sections. the two sides the equation 
not balance with the assumed value another value must tried. 
Usually after two incorrect assumptions, the true value can arrived 
interpolation. Eq. can used for rectangular sections, desired. Calcu- 
lations are facilitated first preparing series curves the several functions 
each section for varying depths. These must prepared for each discharge 
can determined readily. The ratio water (p) can calculated from 
Eq. 24a with the values and for specific depths using the values 
and assumed for the design. From the ratio water the section and 
the computed area, the actual area the cross section can computed: 


For trapezoidal sections, the actual depth air and water can found 
the formula 


(28b) 


which bottom width and side slopes trapezoidal sections. 
Values Ro, and then can readily determined. For use Eq. 
values are calculated from Eq. For practical use, curve the func- 
tion can plotted against the depth The curve can marked off 
indicate increments the depth for changes the nearest 0.0005. 
These values are sufficiently accurate and avoid too frequent change the 
value thereby facilitating the computations. 

Friction loss computations the Manning can made most 
readily first preparing curve the function 1.486 plotted against 
the depth Once value assumed the friction loss can computed 
the left-hand member Eq. from the values and 1.486 
interpolated from the curves. 

The computations flow down the incline must started from the head 
the chute. With smooth entrance conditions, the computations can start 
from the critical section. When the velocity factor used, the formula for 


critical depth usually given hydraulic texts must modified the intro- 


which water surface width trapezoidal sections. 

advantage determining the calculated areas solid water realized 
the design jump basin, which many cases forms the chute terminal. 


Compare with formula for critical depth Hydraulics Open Boris Bakhmeteff, 
ngineering Societies Monographs, McGraw-Hill Book Co., Inc., 1932, 35. 
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Here, the values and together with the actual water velocity 
the toe the chute should used the energy and momentum equations 
the hydraulic jump. the velocity factor has been used the energy 
equation, similar factor based the summation over the cross section 


High Low High Low High Low High Low 


(1) (2) (3) (4) (5) (1) (2) (3) (4) 
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0.8 
Velocity (c) 368 Sec (d) 395 Sec 
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Between these stations the flow did not follow the bottom the flume because the vertical curve 
was sharper than the trajectory the flowing water. 


must applied the momentum equation. This factor smaller than 
and values ranging from 1.04 1.07 are generally applied. (The correction 
the average velocity applicable the energy and momentum equations 
has been determined Coriolis, Saint Venant, and 


See Annales des Ponts Chaussées, 1836; also Transactions, A.S.M.E., Vol. 54, 1932, 57. 
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Calculation the jump basin dimensions this manner will avoid the 
inadequate lengths and depths encountered many constructed designs. 

Another the problems the design chutes the determination the 
path the water around vertical and horizontal curves. Concave vertical 
curves present difficulty design, except for stress created the centrifugal 
force. Convex vertical curves must not made sharper than the trajectory 
that would followed the high velocity water under the action gravity. 
the flume bottom steeper than this trajectory, the water will leave the 
bottom and may rise above the top the side-walls observed Hat Creek 
chute. The trajectory will parabola tangent the slope the chute. 


Dare, 12, 1936 Test Dare, 17, 1937 


Full Run Full Run Partial Run 
Stations Stations Stations 
Ratio, 
Col. 
Col. 


(3) (4) (6) (7) (8) (9) (10) (12) (13) 


2+75 5.135 4.95 .... 


0+25 3.6 1+25 


Observed times have been corrected factors determined for length travel. 
Equivalent time full run. 


The equation the curve, neglecting air resistance, adaptation the 
usual equations the path projectile given many engineering texts.” 

have acquired velocity equal the vertical component, the water 
would have dropped vertical distance 


this same time the flow would cover horizontal distance 


See American Civil Engineers’ Handbook, 5th Ed., 119. 


0+25 | 1+25 2.48 2.48 2+25 1.91 2.02 95 
0+25 | 2+25 4.65 4.52 2+75 2.70 2.865 94 : 
0+25 | 3+75 7.60 7.07 4+09.5) 2.065 | 2.695 77 
0+25 | 1425 2.29 2.29 2+25 1.46 1.51 97 
0+25 | 2425 3.89 3.79 2+75 2.24 2.58 87 
0+25 | 3+75 6.05 5.65 4+09.5) 1.74 1.96 89 = 
0+25 | 4+09.5) 6.68 6.17 
sin? 
(30a) 
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The point tangency the parabolic curve coordinates from 
the point origin. The equation the curve which 


many constructed chutes, found that the water does not follow the 
bottom the channel convex vertical curves. This condition probably 
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Velocity head, 
due the fact that the method calculation did not yield the true velocity 
the stream. Even the true mean velocity has been determined the 
method previously outlined, necessary determine the maximum velocity 
the channel which will occur near the surface the center portion the 
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cross section. short chutes this maximum velocity will practically 
equal the theoretical velocity resulting from the existing vertical drop. 
longer chutes, the maximum velocity will from 15% 20% greater than 
the mean Professor Lane found that the maximum velocity the 
Uncompahgre chute was 15% greater than the 

Aside from the possibility damage the concrete the chute from 
cavitation, from the tendency lift the entire bottom slab, there appears 
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serious objection allowing the stream leave the floor the chute 
convex vertical curves, provided the side-walls are sufficient height prevent 


Engineering News-Record, January 1936, 
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overtopping the spray. The observations the Hat Creek chute indicate 
that the spreading the jet over greater cross section increased the friction 
and aided the dissipation the energy. This fact may useful design, 
where desirable reduce the velocity the stream before discharging the 
jet from the end the chute. 

Since water traveling the chutes always less than the critical depth, 
changes horizontal direction are more difficult make and, possible, 
should avoided. horizontal curve becomes necessary, should prefer- 
ably inserted near the head the chute before the water has gathered much 
velocity. Model tests should made test the design prior actual con- 
struction. 


Run TIME Run Run TIME 


5.90 
3.72 


Observed times have been corrected factors determined for length travel. 


The water surface across the channel chute constructed curve will 
tangent course, this assuming that the channel 
relatively narrow permitting the use the radius the center line the 
average curvature the section. Then, the inclination this tangent with 
the horizontal 


Verification this formula was possible the South Canal chute, where 
the difference elevation the water surface between the inner and outer 
walls approximated the slope model conducted the East Bay 
Municipal Utility District, Robert Kennedy, Am. Soc. E., has demon- 
strated that the flow water around horizontal curve greatly improved 
the bottom the flume superelevated the outer wall accordance with 
this formula. 


The results the depth measurements and velocity observations are given 
Tables 11. The gage readings are recorded for both the highest and 


Text Book Hydraulics, Russell, Henry Holt Co., 4th Ed., 49. 
Western Construction News, April, 1942, 158. 
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lowest point observed each station, together with the mean depth calculated 
from these observations. the case Rapid Flume, readings are given for 
depth the center the flume for the upper part the channel where cross 
waves developed due the curvature the flume. The results for Hat Creek 
chute are given Table for Rapid 
Flume Table and for the South 
Canal chute Table 

making velocity tests, was 
better time the color over rather 
short stretch flume, because, the 
measuring course were too long, the 


TABLE 
Run 
CHUTE 


color cloud spread out over such long Upper Lower Corrected 
distance that was difficult estimate (1) (2) (3) (4) 
the time passing the center 

the color band. the first series 
observations color was added the 
head the chute and the time travel 
ing length having the timers moved 
indicated that the observed time the 23+00 36.36 33.39 
lower stations was too great. 25+00 7.65 7.65 
latter tests the color was dropped 5.93 5.93 
various places along the chute and the 
velocity was timed over courses not ex- 
would permit. the South 14.07 13.72 
chute the great length channel and 12.47 12.22 


without error. The results obtained 

from separate courses these tests 

were totaled give the entire time for the water travel from the inlet 
the outlet the chute. all cases the sum the times taken over the 
short courses was less than the observed time over the total course. This 
anomalous condition can explained readily. 

The color cloud was retarded along the side-walls proceeded down the 
chute and became attenuated with increasing distances travel. observing 
the velocities, the time was taken the approximate center the color cloud 
closely could judged eye. Along 400-ft course the color cloud 
extended between and Analyses the results indicate that probably 
the greatest concentration color was the head the slug; hence, the time 
lapse recorded the center the cloud was too long. shorter courses the 
color cloud stayed together better, and less error, therefore, was introduced 
estimating the location the center eye. The results these observations 
showing the variation time over long and short courses have been tabulated 
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for Hat Creek chute Table Rapid Flume Table and for the South 

Canal chute Tables and 7A. These observations short courses were 

used obtain factors applied correct inaccuracies the observed time over 

the longer courses. Since for purposes the depth readings needed 

recorded for the perpendic- 

TABLE distance, corrections were 

Run Stations, made multiplying the ob- 

CANAL SHUTE served depth taken perpendic- 

chute the cosine the angle 

full run the velocity, particularly near 

the head the chutes, required 

segregation the time over 

shorter courses than was prac- 

tical for timing. Therefore, 

+00 time from the observations for 

putation was made for 25-ft in- 

tervals Hat Creek chute and 

Rapid Flume and for 50-ft intervals the South Canal chute. making 

this subdivision, the times the various tests were compared graphically 

and increments for the 25-ft 50-ft intervals estimated from the elapsed time 
over the longer courses. 

comparison the recorded observations depth and velocity illustrates 
forcibly the difficulty obtaining precise results measuring flows high 
velocities where air entrained with the water. 

The velocity head for the adjusted mean velocity each 25-ft 50-ft 
section was plotted the profile the respective chutes (see Figs. and 
for Hat Creek chute and Rapid Flume), and the energy gradient was drawn. 
From this gradient line the velocity head each station was scaled and the 
corresponding velocity calculated. These velocities were used conjunction 
with the observed areas compute the coefficient retardation. addition, 
the observed velocity and energy gradient were used conjunction with the 


14+00 


area computed the relation and the corresponding hydraulic radii 


calculate retardation factors. The velocity each station was also com- 


puted the equation using for the area the values computed from 


the observed depths, including water plus entrained air. These velocity heads, 
which were smaller than the observed values, were plotted graphs and the 
corresponding energy slopes determined. Necessarily these slopes were greater 
the zone accelerated flow than those obtained use observed velocity 
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heads. These three methods computation were used Fred Scobey, 

his bulletin the computation the retardation factor has been made 
taking the energy gradient parallel the bottom slope. This may correct 
many instances, but the tests described this paper was found that 
the energy slope was seldom parallel the bottom slope, can seen 

These elements the water prism obtained the manner just described 
were used for the determination the coefficients the flow the Chézy 
formula and the retardation factors Kutter’s and Manning’s formulas 
each the three methods described; namely, using (1) the observed ve- 
locities and the observed area, (2) the observed velocity and areas computed 


and (3) the observed areas and velocities computed 
440 


The hydraulic properties the water prism the three chutes tested 
together with the percentage air the water and the value the parameter 


were tabulated for each flow observed. These results are given Table 


Hat Creek chute, Table for Rapid Flume, and Table for South 
Canal chute. 

Bureau Reclamation Tests Kittitas Chute—In 1938 the Bureau 
Reclamation made tests air entrainment the wasteway from the Kittitas 
Main Canal located the Yakima Project western state. This chute 
discharges water from the main canal into the Yakima River. The flow 
water from the canal into the inlet the chute controlled radial gate. 
The total length the chute 1,301 measured the slope, and the total 
drop from the canal the bottom the bucket the Yakima River 340 ft. 
The first 753 the chute are 10° 12’ slope and the initial 318 are 
built taper section. the end this slope there convex vertical 
curve about length connecting with steep slope 33° 10’, having 
length 289 ft. This joins with concave vertical curve about length, 
joining almost horizontal tangent with slight downward slope 04’ 
for distance the end the chute. The jet the end the incline 
discharges into the Yakima River above the maximum high-water level. 

The concrete the chute smooth the side-walls, but judging from 
photographs, the bottom appears roughened slightly, probably due 
sand carried down the chute with the discharge. The width the channel 
decreases from the inlet the end the taper section. From 
this point uniform width carried for the remaining length the 
chute. The side-walls are high, except the inlet and near the outlet. 

Discharge measurements were taken with current meter the canal 
above the chute. Measurements the depth the water surface the 
incline were taken with point gage, reading the distance down from the top 
the side-walls. traverse the channel was made, taking observations 
the side-walls and every foot across the section. The first point gage station 


Flow Water Flumes,” Fred Scobey, Technical Bulletin No. 393, Dept. Agri- 
culture, pp. 89-91 


| 
y? 

| . 

7 
e 
h 
Nn. 
yn, 

m- 

ds, 
the 

q 


1128 AIR Papers 


was located distance 343 from the inlet, short distance below the 
end the taper section. Depth readings were taken approximately 100-ft 
intervals for distance 410 along the 10° 12’ slope, for 210 along the 
33° 10’ slope, and for along the 04’ slope. The depth and area each 
the ten observation stations were calculated from these measurements. 

Velocity observations were taken electrically, using salt conductor, 
balls moistened flour salt being dropped into the chute for each test run. 
The velocity was recorded oscillograph through electrodes installed 
the walls the chute. The oscillograph consisted essentially three moving 
coil galvanometers, optical system, and motor-driven camera. Two the 
galvanometers recorded the passage the salt cloud past the upper and lower 
stations. The third was connected 60-cycle alternating current and 
recorded time scale the film. 

This method has the advantage affording permanent record the 
test. Nevertheless, was difficult determine the center the gravity 
the salt cloud passed each station. Also, was noted the film that 
the flow, particularly the steep grade, occurred surges, with the result 


TABLE Properties, CANAL CHUTE 
AND PERCENTAGE AIR WATER 


Hydraulic 
Depths (Ft) Radius Area (Sq Ft) Percentages 
Station tio, served locity 
Yo Ye Ro Re Ao Ae (p) p 


(2) (3) (4) (5) (6) (7) (9) (10) (11) 


GO 


EEE EEL EN 


CORR 


SCHON 


CORR OOWOO 


e 


3 g R. 
(12) 
0+50 
1+50 
2+00 
2+50 
3+50 
4+50 
+00 
7+00 
8+50 
+00 
9+50 
11+00 
11+50 
12+50 


September, 1942 AIR ENTRAINMENT 1129 


TABLE 8—(Continued) 


AND PERCENTAGE AIR WATER 


Hydraulic 
Depths (Ft) Radius Area (Sq Ft) Percentages 


served puted served puted served puted ter 
Yo Ye Ro Re Ao Ae (p) p 


(2) (3) (4) (5) (6) (7) (9) (10) 


1.60 1.51 1.25 1.20 14.7 94.5] 0.06 22.7 8.0 13.3 
head, corresponding the observed velocity. Width changes from 10.0 7.0 
between Station and Station +50. 
0.4 
0.4 
9.8 that the velocity varied the same course between several tests. Hence 
doubtful this method yields much greater accuracy than the color 
18.9 timing method. 


4 
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Observations velocity were recorded for ten flows, ranging from 
per sec 1,005 per sec. 

cross wave developed the head the chute, due entrance conditions, 
and was reflected from side-wall side-wall, but this had practically dis- 
appeared before the first point gage station had been reached. the convex 
vertical curve the flow left the bottom the chute and caused spray over- 
top the side-walls. This was particularly true the larger flows. This spray 
prevented point gage readings the first station the steep slope for all but 


PROPERTIES AND PERCENTAGE AIR WATER 


AREA PERCENT- 
(Fr) (Sq Fr) AGES Ob- Ve- 


Yo Ue Ro Re Ao Ac (p) p) 


(a) 10° 12’ 


1-3 0.49 0.38 0.44 0.35 3.88 3.11 0.27 28.6 72.6 
3-5 0.50 0.41 0.45 0.37 3.96 3.25 0.22 27.4 63.1 

1-3 0.80 0.79 0.68 0.66 6.35 6.33 0.01 30.5 14.63) 44.3 
3-5 0.81 0.81 0.68 0.68 6.44 6.48 100 0.00 29.8 40.6 

5-7 0.83 0.79 0.70 0.66 6.55 6.33 0.03 30.5 14.46) 43.1 

5-9 0.82 0.76 0.69 0.64 6.48 6.11 0.06 31.6 48.5 

3-5 1.25 1.07 0.96 0.84 9.95 8.54 0.16 42.4 66.5 

5-8 1.24 0.96 0.89 9.85 9.10 0.09 39.8 55.3 

8-9 1.24 0.91 0.96 0.74 9.84 7.31 0.35 49.5 38.1 103.0 

1-3 1.26 1.26 0.98 0.98 10.06 10.66 100 0.00 37.6 44.8 
3-5 1.26 1.07 0.97 0.86 9.98 9.52 0.05 42.1 64.1 

1-3 1.54 1.48 1.14 1.075 12.37 11.66 0.06 42.1 51.3 
3-5 1.55 1.48 1.13 1.075 12.30 11.69 0.04 42.0 51.1 

5-8 1.56 1.38 1.14 1.02 12.36 11.01 0.12 44.6 60.6 

5-9 1.56 1.33 1.15 0.995 12.43 10.63 0.16 46.2 66.7 

3-5 1.80 1.61 1.26 1.145 14.33 13.98 0.10 45.6 32.33) 53.2 

5-7 1.79 1.55 1.26 14.28 12.41 0.15 47.3 62.7 

7-9 1.80 1.42 1.26 1.04 14.36 11.33 0.08 51.8 41.7 80.1 

1-3 2.02 1.93 1.36 1.30 16.11 1§.43 0.04 46.6 51.9 
3-5 2.03 1.95 1.37 1.31 16.12 15.63 0.03 46.0 50.2 

5-7 2.05 1.82 1.38 1.25 16.33 14.55 0.12 49.4 57.9 

7-9 2.07 1.67 1.39 1.175 16.49 13.36 0.23 53.8 76.6 

1-3 2.24 1.86 1.45 1.26 17.87 14.91 0.21 52.1 42.2 67.0 
3-5 2.22 1.90 1.45 1.29 17.68 15.24 0.15 51.0 40.4 62.6 

5-7 2.21 1.91 1.44 1.29 17.58 15.27 0.15 50.9 40.3 62.5 

7-9 2.26 1.89 1.46 1.28 17.97 15.12 0.19 51.4 41.1 64.2 

1-3 2.52 2.14 1.56 1.40 20.02 17.14 0.16 53.8 45.0 64.3 
3-5 2.49 2.14 1.55 1.40 19.79 17.14 0.15 53.8 45.0 64.3 

5-7 2.50 2.2 1.55 1.44 19.84 18.08 0.09 51.0 40.4 56.1 

7-9 2.51 2.13 1.56 1.39 19.96 17.04 0.16 54.1 45.5 65.5 

3-5 2.65 2.31 1.62 1.46 21.11 18.49 0.14 54.4 46.0 63.0 

5-7 2.66 2.25 1.62 1.44 21.15 18.01 0.16 55.8 48.4 67.2 

2.66 2.37 1.62 1.49 21.19 18.93 0.19 53.1 43.8 58.8 


q 
3 
q 
1 


(1) (2) (3) (4) (6) (7) (8) (9) (11) (12) 
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TABLE 9.—(Continued) 


AND PERCENTAGE AIR WATER 


AREA PERCENT- 
(Fr) (Fr) (Sq Fr) AGES Ob- 
served 
Col. ve- 
sec Ob- Com- Ob- Com- Ob- Com- Wa- Gol. locity 
served puted served puted served puted ter 
Ro Re Ao Ae (p) 


(3) (4) (5) (6) (7) (9) (10) 
(b) 33° 10’ 


8.07 
9.04 
9.04 


11.83 
11.83 
12.88 
12.88 


11.75 
12.08 


15.32 
15.32 
15.32 
15.02 
15.02 


16.16 
16.16 


17.42 
17.42 


19.00 
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(c) 


17-20 
18-20 


17-19 
17-20 


17-19 
17-20 


17-19 


17-20 
18-20 


17-20 
18-20 


1.12 0.98 12.03 10.49 34.5 18.50 


1.06 0.89 14.41 9.13 53.8 45.0 
1.06 1.01 14.41 10.81 45.4 32.05 


1.41 1.40 16.95 17.42 .00 33.7 17.67 
1.41 1.09 16.95 12.00 48.9 37.18 


1.50 1.11 18.69 12.50 57.5 51.4 
1.50 1.14 18.69 11.27 63.8 63.3 


1.58 1.24 20.76 14.36 54.1 45.5 


1.68 1.35 22.93 16.46 56.0 48.8 


1.68 1.45 22.98 18.17 55.3 47.6 65.7 


oo 


wt we we wt we 
oo & 


Velocity head, corresponding the observed velocity. Point gage readings made the end 
the section because spray. 


1131 
oc- 
nit’; v2 4 
ead* R. 
(1) (2) (11) (12) 
193 120.0 
222.0 
14-15 226.0 
362 11-14 167.0 
12-14 152.0 
14-15 148.0 
14-16 123.0 
401 11-13 160.0 
13-16 181.0 
491 12-13 130.0 
11-13 164.0 
12-14 150.0 
13-15 179.0 
14-16 147.0 
587 11-13 142.0 
13-15 161.0 
13-16 165.0 
719 11-13 185.0 
13-15 135.0 
13-16 158.0 
777 11-13 126.0 
13-16 296.0 
922 12-13 147.0 
12-14 134.0 
13-15 158.0 
2.68 124.0 
491 1.14 101.0 
4 . 
2.3 
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TABLE SURFACE CONVEX VERTICAL 
Hat 
(395 per Sec) 


Mean 
Mean Maximum Surface Maximum Computed 
tion (ft 

sec 


served served rected per sec per sec per sec 


433.09 
17.44 20.46 
80.23 


2,799.7 


Corrected for slope. Total drop Station 2+75 129.3 ft. The theoretical velocity equals 91.3 
per sec. 


the low flow; and even the second station readings were prevented the 
highest flows. 

From the data furnished the Bureau Reclamation, calculations 

the retardation factor have been made the three methods previously 

described from observations for 

TABLE 10°12’ slope and for the 

ANOTHER PROJECT tions were not taken the 

(Radius Curve, 142 Ft) vertical curves the taper 

section, which prevents the cal- 

profile for this chute. 

The hydraulic properties 
the water prism, together with 
percentage air the water 
and the values the parameter 


were tabulated for each 


the observations Table 


ANALYSES RESULTS 


computed tangent agrees reasonably well with the 
the P.C. the curve due cross waves resulting from tion factor conform, general, 
the curvature the channel upstream. 
with the range discussed under 
the heading, Flow 
Exceptions were noted near the inlets where high values were 


recorded result entrance conditions. Similar high values occurred 


| : Ne- 
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below the convex vertical curve the Hat Creek chute due the spread- 
ing the cross section the water left the bottom the chute. Also, 
the sharp horizontal curve the South Canal chute increased fric- 
tion loss was noted. each chute the most representative values both 
the Manning and Kutter formulas have been tabulated, leaving the excep- 
tions for separate comment. 

The retardation factors obtained from the observed areas and observed 
velocities are equivalent those observed normal channels. This value, 
therefore, has been designated The retardation factor obtained using 
the observed velocity and the area computed designated and 
the values obtained using the observed areas and the velocity computed 


has been designated The use observed velocities and observed 
flo 


areas gives results for the retardation factor comparable with those obtained 
wood and concrete channels under normal flow conditions contrasted with 
the shooting flow being tested. The values are less than normal values 
and correspond those which could computed Eq. 26. The value 
considerably larger than the results obtained the first two methods. 
The results obtained for all chutes (see Table 12) are the most in- 


CHUTE 


Formula Slope, 10° 12’ Slope, 33° 10’ 


Manning 
Kutter 


Kutter 0.008 0.010 0.008 0.011 0.012 0.014 0.015 0.015 


Manning 0.024 0.031 0.014 0.019 0.015 0.016 0.016 0.018 0.030 0.055 
Kutter 0.023 0.030 0.014 0.018 0.015 0.016 0.016 0.018 0.030 


consistent the three methods used, largely due the vagaries the computed 
energy gradient. Therefore, this method has little either 
The values for the Hat Creek chute include the tests for 115, 163, 368, and 
395 per sec observed the straight section the channel between 
Stations 0+50 and 2+75. The retardation factor materially increased the 
steep grade below the convex vertical curve and the first portion the 
bucket. Using the observed values velocity and area, the values Kutter’s 
and Manning’s average about 0.020 this section the channel. The 
values and correspondingly increase over values previously cited the 


0.010 0.013 0.010 0.012 0.012 0.013 0.015 0.016 0.019 0.021 
0.010 0.013 0.011 0.012 0.012 0.013 0.015 0.016 0.019 0.021 
a- 
al, 
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other methods computation. The increase the retardation factor this 
section the channel due the increased friction resulting from the path 
taken the trajectory the flow around the vertical curve. This trajectory 
has been calculated for 395 per sec means Eq. 31, using the 
calculations the mean velocity and also velocities 10% and 20% above the mean 
velocity. These results are given Table 10. 

The mean velocity the start the vertical curve was per sec, but 
comparison the results Table indicates that the water surface follows 
trajectory corresponding velocity per sec. This equal the 
theoretical velocity the water falling from the inlet the chute the start 
the vertical curve. short channel the water surface the center 
the chute, therefore, not retarded friction, but travels the theoretical 
velocity. compensation, the velocity the side and bottom the incline 
necessarily must travel less than the mean velocity, with the result that 
vertical curves the bottom layers should adhere the floor the chute, whereas 
the upper layers separate from the water prism, throwing the stream into 
spray. Air resistance ultimately slows the surface velocity comparison 
the calculations with the observations shows that the surface curve drops 
suddenly Station 3+25 after following the trajectory for velocity 
per sec for distance ft. The increase the retardation factor this 
section the chute indicates, previously mentioned, that there exists 
possibility dissipating energy this means. The same principle utilized 
some extent discharging water from chute bucket above the water sur- 
face the discharge pool. The throwing the jet into spray assists dis- 
sipating the energy, addition further dissipation energy eddies and 
turbulence the pool itself. 

The results Rapid Flume correspond large measure those obtained 
from the tests the Hat Creek chute. this case the convex vertical curve 
afforded uniform retardation factor throughout the length the chute after 
the turbulent conditions the inlet had been passed. 

the South Canal chute the computed retardation factors are more 
erratic. Some the apparent variation probably due faulty observa- 
tions area velocities. The average results for the entire length channel 
are reasonable, although when the roughness the bottom considered, 
value between 0.014 and 0.015 would have been expected, rather than the 
lower value actually obtained. 

the steep section the chute constructed the 142-ft radius, located 
between Station 23+00 and Station 25+00, the value averaged 0.021. 
This was due the fact that the increased insufflation air this portion 
the channel increased the cross section and correspondingly reduced the 
computed values the velocity. 

The superelevation the water surface around the 142-ft radius curve has 
been tabulated and compared with the theoretical values computed according 
Eq. results are given Table 11. The agreement between the 
actual and theoretical results sufficiently close indicate that evident 
that the flow around this curve would have been improved the supereleva- 
tion the bottom the channel accordance with the theoretical values. 


i 


— 


September, 1942 AIR ENTRAINMENT 1135 


The flow, thereby, have been kept within the channel walls and the 
overtopping the outer wall spray would have been avoided. the 
other curves along the channel the radius curvature was not sufficiently 
sharp cause overtopping the wall. Hence, the necessity for supereleva- 
tion the bottom not essential. 

the Kittitas chute the slope the energy gradient could determined 
only from velocity observations taken over relatively short detached portions 
the chute. The energy gradient computed from these data proved 
erratic, but the average the results approximated the bottom slope the 
chute, this slope was used for the energy gradient. The values the retarda- 
tion factor are larger than those obtained the tests the other three chutes, 
but the difference may due part the fact that the slope the energy 
gradient actually less than the bottom slopes, from unusual roughness 
the bottom the channel not apparent from the information available. The 
high value the 33° 10’ slope based observed areas and velocities 
undoubtedly due part the stream having been thrown into spray the 
convex vertical curve immediately above this incline (see Table 12). The same 
effect was noted the test the Hat Creek chute. 


The relation the ratio air water and the parameter has been 


plotted graphically. All the four tests described have been shown Fig. 10. 

Although there some scattering the plotted data, the results are 
reasonably consistent, considering the difficulty obtaining the velocity and 
area measurements. The graphs afford experimental demonstration the 
linear relation between the air-water ratio and the Froude number expressed 
Eq. values and for the various chutes are follows: 


Kittitas chute (upper end)......... 0.0069 0.276 
Kittitas chute (lower end)......... 0.00535 


The values are reasonably consistent with the general condition 
roughness the respective channels. Rapid Flume the wood construction 
results the lowest value the coefficient. the concrete channels Hat 
Creek has the lowest value, although the roughness the channel the 
Kittitas chute (judged from the photographs the case the latter) not 
great than Hat Creek. However, the results from these two chutes 
are close accord. The South Canal chute another project with its rough 
bottom and sinuous alinement has the highest value 

The selection the proper value left the judgment the engineer 
using the presented data guide. concrete chutes, the water dis- 
charged down the channel relatively clear, the case with spillway chutes 
from reservoirs from canals normally carrying clear water, the use value 
0.006 recommended. This coefficient slightly larger than the 
average the tests and affords some factor safety. the chute serves 
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inclined drop structure convey water irrigation canals from bench 
lands lower levels, sand carried the water will roughen the bottom, and 
the use higher value comparable that determined from the South 
Canal chute more desirable. 

More experiments over wide range channels would give added con- 
fidence the proper selection the factor and The straight-line rela- 
tion the air-water ratio dimensionless parameter indicates its application 
all ranges velocity and channel sizes from the smallest the largest. 

will noted that the application the formula large channels will 
result smaller air-water ratio for given velocities than computed small 
channels. This conforms with experience observation this phenomenon. 

One word caution must given: Eq. 24a will not apply very short 
channels which the distance travel too short entrain air; neither will 
apply streams such small size that surface tension important factor. 
general, distance travel must 100 more, unless air entrained from 
turbulent entrance conditions. Despite these minor limitations the formula 
should value designing spillway chutes. 

comparison the actual maximum depths the chutes with the height 
side-walls disclosed that large factor safety has been allowed the design 
cover the effect the little known factors under discussion. horizontal 
curves sharp radius (as, for example, Canal chute) sharp 
convex vertical curves Hat Creek and Kittitas chutes, the additional 
allowance proved insufficient. However, straight stretches 
moderate curves, the allowance for freeboard very generous. This in- 
dicated comparison maximum depth the chutes, except the afore- 
mentioned curved sections, with the height side-walls. Depths near the 
inlet before the water had accelerated have also been omitted. The average 


depths the remaining portion the chutes were well below the following 
maximum values: 


Name chute Maximum depth Height walls 
(ft) (ft) 
South Fork (10-ft width)........... 2.0 4.5 
South Fork (7-ft 3.0 3.5 5.25 


most instances more exact knowledge the performance ex- 
pected the water high velocities would have permitted reduction the 
freeboard allowance with corresponding economy the cost the structure. 


CONCLUSIONS 


The flow water steep gradients essentially different from that found 
channels ordinary slope. Velocities are greater than the critical, and 
range from minimum value per sec maximum values from 


100 per sec. The kinetic energy greatly exceeds the static pressure 
the water prism. 


‘ 
5 
4 
4 
> 
4 


K- 


4 


September, 1942 AIR ENTRAINMENT 1139 


The initial entrainment air depends great extent entrance condi- 
tions. With well-designed entrance, the water accelerates uniformly, gradu- 
ally drawing air into the water prism. This insufflation develops first along the 
side-walls and spreads gradually into the section until the admixture air and 
water covers the entire wetted area. the water drops over weir the 
inlet the chute, turbulent conditions with entrained air will exist from the 
head the chute throughout the length the channel. Similar conditions 
exist when the approach not alinement with the chute channel, when the 
inlet the chute narrowed too rapidly. Under these latter 
high, rough, cross waves develop which are reflected from side side the 
chute until the waves are dissipated the end the incline reached. 

The flow water down steep channel with entrained air does not follow 
The observed area the stream, before application the equation, must 
reduced the percentage air the admixture. The retardation factors 
the Manning Kutter formula, obtained using observed velocities and 
areas, agree with those obtained normal channel flow similar materials. 
This relation demonstrated correct, both theory and from actual 
observations. For practical purposes computations are simplified the design 
chutes the use the computed area, assuming This 
necessitates the use value smaller than that normally used, and one 
that depends for its value the velocity the water, indicated Eq. 26. 

Bernoulli’s equation then can written, starting from the critical section 
the channel the inlet and extending down the chute the outlet. For 
rectangular chutes, increments velocity are assumed and the length travel 
required gain this increment calculated Eq. resulting veloc- 
ities are then true water velocities, but such calculations yield wetted areas 
less than those actually existing. The actual cross section air plus water 
the channel can determined Eq. with the selection coefficient 


based the data presented Fig. 10, which give the ratio air water 
for values the parameter For channels other than rectangular, Eq. 


and the procedure calculation outlined afford rational solution the 
flow steep chutes. 

Convex vertical curves must not made sharper than the trajectory the 
flowing water falling under the action gravity determined from Eq. 31, 
the stream expected adhere the bottom the channel. The portion 
the water prism near the surface and the center the channel travels 
velocity from 15% 20% greater than the mean velocity the cross section 
and this maximum velocity must used computing trajectories. chutes 
short length the maximum velocity will found equal the theoretical 
velocity, resulting from the vertical drop between the inlet and the point 
vertical curvature. the side-walls are sufficient height prevent over- 
topping, serious consequence results from the jet leaving the bottom the 
channel. The evidence indicates that more energy dissipated when the 
stream cross section expanded leaving the floor the chute than would 
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otherwise the case. Under extreme conditions damage the concrete floor 
the chutes may result from cavitation. 

Flow around horizontal curves should avoided whenever possible. 
used, they should installed the head the chute before the water has had 
opportunity gain great momentum. large channels they should never 
used unless the performance the water has been verified model tests. 
Model tests are also desirable check entrance conditions. The supereleva- 
tions the water surface curves correspond with that computed the 
centrifugal theory indicated Eq. cases where the channel too 
wide with respect the depth, large percentage the water will either follow 
around the outside the curve, leaving the inner portion the floor exposed, 
some instances the slope the superelevation will not uniform across 
the wetted section, due the variation velocities across the water prism. 
sharp curves, the computed superelevation the water surface exceeds 
value tan 0.120, the superelevation the bottom recommended. 

More experimental data over wider range channel sizes and with varying 
degrees roughness would add the engineering knowledge the phenom- 
enon flow high velocities. 
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CLOSED CONDUIT FLOW 


JUN. AM. Soc. 


The general problem treated this paper the study air entrainment 
flowing water pipes, particularly reference the removal air pockets 
from water supply lines. The formation hydraulic jump pipe 
important method which air can removed, and this paper presents data 
relating the ability hydraulic jump entrain air and have carried 
away the flowing water. Quantitative data are presented relating the 
characteristics the jump with the rate air removal, and also data regarding 
the hydraulic features jumps closed conduits. 


GENERAL PROBLEM 


When air pocket formed the summit pipe line beyond 
partly open gate valve, the conduit flows only partly full, and the end 
the pocket hydraulic jump usually occurs. The air may removed, 
course, use air valve; but this not feasible the pipe line under 
sub-atmospheric pressure. any case, interest know under what 
conditions air valve may necessary or, air valve cannot not 
used, whether the air can removed the flowing water. 

Since the hydraulic jump the end air pocket common phe- 
nomenon, was decided that the first phase this investigation should 
concerned with the ability hydraulic jump entrain air. This paper 
specifically concerned with this particular item. course, the entrainment 
air the hydraulic jump not the sole item consider the general 
problem air removal from pipe lines; but important feature the 
general problem. 

The violent eddies hydraulic jump the end air pocket will 
entrain air and pump this air into the conduit beyond the jump. Whether all 
the air entrained the jump will carried along the flowing water be- 
yond the jump another problem considered. Also the problem the 
movement finite bubbles air the pipe line beyond hydraulic jump 
must considered. This may the entire problem involved the conduit 
has indicated that the slope the conduit not greater than the hydraulic 
gradient, such bubbles will moved along without any difficulty. also 
discusses the use the hydraulic jump remove air from conduits. 


Asst. Prof. and Research Engr., Iowa Inst. Hydraulic Research, Univ. Iowa, Iowa City, Iowa. 


Care, David Taylor Model Basin, Bureau Ships, Navy, Washington, C.; formerly Research 
Fellow Hydraulics, Univ. Iowa, Iowa City, Iowa. 


“The Design Pipe Kennison, Journal, A., Vol. 47, 1933, pp. 27-39. 
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EXPERIMENTAL APPARATUS AND PROCEDURE 


study the entrainment air hydraulic jump circular pipes, 
pipe line made 6-in. transparent was set the laboratory. 
The line was about long with inside diameter, 0.49 
could placed any slope 30%. Water was discharged into through 
sliding gate, the head which could controlled; thus the depth and 
velocity could varied independently. Fig. photograph the pipe 
line 30% slope. Note the hydraulic jump and the large volume air 
beyond it, indicating the jump pumping large quantity air. 


View 6-IN. TRANSPARENT INSTALLATION 


The rate air pumped the jump was measured means standard 
orifices. maintain steady conditions air was supplied from independent 
source the pocket above the jump rate equal the rate which the 
jump was pumping air. the lower end the pipe line the water flowed 
into rectangular vertical box open the top permit the escape the air. 
gate this box controlled the water level and thereby the pressure 
the pipe line. The rate water flow was measured calibrated orifices. 

Studies were made with the line horizontal and slopes 0.2%, 2.0%, 
5.0%, 10.0%, 20.0%, and 30.0%. For any particular slope various gate open- 
ings were used thus determining the depth flow ahead thejump. number 
different water discharges were run for each gate opening. Water depth and 
pressure readings were taken along the pipe line for each run. The water flow 
approaching the jump was always kept under atmospheric pressure. This 


q 
| 
if 
\ 
] 
‘ 
‘ 
= 


September, 1942 AIR ENTRAINMENT 1143 


item did not have any influence the results since higher lower pressure 
could have been maintained simply raising lowering the pressure the 
entire system. 


NOTATION 


The letter symbols used this paper are defined where they first appear 
and are assembled for convenience reference the Appendix. 


AIR DEMAND JUMP 


The results reported herein are for conditions which all the air entrained 
the jump carried out the pipe line. The general flow conditions ob- 
tained are typically illustrated Fig. 12, where the initial depth, 
the water discharge, the air discharge, and the Froude number 
the approaching flow defined where the approach velocity, 
and ye, the effective depth, the water area divided the surface width. 
The air pumped into the flowing water the jump forms large pocket 
beyond the jump which extends the point where the air leaves the pipe line. 
For low water flows and especially the low slopes single long pocket did 
not exist, but rather series relatively large bubbles was present, which 
moved downstream. Where single long pocket did occur the water flow 
this pocket was the normal depth for that particular slope. The air this 
downstream pocket was under higher pressure than the air ahead the jump, 
since the pressure head the lower stage the jump was always greater than 
the pipe diameter. 

For any given initial depth flow and slope pipe, the water discharge 
was decreased, point was reached where the air pumped into the jump ex- 
ceeded the rate which air was carried along the pipe beyond the jump. 
Visually this critical point was quite apparent. The air pumped the jump 
would form large bubble just beyond the jump. This bubble did not extend 
the end the pipe, but some the air from would carried downstream 
smaller bubbles. However, periodically this large bubble would blow out 
over the jump causing the jump move downstream; then another jump would 
form and the process would repeated. appeared that for discharges 
below this critical value the average rate air removal from the pipe line was 
controlled the flow conditions below the jump and not the jump itself. 

Considering the rate air entrainment hydraulic jump circular 
pipe, general analyses indicated that for geometrically similar conditions 
should depend the water discharge and the intensity agitation turbu- 
lence the jump. The relative intensity violence the jump dependent 
the initial Froude number, which open channel flow found 
related the energy loss the jump. The slope the pipe line, might 
also affect the air entrainment the jump; but this did not seem likely. 
Dimensional considerations thus lead the following functional relationship: 
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plotting the data obtained did not reveal any particular variation the 
air-water ratio with The only consistent variation was with 
(see Fig. 13). The air-water ratio plotted against 1), since 


y,=0.09; 


(a) 


hydraulic jump would occur when From this appears that the 
functional relationship Eq. can expressed the equation 
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should emphasized that the data Fig. are for conditions which 
the pipe line carried along and discharged all the air the jump entrained. 
For any value there was value the Froude number for which the 
pipe line would carry only part the air pumped into the water the 
jump. This critical value for the various relative depths and slopes 
studied shown Fig. 14. 


0.70 


30% 
20% 
10% 


Values 


all the runs, water depth and pressure readings were taken along the 
pipe line (see Fig. 15). most the runs the higher Froude numbers, 
the conduit beyond the jump had one long pocket air under pressure. The 
water depth this part corresponded the normal depth flow for that 
particular water discharge and slope. Since the velocity the air flow 
this pocket was relatively low the frictional losses due the flowing air were 
negligible. other words, the air the pocket was under constant pressure 
and thus the pressure line was parallel the water surface the pipe line 
which gradually became parallel the pipe line. This well shown the 
data Fig. 15. 

some the runs made the steeper slopes several consecutive stable 
jumps were obtained. Obviously they remained stable each jump had 
pump the same amount air. The pressure variation with several jumps 
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Horizontal Distance Along Pipe, Feet 


| Qu Qa 
(a) 0.147 0.720 0.060 


(b) 0.127 
0.290 
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the line indicated one the pressure plottings Fig. 15. Note how 
each jump increases the air pressure each succeeding pocket; nevertheless, 
the pressure gradient still has downward trend, must. The criterion 
which determines whether two more stable jumps can form series 
depends whether the normal depth flow, for any discharge and pipe 
slope less than y;. the ratio decreases, the closer together will 
the jumps. Since each jump must pump the same quantity air, the 
Froude number the flow approaching each jump must the same, and 
since constant, this means that the depth flow ahead each jump 
must equal y:, the depth the first jump. one jump has lower 
initial depth than the jump next below it, the first will send more air into the 
intervening air pocket than the second will take out. The two jumps will 
then move apart until the for the second jump the same that the 
first. the second jump has the smaller the action will reversed and 
the jumps will move closer together. 


VERIFICATION MOMENTUM PRINCIPLE 


The jump closed conduit such the one under consideration here can 
analyzed the use the momentum principle much the same manner 
done with the jump open conduits. Knowing the initial velocity and 
depth, the pressure the pipe when the conduit flows full just beyond the 
jump can obtained. The use the momentum principle for the hydraulic 
jump horizontal closed circular conduits was Lane, 
Am. Soc. E., and Kindsvater, Jun. Am. E., 1938. How 
the slope must considered when applying the momentum relationship 
hydraulic jumps open channels has been indicated” Bakhmeteff, 

Referring Fig. 12, the general relationship for the present problem 
follows: 


which, and are the total pressure pounds sections (1) and (2), 
respectively; and are the total momenta the water sections (1) 
and (2), respectively; the weight the jump, and the angle the 
conduit makes with the horizontal. The frictional effects the conduit 
walls are neglected. The weight the jump was computed from observations 
the length the jump, the area the pipe A,, and the unit weight the 
water corrected for entrained air. The method calculating the various 
terms Eq. was follows: 

(1) the product the water area section (1) and the average 
unit pressure this area which for sloping pipe cos which 
the distance, perpendicular the conduit, the centroid and the 


unit weight water. 
Jump Enclosed Conduits,” Lane and Kindsvater, Engineering News- 
Record, Vol. 121, 1938, pp. 815-817. 


Hydraulic Jump Sloped Boris Bakhmeteff and Arthur Matzke, Trans- 
actions, E., Vol. 60, 1938, pp. 111-118. 
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(3) computed from the measured value the maximum pressure rise 
below the jump, The product the pipe area and the average unit pres- 
sure this area gives thus: 


c 
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the entire conduit not filled with water should taken into account; this 


can best done letting the area section (2) equal which 
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the air-water flow (or volume) ratio. Therefore, 


2 
(5) equal the product the average length the jump L’, the 


specific weight the air and water mixture which can expressed 


and the pipe area: 


The value was taken being approximately equal the total ob- 
served length the jump, minus 

Fig. the two sides Eq. are plotted, and the fact that the points 
group around 45° line indicates that the data conform this expression. 
The scatter the data due difficulty measuring values like depth 
flow and pressure the high velocities with any great degree refinement. 
any case the use Eq. verified and especially the use the term 
sin since for the steep slopes this term was much 75% 


CONCLUSIONS 


The results the experiments indicate that above certain critical condi- 
tion the rate air removal from air pocket pipe line will depend the 
ability the hydraulic jump which formed entrain air. This critical 
condition, for any slope pipe and for any relative depth flow the air 
pocket, depends the value the Froude number the flow ahead the 
jump. Below this critical value the Froude number the flow beyond the 
jump will not able handle all the air entrained the jump, and thus the 
air removal will not function the jump characteristics but rather the 
hydraulic features the flow beyond the jump. 

Air flow measurements indicated that the rate air entrainment 
hydraulic jump depended largely the water discharge and the Froude 
number the flow ahead the jump. empirical relationship presented 
for the air flow. 

Data the hydraulic characteristics the hydraulic jump closed 
sloping circular conduits are presented. That the data can correlated 
the momentum relationship for the jump verified. 
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APPENDIX 


The following letter symbols, adopted for standard use this Symposium 
and subsequent discussion, conform essentially American Standard Letter 
Symbols for prepared Committee the American Standards 
Association, with Society representation, and approved the Association 
1942. 


area: 
cross section defined wetted perimeter due water 
and entrained air; 
cross section defined wetted perimeter due water 
observed cross-section area; 
cross-section area defined the inner surface closed 


conduit; 
cross section sections (1), (2), (3), are designated 
As, respectively; 


elemental area cross section; 
width: 


water-surface width trapezoidal section; 
inside diameter conduit; 

the area section (1); 
Froude’s number, subscripts referring given sections; thus 
the Froude number the flow approaching section (1); 
force; total pressure section denoted appropriate subscripts; 
West 39th Street, New York, 
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frictional resistance the channel walls for cross-section area 
elemental length; subscripts refer the corresponding 
sections; 
acceleration due gravity; 
head: 
observed head; 
velocity head, computed from the mean velocity; 
potential energy head; 
coefficient air entrainment: constant for air entrainment; 
coefficient frictional resistance; 
total observed length jump: average length the jump; 
(see Fig. length measured along channel bottom closed conduit; 
total momenta the water, subscripts indicating the section; 
g 


parameter parabola for trajectory jet 


Kutter’s Manning’s roughness coefficient: 
computed roughness coefficient used with and 
computed roughness coefficient used with and 
wetted perimeter for mixture air and water: 
computed wetted perimeter section water only; 
observed wetted perimeter section air and water; 
unit (see Eq. 36); 


air discharge; 
water discharge; 
unit discharge; flow per foot width; 
hydraulic radius A/P: 


computed hydraulic radius, 


observed hydraulic radius; 
radius the center line channel; 
sine angle between energy gradient and horizontal: 
sine angle between energy gradient and horizontal for 
computed from 
sine angle between energy gradient and horizontal for 
computed from 
sine the slope the channel bed; 
side slopes trapezoidal section; 
velocity, average section: 
average observed velocity section; 


o 


observed velocity; 

approach velocity, velocity water entering section (1) 
(subscript denoting appropriate 

velocity through element area; 
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weight the jump; 
horizontal lengths referred plane: horizontal compo- 
nent trajectory jet; 
depths, numerical subscript denoting section; vertical distances re- 
depth stream comprising mixture air and water; 
computed depth; 
effective depth; 
normal depth flow; 
observed depth; 
vertical component trajectory jet; 
initial depth depth section (1); 
mean velocity-head coefficient used obtain the average kinetic 
energy section; 
air-water flow ratio; 
unit weight water used Mr. Hall); 
angle channel bottom pipe with horizontal, degrees; 
ratio water mixture air and water the relative density 
and 
transverse slope the water surface between the inner and outer 
walls channel. 
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the practice engineering, depreciation not the only subject 
much misunderstood due the lack clear definition meaning. 
among the most important such matters, however, and one that inevitably 
will demand increasing amount attention. The appearance this 
Symposium therefore welcomed, especially since intended—in the 
words Mr. make new beginning toward formulating 

The three papers which constitute the Symposium make somewhat unequal 
contributions toward the objective set Mr. Scharff. Although the writers 
agree with many the individual statements Mr. Walls, the general effect 
his paper but confirm the complexity the subject and the divergent 
views that are possible upon many aspects it. will generally agreed 
that may well subject for engineering investigation and 
careful exercise The fact that the Symposium presented 
the Society tribute this view; but must also admitted that factors 
other than those which can classed strictly engineering matters affect 
depreciation marked degree. Therefore, although engineers have defi- 
nite contribution make what may called “the depreciation problem,”’ 
theirs not the only contribution, Mr. Walls’ paper would seem infer. 


Symposium was published November, 1941, Proceedings. Discussion this Sym- 
posium has appeared Proceedings, follows: November, 1941, Messrs. Edwin Wendt and 
Fleming, and Anson Marston; December, 1941, Messrs. William Atwood, George Goldthwaite, 
Nathan Jacobs, Flagg, Nelson Lee Smith, and Beverley Benson; February, 1942, Messrs. 
Paul Norton, Jr., Ripley, and Carroll Farwell; March, 1942, Messrs. Thomas 
David Kosh, Walker, and Lee Hyder; April, 1942, Messrs. James Ryan, John Page, 
and John Worley; and June, 1942, Messrs. John Miller, and Fred Asa Barnes. 
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The view expressed Mr. Walls that accom- 
plished fundamental. contrasts somewhat strangely with the sug- 
gestion under item (d) that the capacity enterprises make provision for 
depreciation will vary. Since depreciation accomplished fact any one 
time, the provision for not, therefore, connected with the profitability 
the enterprise whole. Provision for the latter is, rather, the contingency. 
Correspondingly, since depreciation fact, then (using this term 
generally, does Mr. Walls) will serve indication improvement 
the art currently taking place, and will guide least estimates future 
life. 

The writers are complete agreement with the views expressed Messrs. 
Crum and Winfrey. Their only regret that, due the lack interest 
the part the average citizen all governmental financial matters other than 
his they see little hope the wide adoption the near future 
the sound principles advocated this third paper. The appearance the 
paper, and the indications gives the study now being given this impor- 
tant aspect public finance, give much encouragement all who have 
heart the interests good government. One factor only does not appear 
receive full consideration this paper—namely, the fact that there are now, 
and there will increasing numbers, public works for which the benefits 
cannot expressed wholly financial terms. For such cases, some modifica- 
tion general depreciation practice will have evolved. 

Professor Grant’s paper interesting, provocative, and challenging. His 
examples are well chosen illustrate the points makes and the definitions 
suggests. the writers question certain features his paper, they 
must, they wish preface such questionings with this genuine tribute. 

Professor Grant starts his paper with three explanations terms common 
amortized cost; and depreciation,” 
the difference value between the existing old asset and hypothetical new 
asset. proceeds illustrate the use these terms the interesting cases 
that cites his paper. Example particularly good illustration the 
use the three terms; referred the following note. 

The writers wish suggest that applying the word depreciation these 
three values somewhat confusing and represents one the difficulties that 
must cleared before the matter can brought proper basis. The 
values cited for popular depreciation and appraisal depreciation differ funda- 
mentally from the value given that determined accounting depreciation, 
that both them are calculated the actual time that the purchaser 
considering buying the house. There is, therefore, little uncertainty about 
them except far the appraisal depreciation makes deduction for shorter 
life expectancy, etc. the other hand, the value given for accounting de- 
preciation was determined the basis calculations made the year 1932 
when one could tell with certainty what the conditions with respect prop- 
erty would the year 1940. 

The writers believe that compare these three values equal footing 
scarcely fair procedure. They suggest that all discussions deprecia- 
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tion clear distinction should made between actual depreciations that can 
calculated with respect past time and estimates depreciation that have 
made with respect future time. 

seems the writers that fundamental understanding depreciation 
obtained looking the trend annual unit costs which the most 
efficient and modern machine designed for use some particular field provides, 
using the lowest attractive interest rate that the class investment involved 
may direct. general, two things will observed for most products: 


(1) The price obtainable for the product will decline continuously—that 
is, the amount available from the sale the product service will decline 
from year year, and thus the sum of: Operations (O) maintenance (M) 
depreciation (D) interest (J) total must decline. This means that 
the sum the last three which must decline, the first 
item may taken constant. 


(2) The trend maintenance costs will upward the asset becomes 
older. 


Thus there are two limiting conditions: 


(a) Due the declining allowable total there will come time when the 
allowable will equal the total regardless how slowly may 
rise. this point the asset will cease have value operating unit 
because its operating costs are relatively too high, and thus has become 

(b) Due the gradually rising there will come time when the 
total alone will equal regardless how slowly may fall. This may 
defined statement the extent dilapidation. 


Now any formula arrangement for estimating either value deprecia- 
tion allowances must satisfy the series figures associated with the two trends 
(1) and (2) and the two limits associated with (a) and yields basic 
rule that the total the charges must general steady, 
decline from year year, and must never rise. 

Now apply this law Professor Grant’s formula under “Engineering 
Economy Studies Relative Proposed New Investment.” Although, 
Professor Grant indicates, may possible arrange series annual 
payments depreciation account, which would end establishing the same 
capital value given this formula 1), and not involve uniform annual 
paymertts, such series would extremely limited form and not all 
likely fit very many cases found practice. Eq. thus essentially bound 
with the idea equal annual payments. The aforementioned rule thus 
limits its application those cases where, for practical purposes, the curve 
straight, horizontal line. 

Eq. used with asset having other than horizontal, straight-line 
for then the totals must averaged, which case the 
result will only approximate. the curve not horizontal and also 
not straight line, the averaged result may considerably error. The 
latter condition frequently found. such calculation thus will 
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essential determine the curve and take average value that too. 
This average will open the same risks that itself. Professor 
Grant has apparently taken into his value for annual disbursements 
Examples and This equal and presumably averaged figure. 
general, this formula would seem applicable those assets having 
horizontal, straight-line relationship for and the same time having 
equal zero, forming horizontal straight line. thus not generally 
applicable accurate formula. 


TABLE IRREGULAR MAINTENANCE TREND AND DECREASING 


Cost 


Depreciation Maintenance Interest Total 
Beginning End 
year year 


2,000 
1,600 
1,300 
1,100 
950 
800 
650 
525 
400 
275 
175 


Total 

Annual 
Maintenance deducted ($171) 
Sum added for interest the beginning the year 
Net cost the beginning the year 


(operating costs, being constant, have been omitted), and which irregular 
trend maintenance costs will noted that this series 
capital recovery factor 8.576. The average annual net total $216 inter- 
est taken the investment the beginning the year, $211 interest 
taken the average investment during the year, will support valuation 
$1,810 the latter. These values are 7.5% and 9.5%, respectively, below the 
actual worth the asset. The formula assumes the former case. thus 
would seem possible obtain errors much 10% more with this 
formula. will noted that this error inherent the formula, and 
consequently additional any errors the basic data, upon which the 
lations were made. The writers realize that Eq. provides most valuable 
means determining the value many items property. The object 
this criticism indicate that errors are possible, nevertheless. 


and Its Influence Unit Hart, undergraduate thesis (unpublished) 
presented the Faculty Applied Science and Engineering, Univ. Toronto, Toronto, Ont., Canada, 
1939. 


1,600 1,800 400 108 508 
1,300 1,450 300 462 
1,100 1,250 200 150 425 
950 1,025 150 175 387 
800 875 150 175 378 
650 725 150 175 369 
525 588 125 200 360 
400 125 200 353 
275 338 125 200 340 
175 225 100 225 339 
125 100 225 332 
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Discussion 


capacities single piles and pile groups has been controversial 
subject among foundation engineers for many years. The experiments de- 
scribed this paper might have served add materially the available 
knowledge this problem. Unfortunately, appears though too much 
effort was devoted the development unsound mathematical theory and 
too little effort the development carefully thought out test program that 
would include field measurements all the quantities that are important 
problems this kind. 

Referring the mathematical theory, the effort use the Bous- 
sinesq point-load equations simply adds another chapter the long story 
the incorrect application these equations all kinds foundation problems. 
Several years ago, the writer the original derivation the Bous- 
sinesq equations together with statement the fundamental assumptions 
which the derivation based. Whenever these equations are applied 
actual problem, necessary keep these basic assumptions mind and 
make sure that the conditions the problem not conflict with the assump- 
tions under which the equations were derived. 

For example, one Boussinesq’s fundamental assumptions that the load 
applied the plane boundary the semi-infinite, elastic, isotropic, solid and 
not some point below the surface the solid. The author has used the 
Boussinesq equations for loads that are applied below the surface the ground 
although has already been explained Raymond Assoc. 
Am. Soe. E., why the Boussinesq equations cannot used this manner 


paper Frank Masters, Am. Soc. E., was published November, 1941, 
Proceedings. Discussion this paper has appeared Proceedings, follows: February, 1942, 
Connor, Am. E.; April, 1942, Messrs. Glenn Woodruff, Jacob Feld, Greulich, 
and Paxson; May, 1942, Raymond Mindlin, Assoc. Am. Soc. E.; and June, 1942, 


John Solomon, Esq. 
Dist. Mgr., Raymond Concrete Pile Co., Chicago, 
Received the Secretary August 1942. 
Transactions, Am. Soc. E., Vol. 103 (1938), 324. 
the Distribution Stress Around Mindlin, Proceedings, International Con- 
Soil Mechanics and Foundation Eng., Cambridge, Mass., Vol. IIT, 1936, 71. 
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determine the distribution stresses around pile. Furthermore, the 
author does not appear have given any consideration some the other 
fundamental assumptions the mathematical theory elasticity which 
the Boussinesq equations are based. particular, the paper says nothing 
about strains displacements, although these are just important the 
stresses problems this kind. Also, the extension his analysis 
from the single pile the pile group presupposes the validity the principle 
superposition. considerable interest examine the proposed 
theoretical solution from the standpoint displacements and the principle 
superposition. 

well known, the entire mathematical theory elasticity based 
law concerning the proportionality stress and strain. Whenever 
force applied elastic body, some kind displacement occurs result 
that force. The problem under consideration exception this rule 
and the stresses determined the author would necessarily accompanied 
displacements the soil surrounding the pile. Since the author has as- 
sumed the soil incompressible (Poisson’s ratio 0.5), the horizontal 
tangential stress vanishes everywhere and the two normal stresses given 
Eqs. and will determine the horizontal radial displacements means 
the equation 


which the radial displacement, and Young’s modulus the soil. 


The vertical displacements, may obtained similar manner from the 
author’s stress equations. 


Fig. represents any meridional plane the soil through the vertical 
axis the pile The solid lines represent imaginary rectangular network 
existing the ground after the pile driven and before any load applied 
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the pile head. When the load, placed pile stresses are set the 
soil around the pile, and these stresses are accompanied necessarily certain 
vertical and radial displacements the soil. The dotted lines Fig. are 
qualitative representation the distortions that would occur the rec- 
tangular network due the displacements determined the author’s theo- 
retical stresses. The actual amounts the displacements would depend 
the magnitude the load, and the numerical value Young’s modulus, 
the soil. All points the ground surface would displaced vertically 
downward and there would radial displacements the ground surface 
because the assumption incompressibility (Poisson’s ratio 
example, soil particle originally would move Within the soil, all 
vertical cylindrical surfaces concentric with the vertical axis the pile would 
bulged radially outward and all horizontal planes would bent downward. 
soil particle originally point would move point 

This system displacements necessary consequence the author’s 
stress system. the stresses the soil were those determined the author, 
the displacements would have the type shown Fig. 18. con- 
ceivable that such system stresses and displacements could exist the soil 
the vicinity single pile. However, the soil around the pile were sub- 
jected restraints that would affect the displacement pattern, these restraints 
would also affect the stress pattern. This follows inevitably from the fact that 
the entire mathematical theory elasticity based definite set stress- 
strain relationships. Boussinesq’s stress equations were derived with the aid 
the basic differential equations which relate the stresses the space deriva- 
tives the displacements. Accordingly, the Boussinesq stress pattern can 
exist only when accompanied certain definite pattern displacements. 

For the case illustrated Fig. 18, the question arises what will 
happen the stresses and displacements around pile when that pile sur- 
rounded group other piles such and The elastic properties 
the piles will certainly quite different from those the surrounding soil 
regardless whether the piles are made wood, concrete, steel. effect, 
the piles will act vertical reinforcement the soil. They will distort the 
displacement pattern which determined the author’s stress equations and 
this turn will change the stresses themselves. Accordingly, not 
expected that the stress distribution around single pile can extended 
pile group simple superposition done the author his theoretical 
analysis. Furthermore, should noted that the combination soil and 
piles within the pile group means the isotropic material postulated 
Boussinesq his analysis. The presence the piles produces compound 
material consisting hard spots (piles) embedded matrix much softer 
material (soil). 

The principle superposition sometimes used determine the stress 
distribution foundation problems. For example, large mat foundation 
often divided into small segments and the stress distribution for each segment 
computed the segment were point The stresses produced 
Proceedings, Am. Soc. E., May, 1933, 784. 
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each segment are then added simple superposition obtain the stress dis- 
tribution for the entire mat. This application the principle superposition 
mat foundation something quite different from the application that 
principle pile group. the case the pile group, each additional pile 
affects the elastic properties the underground. Even though the soil around 
single pile might considered having the properties elastic isotropic 
material, the introduction additional piles destroys the isotropic properties 
and materially changes the elastic properties. Accordingly, the 
method applying the principle superposition pile group incorrect 
procedure. 

Another factor this problem which the author has paid attention 
the contraction the pile the pile place the ground and 
before any load has been applied the pile head, there shearing stress 
the soil such that given Eq. 2c. When load applied the pile head, 
there downward movement the upper end the pile due the elastic 
shortening the pile. the pile starts move with respect the surround- 
ing soil, the movement mobilizes the shearing resistance the soil. Since the 
movement begins the top the pile, the mobilization shearing resistance 
the soil also begins the top the pile near the ground surface. 
small load applied long friction pile, all the reaction for the load 
developed the soil around the upper part the pile. The load produces 
stress the lower part the pile and consequently elastic contraction and 
mobilization shearing resistance the soil surrounding the lower part 
the pile. the load increased the stress the pile approaches the pile 
point but does not actually reach the point until the load attains cer- 
tain value which determined the elastic properties the pile and the 
surrounding soil. 

Fig. the author has computed the shearing stresses around single 
pile long carrying 10-ton load. value 1.5 used the 
calculation. Actually the value would very much greater than 1.5 for 
the conditions shown the figure because little load would delivered 
the soil the bottom the pile. load delivered the soil the 
bottom the pile, the value becomes infinite. matter fact, the 
determination the numerical value should the starting point for any 
theory friction pile behavior. The real problem determine how the 
delivers its load the surrounding soil. author states that function 
pile length, pile taper, and internal friction the soil. has overlooked 
the fact that also function the elastic properties the pile well 
function the load that applied the pile head. the pile were per- 
fectly rigid that elastic shortening the pile would occur under any load, 
might possible have 1.5 for all loads and have 
distribution similar that shown around the single pile Fig. 
actual piles are elastic and will therefore shorten under compressive stresses, 
this elastic shortening will influence the transfer load from the pile the soil. 
The value therefore, will depend the elastic properties the pile 
the amount the applied load well the three factors mentioned 
the author. 
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The usual procedure the case pile load test place the load top 
the pile and then measure the downward movement the head the 
pile for each load increment. not possible determine the elastic shorten- 
ing the pile from such measurements because nothing known about the 
movement the pile point. Neither can such measurements furnish any in- 
formation the manner which the transfers the load the surround- 
ing soil. the other hand, the actual elastic shortening the pile were 
determined suitable measurements, would possible determine from 
these measurements approximately how the pile was transferring its load the 
soil. For example, simple matter compute what the elastic shortening 
the pile would were free-standing strut with part the load sup- 
ported friction the sides the strut. When any part the load 
supported friction along the sides the pile, the elastic shortening always 
less than would for the same pile carrying the same load free-standing 
strut. When all the load carried side friction the author’s prob- 
lem, possible compute the elastic shortening the pile for any assumed 
distribution friction along the sides the pile. This can done even the 
friction distribution assumed discontinuous due soil stratification. 
Conversely, the elastic shortening the pile carefully measured the field, 
possible deduce the friction distribution from the field measurements. 

Experiments this kind have been made The Netherlands 
and others. The piles tested The Netherlands were pre- 
cast concrete piles. When they were cast, piece 1.5-in. pipe extending the 
full length the pile was placed the vertical axis. 1-in. rod was inserted 
this pipe and fastened the pile point. The test load was applied the 
driven pile such manner that was possible measure the deflection 
the pile head and also the movements the 1-in. rod which gave the deflections 
the pile point. The difference between the two readings for any given load 
represented the elastic shortening the pile that load. Fig. 19, curves 
and are reproductions pair load-settlement curves measured 
Boonstra. The solid line shows the settlement the pile head and the broken 
line that the pile point. The specimen was reinforced pre-cast concrete 
pile octagonal cross section, in. diameter and was driven 
through soft clay and peat, fine sand, and then into sand and 
gravel for depth about ft. The condition was similar that pile T20 
the author’s tests. 

can seen from Fig. 19, there was practically movement the pile 
point until the load amounted tons although the head the pile 
began move soon the loading was started. The vertical distances 
between the two curves, and represent the elastic shortening the pile, 
but order show this relationship more clearly, curve Fig. 19, was plotted 
with loads abscissas and elastic shortening ordinates. easily seen 
from this figure that the relationship between load and elastic shortening not 
the distribution load transfer from pile soil were constant 


and van Heipalen Franx, Polytechnisch Weekblad, 1934, Nos. 23, 
oO. 


Beschouwingen over den Puntweerstand van Palen,” Boonstra, Ingenieur, 
ebruary 23, 1940, Vol. 55, No. B-33. 
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throughout the entire load range constant), and the pile material were 
perfectly elastic, curve Fig. 19, would straight line. the elastic 
properties the pile, could expected behave under load like reinforced 
concrete column. The stress-strain curve for such column might slightly 
curved but the curve would concave toward the strain axis, whereas most 
curve Fig. 19, concave 

toward the stress axis, 

Accordingly, the only con- 

Settlement clusion that can drawn 


Pile Point 
distribution load trans- 


fer from pile soil varies 
throughout the load range. 
constant. the load 
increased, the transfer 
load from pile soil dis- 
tributes itself various 
ways depending the 
elastic properties the 
pile and the soil and the 
amount the load. The 
sudden break curve 
about 150 tons indicates 
that the upper part the 
pile probably broke loose 
from the surrounding soil 
entirely that practically 
all the applied load was 
resisted the lower end 
the pile which was 
the sand and gravel. 
Load, Tons Consideration this 
type test data raises 
serious question the 
validity the figures presented Table All the pile loads given this 
table are computed with 1.5 the assumption that the distribution 
load transfer from pile soil constant over wide range loads and pile 
lengths. The author’s computation implies that independent variable 
which numerical values can assigned for the purpose computing pile 
loads. Actually, the situation just the reverse. The amount the load 
influences the manner which the pile transfers load the soil that the 
numerical value depends the load. 

This aspect the problem something that could easily have been in- 
vestigated the Morganza tests. would have been possible set three 
four rods hollow concrete pile, such pile C9, that settlement observa- 
tions could have been recorded not only the pile head but also the pile 
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point and several places along the length the pile. the strains the 
pile itself had been measured, would possible determine how the pile 
was transferring load the soil each load increment. 

Because has completely neglected the entire question strains and dis- 
placements, the author has developed mathematical analysis that theo- 
retically unsound and practical value. The analysis based the 
mathematical theory elasticity but nowhere the paper there any mention 
the elastic properties the piles themselves. The author has treated the 
piles though they did not exist except imaginary vertical line loads the 
soil. Unfortunately, this method has led theoretical stress distribution 
which, the writer’s opinion, must also classified imaginary. The piles 
are actual physical bodies having certain elastic properties, and neither the 
single pile nor the pile group can analyzed without taking these elastic 
properties into account. 

Throughout his theoretical analysis pile groups the author has assumed 
that each pile the group would subjected the same load all the 
others. has made this same assumption connection with the groups that 
were actually tested. This question load distribution over the pile heads 
important factor pile group behavior. The author has referred this 
subject his conclusions remarking that the piles under the outer edges 
stiff pier base would subjected greater loads than those under the center. 
This remark perfectly true but the question arises what was happening 
this respect during the Morganza group tests. 

Fig. 2(c) shows the load arrangement that was used for the 16-pile group. 
each pile the group were subjected the same load, the center piles could 
expected settle the most, the edge piles would settle somewhat less, and 
the corner piles would settle least all. The photograph indicates that the 
applied load was quite rigid that could scarcely expected warp 
manner that would conform with the settlements that would produced 
equal pile loads. Instead, the rigidity the test apparatus would require 
equal settlement the pile heads and this turn would affect the distribution 
load the piles. perfectly rigid structure would put the greatest load 
the corner piles and the least load the center piles. 

This another aspect the problem that could have been investigated 
the Morganza group tests. Fig. indicates that the pile heads were several 
feet above the ground surface and would have been possible make strain- 


measurements the upper ends the piles determine how much load 


each pile was carrying. There justification for the assumption that all 
the piles group were carrying equal loads since effort was made pro- 
duce this condition. insure equal loads all piles group would 
support the test load independently the pile group and then 
load all the piles simultaneously with individual hydraulic jacks reacting 
against the under side the load platform. should kept mind that the 
distribution pressure the contact plane between any two solid bodies will 
depend the relative elastic properties the two bodies. This applies 
loaded pile groups well loaded masses earth other materials. 


| 

lly 

the 

pile 

ble 

pile 

oad 
the 

pile 


1166 CUMMINGS FRICTION PILE FOUNDATIONS Discussions 


Furthermore, should noted that this phenomenon could explain the 
behavior the Morganza pile groups without reference the theory 
stress distribution. The practically rigid structures built top the pile 
groups could have had the effect overloading the corner and edge piles 
that the groups would have failed leaning toward corner side described 
the author. theory pile group behavior can considered being 
complete even approximately correct unless this phase the problem 
taken into account. 

demonstrate the accuracy his theory, the author has prepared Table 
which indicates remarkably close agreement between theoretical and actual 
failure loads. However, the writer finds difficult follow the line reason- 
ing used the author his selection shear values which the comparisons 
are based. average value used for four the group tests, and the other 
three are compared with values determined for certain individual 
ferring Fig. these two profiles represent cross sections the floodway 
miles wide. The profiles themselves are some miles apart along the length 
the floodway. The single piles, well the pile groups, are therefore 
over distances that are measured miles, and the first question 
decided whether the comparisons should based averages whether 
each group test should compared with the results the nearest single pile. 

considerable interest examine Table with these ideas mind. 
check the group, the author used pile which was 130 from the 
group. However, pile T10 was only 294 from the group. Accordingly, 
the proper way handle this problem means averages, would seem 
reasonable average and T10 order check this group. Such pro- 
cedure would not check the theory closely because the average friction 
developed and T10 was 792 per ft, whereas the author used 685 
per from alone. the other hand, the proper way handle this 
problem comparing each group with the nearest single pile, groups 
and T37—T52 should compared with T20 alone and not 
with average that includes piles that are miles away from these groups. The 
distance from group group 1,958 and the third group 
well pile T20 are between the two. However, comparison these 
groups with pile T20 would not check the theory well because pile 
developed friction value only 680 per ft, whereas the author 
these groups using the general average value 778 per ft. 

The data given Table indicate that the tests provided almost perfect 
check the author’s theory. The agreement between theoretical and actual 
loads really too perfect view the fact that the theory unsound from 
the point view both mathematics and theoretical the 
writer’s opinion, the theory analogous the old Ptolemaic theory 
geocentric solar system. For 1,400 years the Ptolemaic system was 
cepted generally because appeared able account for the motions 
the planets. After Galileo invented the astronomical telescope, 
were able stop speculating and start measuring. The very first measure 
ments caused the abandonment Ptolemy’s geocentric system favor the 
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modern heliocentric system. The fact that Galileo got into trouble for measur- 
ing things and upsetting ancient theory should not discourage modern 
engineers from making all the measurements that can possibly made 
given experiment. theory based insufficient data can appear prac- 
tically perfect and yet completely wrong. 

The writer joins the author the hope that engineers interested pile 
foundations will continue make tests determine how such foundations 
behave under load. However, the writer considers necessary supplement 
that hope with several suggestions. develop adequate theory for the 
behavior single friction pile, will necessary determine, means 
suitable field measurements, how the pile transfers its load the surrounding 
soil. order develop adequate theory for the behavior group 
friction piles, field measurements must made determine the distribution 
load the pile heads addition the measurements that are required 
determine the distribution load transfer from pile soil. Field measure- 
ments that include only record the settlement and rebound the pile head 
are inadequate for the development theory friction pile behavior. 


Correction for Transactions: June, 1942, Proceedings, Eq. 7b, change 
page 1049, line 10, change Table 16, change 
“0.035” and and respectively. 
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DISCUSSIONS 


DESIGN ST. GEORGES TIED ARCH SPAN 


Discussion 


arch rib tied with stiffening girder extremely complicated problem. 
the technical literature very little can found this subject. com- 
puting the stresses the St. Georges arch bridge, Professor Garrelts has adopted 
various assumptions simplify the computations. The writer proposes 
clarify, further, the limitations assumptions adopted the light economic 
aspects tied arch bridges with stiffening girders. 

From Fig. evident that the stiffening girder the St. Georges bridge 
connected the arch rib rigid joint, whereas, the development 
Eq. for the computation stresses the arch rib, the author assumed 
that the stiffening girder connected the arch rib frictionless hinge. 
This assumption has negligible effect computation the horizontal 

ever, the equations developed with this assumption are not adequate for the 
design the arch rib sections the ends and their connections the stiffening 
girder. 

The axial stresses the stiffening girder, produced the arch rib thrust, 
are generally small compared with the bending moment stresses. Therefore, 
sections stiffening girder are usually varied according variations the 

bending moment stresses along the girder axis. developing Eq. for the 
horizontal thrust the arch rib, Professor Garrelts, however, assumed that 
the stiffening girder section constant for the entire length. Evidently with 
this approximation the refinement considering that the arch rib section 
varies directly cos (Eqs. 16) has doubtful meaning. 

the design tied arch rib bridges, important consideration the 
design the lateral stiffening system. Often economy material gained 

paper Garrelts, Assoc. Am. Soc. E., was published December, 194! 
Proceedings. this paper has appeared Proceedings, follows: February, 1942, 


Abbett, Am. Soc. E.; April, 1942, Jacob Karol, Esq.; May, 1942, Alexander Dodge, 
and June, 1942, Carlos Bejarano, Jun. Am. Soc. 


Associate Bridge Engr., Bridge Dept., Div. Highways, State Dept. Public Works, 


10a Received the Secretary June 22, 1942. 
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careful computations the stresses the arch rib and stiffening girder may 
easily lost the wasteful approximate design the lateral stiffening system. 
hoped that, his closing comments, the author will discuss method 
designing the lateral stiffening system and will present some principles 
govern the selection panel lengths. 

The author commended for making valuable contribution this 
subject. 


formulas for the design the tied arch derived Professor Garrelts, the 
question immediately arises whether the validity these formulas not 
subject the same limitations the formulas the elementary theory 
the suspension bridge, particularly taking into consideration the relatively large 
span the structure actually designed. High deformability common 
feature all structures which shallow girder suspended from, supported 
by, flexible member provides the bending resistance for unevenly distributed 
live load. this deformability that renders the methods the elastic 
theory structures ineffective the design suspension bridges well 
arches, particularly those with flexible ribs long 

general, designers are quite familiar with the fact that considerable 
differences are expected between the bending moments the stiffening 
truss suspension bridge, computed the elastic theory, and the actual 
moments. They are less familiar with the effects the deflections upon the 
bending moments arches. Therefore, few quite general remarks regarding 
the limitations imposed upon the validity elementary formulas the flexi- 
bility structure may not devoid interest. respect these 
comments deprecate the method design presented the author which, 
incidentally, practically unaffected considerations deflections. 

non-tied flexible systems (suspension bridges and arches) the bending 
moment due live load can expressed 


which and denote statical values due live load and dead load, respec- 
tively, and denotes the deflection the cross section the girder considered. 
Hence, the moments the stiffening girder computed the elastic theory 
are either reduced (suspension bridges) increased (arches) the amount 
Owing the considerable value the dead-load thrust for 
large-span structures, the deviation from the elastic theory may attain com- 
paratively high proportion the total bending moment. 

can easily shown that self-tied structures the bending moment 
the girder (Fig. 14), 


Engr., Marine Trust, Tel Aviv Port, Tel Aviv, Palestine. 
Received the Secretary August 1942. 


Bridge and Structural Engineering, Vol. 1935, 100. 
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not affected deflection the flexible member (the arch rib the cable) 
has bending resistance. However, the bending resistance cannot 
neglected (as most arches), the effect the deflections upon the value the 
bending moments makes itself felt the rib moments. The real values 


LIT 


these moments, therefore, may exceed, considerably, those obtained simple 
distribution the total bending moments between rib and girder according 
the respective moments inertia, proposed Professor Garrelts. 
particularly the buckling resistance the rib that may impaired 
bending moment exceeding the design value. 
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DISCUSSIONS 


HYDRAULIC DESIGN DROP STRUCTURES 
FOR GULLY CONTROL 


Discussion 


Messrs. CHRISTENSEN AND DWIGHT GUNDER, AND 
Boris BAKHMETEFF AND NICHOLAS FEODOROFF 


exceedingly large number drop structures needed gully- 
and irrigation work reason enough for laboratory research this 
kind. The actual value the research becomes very evident when one con- 
siders the large number failures drops already built and also the large 
number which have been overbuilt order insure safety. The value this 
research has already been demonstrated; but, another wave conservation 
and reclamation sweeps the United States after World War II, these researches 
will pay unusually large dividends. 

The “half technique used these experiments deserves further 
elaboration. The two principal advantages this technique are: 


(1) permits accurate observation conditions the center line the 
channel; and 


(2) With given model size, the quantity water required reduced 
half. 


The objections which might raised are: 


(1) Extraneous friction forces are introduced the glass plate; and 
(2) The turbulent action normally found the center the stream 
eliminated the rigid glass plate serving boundary. 


The two advantages enumerated are obvious and need further discus- 
sion. The following remarks will confined the disadvantages mentioned. 


was published January, 1942, Proceedings. Discussion this has appeared Proceedings, 
follows: April, 1942, John Hedberg, Assoc. Am. Soc. E.; May, 1942, Messrs. Standish 
Hall, and Christiansen; and June, 1942, Walter Wilson, Assoc. Am. Soc. 


Dean Eng., and Chairman Eng. Div., Experiment Station, Colorado State College Agriculture 
and Mechanic Arts, Fort Collins, Colo. 


Assoc. Prof. Civ. Eng., Colorado State College Agriculture and Mechanic Arts, Fort Collins, 
O10, 


Received the Secretary June 29, 1942. 
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The friction forces introduced the glass plate are source error, but 
every model study this kind friction forces are always source error 
because Froude’s number used criterion for similarity. With this num- 
ber criterion, one assured the proportionality the gravity and inertia 
forces the prototype and model. The friction forces, the other hand, are 
entirely neglected and considered insignificant far results are concerned. 
The smaller the model, the more serious the error introduced this assump- 
tion. quite likely that the increase scale made possible the 
technique more than compensates for the extra friction force introduced 
the presence the glass plate the center line. 

order appraise the importance the second objection the 
model” technique, attention called natural channels which drop struc- 
tures have been erected. From study these natural channels, obvious 
that the excessive scour which sometimes occurs just below the structure 
caused lateral currents the stream itself and not the haphazard 
turbulence which present most streams. normal stream turbulence 
were the cause the excessive scour, would not localized just downstream 
from the drop structure but would occur equally along the entire reach the 
stream. These lateral currents and excessive scour are invariably associated 
with unlevel water surface and the pressure gradients within the body 
water which cause continuous currents having lateral velocity components. 
most drop structures, vertical plane symmetry exists, which includes 
the center line the structure. Not only the structure itself symmetrical 
about this plane, but also the mean flow and the pressure field within the flow. 
such structures vertical plane may inserted the center line, and the 
mean flow will essentially unchanged. true that the turbulent flow 
components will reduced, but these components play only secondary part 
the excessive erosion that sometimes takes place just downstream from the 
drop structure. 

The requirements necessary for the “half model’’ technique would seem 
be: 


(a) Geometric symmetry about the center line the structure; 
(b) Kinematic and dynamic symmetry the mean flow; and 
(c) Stability the flow pattern normal the plate the center line. 


Whenever these three conditions not obtain, the model’’ technique 
should not used. 

the paper, inference was made that gully-control structures designed 
according the principles set forth would require very little maintenance. 
For series drop structures permanently and completely successful, 
each structure must not only soundly designed, but the reaches between 
the structures must also completely stable. The question what 
stable grade for the gully must first answered before the height and spacing 
the various drops can determined. The importance accurate deter- 
mination stable grade under the varying flow conditions the gully 
should emphasized part the over-all problem gully control means 
structures. This problem further complicated the fact that the flow 
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conditions gullies may change rather rapidly due various erosion-control 
practices. 

The writers would like request elaboration certain parts the paper. 
The statement made beginning nine lines from Eq. height the 
end sill the stilling pool necessary develop energy-dissipating action 
apparently determined the depth and energy content the flow the 
somewhat inconsistent with Eq. 7c. This formula would 
for the same sill height regardless the height drop. The height 
drop, the other hand, has much with the energy content the flow 
the pool. seems indicate that the height the sill determined 
Eq. concerned, the question energy does not enter. 

The last sentence the third design rule (see heading Design 
Rules and Formulas’’) seems somewhat inconsistent with the material 
contained the article proper. The discussion involving Eqs. well 


Fig. 14, seems indicate that Eq. would apply for values smaller 


than the use this formula hazardous for values smaller than 


then certain part Fig. should deleted. 

should noted also that, whereas the derivation Eq. starts with 
what appears rational attack, the final equation almost completely 
empirical, especially that contains purely empirical term. Further, 
the expression for the thickness the falling sheet (given proportional 


relatively small. can seen once that this expression, for the 


thickness approaches infinity becomes smaller. One thus inclined 
suspect that any equation using this thickness will inaccurate the region 


where large. This just what was found the authors. 


slightly different attack this problem will now given. the depth 
the jet the floor, and the other assumptions are given the paper, 
found that the horizontal distance from the crest the middle the im- 
pinging sheet The thickness the falling sheet under 
the same conditions given the continuity equation 


and the horizontal width this sheet given 
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which the angle between the vertical and the direction the sheet 


filaments. 


one now assumes that the length the apron linearly dependent upon 


each these equations and are empirical 


Since the writers did not have the original data available for computations, 


rough comparison the apron lengths given Eq. and Eq. was 


made assuming This gives 


paper. Using 3.5, the values C’; were computed and are compared 


TABLE COMPARED WITH 


Coefficient 
0 0.2 0.4 0.6 0.8 1.0 1.26 
2.5 2.7 3.0 3.3 3.75 4.3 5.0¢ 


These values are the range where the authors claim 
that should not used. The writers not know what 
values should used. 


with the values estimated 
Fig. (see Table 4). These 
results show maximum varia- 
the values (excluding the 


doubtful value 1.2). Since 


this accomplished without the 
inclusion purely empirical 
3 
term such the 0.7 
6a, might considered 
somewhat better than 6a. 


should noted that the foregoing results were obtained with just one em- 
pirical constant. Perhaps even better results could obtained adjusting 


both and 


The writers not wish claim that this rigorous rational derivation 
the expression for apron length but felt that may helpful 
further study the problem the design gully-control structures, which the 


authors have admirably attacked. 
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paper presents interesting endeavor proportion structural 
design basing the dimensioning data derived from detailed observations 
the physical aspects the flow. evident from Fig. that the dimensions 
the stilling pool are principally conditioned the horizontal distance from 
the wall which the falling live vein strikes the apron, and implicitly the 
angle which such striking effected. The latter seems particular 
determine the direction followed the rebounding vein, and thus the length 


(a) r=1 Inch 
Outflow Without Sill 


0.573 
0.343 
0.185 
0.057 


Sill 


r=3.5 Inches 
Outflow Without Sill 


(for 1.02, Sill 


102 0.319 1.09 2.73 
0.343 0.154 0.530 5.65 
0.185 0.102 8.52 
0.159 188 


— 
q 


ry 


> 
Feet Below Crest 


Horizontal Distance Feet, 
PROFILES WITH AND WITHOUT SILL 


Prof., Civ. Eng., Columbia Univ., New York, 
Research Associate, Civ. Eng. Dept., Columbia Univ., New York, 
Received the Secretary July 1942. 
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the downstream channel required for the flow assume more less normal 
features. 

stands reason that economy construction, well more favorable 
outflow, would secured the falling vein were caused strike the apron close 


the wall and thus larger angle. This could achieved resorting 
the so-called type nappes, which would follow the surface the 
wall and which would strike the floor the very toe the structure. The 
general disadvantage such nappes, course, their instability. fact, 
when the brink the fall has sharp wedge, Figs. and the nappe, 
except for extremely low heads, will tend detach itself from the wall, changing 
the customary outline. 

The situation can remedied, however, properly rounding the exit 
wedge. Indeed, observations show that, the brink the fall shaped 
inserting arc circle, the adherent nappes become quite stable. 
relatively small rounding seems serve the purpose securing dependable 
vertical falling nappe, and that comparatively high rates flow. 
tion illustrated Fig. 20, which depicts flow profiles obtained the Fluid 
Mechanics Laboratory Columbia University New York, The 
apparatus used the experimenting has been well described elsewhere.” 
Fig. 20(a) the wedge rounded radius in., whereas, Fig. 
both figures the solid lines show the downstream surface profiles 
free undisturbed outflow. The broken lines demonstrate the eventual effect 


Proceedings, Am. Soc. E., April, 1942, 637. 
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interposed 2-in. sill forming stilling basin. The unit width discharges 
are reflected the parametric values the critical depth and the last two 


columns present the dimensionless ratios and being the height the 


fall. Obviously, the foregoing ratios govern the flow forms. exhaustive 
studies were made establish the exact relationship which the adherent 
nappe changes its form. both instances, however, the profiles plotted 
were close the safe limit, further increases the flow brought about signs 
instability. 

characteristic feature the comparatively small rebounding effect ac- 
companying vertically dropping veins, especially the presence asill. Also, 
the latter case, the outflow into the downstream channel was relatively calm, 
where the size the stilling pool was obviously insufficient. fact, the general 
impression that, with adherent vertical nappes, the quieting effect exercised 
the stilling pool more decisive than with detached nappes. The general 
physical aspects are partly illustrated Fig. 21. The underpressures de- 
veloped next the rounded surface the wedge were relatively small. 
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DISCUSSIONS 


THE GREASE PROBLEM SEWAGE TREATMENT 


Discussion 


DENIS DICKINSON, AND CHARLES VELZY 


used for the determination grease sewage associated closely with the 
definition the determination made order estimate 
the difficulties that may arise the purification the sewage and the subse- 
quent disposal the sludge, then there would appear sufficient reason 
for the adoption chloroform the solvent, since affirmed that this 
liquid also extracts resinous substances that are likely more difficult 
break down biologically than the less complex fatty acids. Especially does 
this apply those cases which disposal left auto-purification the 
sea water course. the same grounds, chloroform should adopted 
for testing the efficiency any grease removal plant. the other hand, 
desired assess the practicability grease recovery, then more 
selective solvent preferable. 

interesting note that the ratio the grease determined 
chloroform that determined petroleum ether frequently constant for 
sludges from any one sewage, although the data available are not sufficient 
justify any generalization this point. would appear that there much 
learned from study the characteristics the greases extracted 
the various solvents. 

The recovery grease has been the practice several works England 
for many years, and there are other cases where floating fat skimmed off 
(usually scoops), drained, and sold for the manufacture wagon grease, 
etc. The process for the recovery grease solvent extraction the sludge 
has further advantage that the sludge finally grease-free and weed-free, 
all seeds having been destroyed during the extraction The latter 


paper Almon Fales and Samuel Greeley, Members, Am. Soc. E., was published 
February, 1942, Proceedings. Discussion this paper has appeared Proceedings, follows: 
1942, Arthur Weston, Am. Soc. E.; May, 1942, Harvey Ludwig, Jun. Am. Soc. 
and June, 1942, Mahlie, Esq. 
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Halifax Sewage Disposal Works,” anonymous, Journal and Proceedings, Inst. 
Purification, 1937, 60. 


1178 


‘ 
' 
‘ 
c 
7 


we 


September, 1942 VELZY GREASE SEWAGE 1179 


very important point increases considerably the market value the 
sludge fertilizer. 

The practice burying grease, which mentioned the paper, cannot 
condemned too strongly, feared that many cases merely 
defers the day which final disposal must undertaken. 

very large amount grease often discharged from slaughterhouses and 
meat products factories. Much this can eliminated intelligent re- 
covery the factory. There tendency make factory grease-traps too 
small; wherever possible they should size such that the intermittent 
discharge hot effluent cannot and emulsify the fat which normally 
floats the surface the outlet end the trap. 


the general problems resulting from grease sewage presented this 
paper, which includes methods grease determination and removal, the effect 
upon receiving waterways, and brief statement regarding methods disposal. 
sewage treatment works, grease disposal always problem the operator. 
Its magnitude depends upon the size the plant, the quantity and character- 
istics the grease the sewage, and the facilities the treatment plant. 
Brief comments the handling and disposal grease sewage treatment 
works are offered. 

Although skimming grease from sedimentation tanks reasonably suc- 
cessful, the skimming equipment becomes coated with grease, unsightly, 
and there tendency produce odors. The equipment commonly 
cleaned scraping. larger installations special cleaning equipment using 
combination steam and cleaning compound may used effectively and 
with saving labor. Frequent cleaning will keep the equipment satis- 
factory condition. 

After grease has been removed from skimming tanks sedimentation tanks, 
must transported the point disposal. small plants this can 
done hand. large plants pumping probably the most satisfactory 
method. The difficulties involved result from the characteristic grease 
stick whatever comes contact with, and from the presence floating 
trash which comes off with the grease. This combination defies screening, 
tends clog pipes and pump valves, and forms stiff mat when allowed 
stand. 

Although these conditions cannot eliminated, number things can 
done minimize the difficulties. Trash can reduced good screening 
raw sewage. Grease can kept more less liquid moving frequently 
with sufficient water, using jets sprays necessary. 

some cases grease will coat the interior piping through which flows 
such extent increase friction head greatly and necessitate removal. 
This has occurred the Buffalo, Y., plant and there the application 
steam inside the sludge line has kept the line open. has been necessary 
apply steam intervals from four six months. has been introduced 


Asst. Engr., Greeley Hansen, Buffalo, 
Received the Secretary July 30, 1942. 
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through injector with water that the grease actually melted and 
removed hot water. Temperatures approximately 180° have been 
required. 

The final disposal grease will vary with the size the treatment plant 
and other conditions. the small plant where grease handled manually 
can buried. the moderate size plants sometimes piped low areas 
where drained and burned. common practice pump with sludge 
digestion tanks. This has the advantage reducing the quantity, some 
types will digest, and also adding the amount gas produced. Some grease, 
however, does not digest, least very slow digest, which case 
accumulates the top the digestion tank. the accumulation becomes 
several feet thick, tends form stiff mass, and takes volume the 
tank needed for digestion. This mass floating material consists trash 
with varying proportions grease and frequently called scum. 
excessive depth, may become necessary remove it. the treatment plant 
employs incineration for sludge disposal, the obvious procedure burn the 
grease. 

There are least two problems involved this method disposal. One 
that.of handling the grease this condition and the other the control 
temperature the incinerator. Buffalo the scum grease the sludge 
tanks softened with spray supernatant liquor until will flow and can 
pumped. then drawn off and mixed with sludge for filtering and 
incineration. the incinerator controlled regulating the 
the proportion scum and sludge. 

general, problems handling and disposing grease must solved 
each sewage treatment plant according the facilities the plant and the 
characteristics the grease obtained the plant. Some procedures that 
have been used with varying degrees success have been indicated. Grease 
sewage, from the standpoints determination, removal, handling, and 
disposal, constitutes one the problems sewage treatment requiring con- 
tinued study and research. 
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AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 


1852 


DISCUSSIONS 


TECHNIQUE DETERMINING SHEARING 
STRENGTH SOILS 


PROGRESS REPORT SPECIAL COMMITTEE 
THE 

SOIL MECHANICS AND FOUNDATIONS DIVISION THE 

TECHNIQUE SOIL TESTS 


Discussion 


presented concisely the essential features the several devices now common 
use for the testing soils shear and has noted many the problems still 
unsolved their uses and the applications shear test results. The sub- 
ject matter the committee report divided into three parts: (1) Equipment 
used for shear testing; (2) testing procedure including sample preparation; and 
(3) methods interpreting and applying test data. The writer proposes 
follow this division his discussion. 

Equipment Used for Shear committee has advisedly empha- 
sized the merits and limitations each type shear testing device the 
basis the distribution stresses within the specimens being tested. Con- 
sideration stress distribution within shear test specimens can well sepa- 
rated into two substantially independent parts: 


(1) The distribution stresses imposed the entire soil mass the 
test specimen the mechanical features the testing device; and 

(2) The distribution stresses within the soil mass between its solid, 
liquid, and gaseous phases. 


The second these two parts involves study the true “intergranular 
referred the committee and primarily problem testing procedure 
and interpretation test results. The first part specifically concerned with 


Report was published February, 1942, Proceedings. Discussion this Report has 
appeared Proceedings, follows: February, 1942, Glennon Gilboy, Assoc. Am. Soc. E.; April, 
1942, Cloyd Beerup, Assoc. Am. Soc. E.; and June, 1942, Jacob Feld, Am. Soe. 


Prof., Soil Mechanics, School Civ. Eng., Purdue Univ., West Lafayette, Ind. 
the Secretary June 1942. 
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the testing devices, which the committee has separated generically into trans- 
latory shear and cylindrical triaxial compression types. important 
feature the distribution stresses each these types has not been men- 
tioned the committee. 

the first place, translatory shear devices attempt measure directly 
the tangential shearing stresses acting along one two planes within the test 
specimen. Shearing stresses result from the shearing strains induced trans- 
latory motion between adjacent rigid parts the device holding the specimen. 
Shearing strains are angular distortions material within zone shear 
finite thickness. For given translatory motion between two parts speci- 
men, shearing strains will vary directly with the thickness the zone shear. 
round rectangular box-type translatory devices the thickness the zone 
shear controlled the ends the specimen the space between the 
adjacent rigid translatory parts the device and the interior the specimen 
its over-all dimensions. the specimen thicker than the opening be- 
tween the rigid parts the device holding it, the thickness the zone shear 
will greater the interior than the ends, and shearing strains and stresses 
will nonuniform over the cross section the direction translation. There- 
fore progressive failure, starting points minimum thickness the zone 
shear, will occur all except perfectly plastic materials. This characteristic 
Soc. E., 1938. 

reference cylindrical triaxial compression testing devices, the com- 
mittee has emphasized the unknown effect internal stresses the lateral 
restraint friction the rigid end-bearing surfaces through which axial load 
applied the specimen. Lateral end restraint the common condition 
all tension and compression tests all engineering materials. Within the 
range elastic deformations, Saint-Venant’s principle states that the 
forces acting small portion the surface elastic body are replaced 
another statically equivalent system forces acting the same portion 
the surface, this redistribution loading produces substantial changes 
the stress locally but has negligible effect the stresses distances which 
are large comparison with the linear dimensions the surface which the 
forces are Therefore, specimen relatively long compared 
with its end diameter, lateral end restraint together with uneven application 
the axial loads will cause negligible deviations from uniform stress its 
central portion, where failure usually occurs. The length-to-diameter ratio 
compression specimens limited approximately four the tendency 
toward buckling slender specimens. Thus, the distance from the ends 
the zone failure limited ‘‘in comparison with the linear dimensions the 
surface which the forces are This limitation offset some 
extent more rapid redistribution toward statically equivalent uniform stress 
during semi-plastic deformations, because large strains cause relatively small 


Cambridge, Mass. 


“Shearing Properties Ottawa Sand,” Taylor and Leps, Proceedings, Soils and 
Foundations Conference the Engr. Dept., Boston, Mass., 1938. 


Elasticity,” Timoshenko, New York, Y., 1934, 31. 
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changes stress. For these reasons the writer cannot share the concern the 
committee about nonuniform distribution stresses cylindrical compression 
specimens which have length-to-diameter ratio great three. 

Finally, regard testing devices, the committee deplores loading equip- 
ment which amateurish appearance and consist- 
ing cans, wires, bicycle wheels, hanging weights, levers, pulleys, and 
wooden framework Would not better judge testing equip- 
ment the bases its usability and the accuracy and reliability the 
results produces, and not its appearance? 

Testing testing procedure the committee’s 
types loading equipment and the quick,” and 
shear tests cohesive soils. The writer believes that testing procedure 
fully important the testing device used. tests cohesive soils the 
results are almost completely dependent the testing procedure and the con- 
dition the specimens being tested. Such tests can divided into two gen- 
eral groups: 


(1) Tests ‘‘undisturbed” specimens determine the strength the 
soil exists nature; and 

(2) Tests remolded and compacted specimens determine the strength 
the soil after placing earth dam, levee, embankment, foundation fill. 


The results unconfined compression tests, triaxial compression 
tests, and translatory shear tests under zero normal pressure natural co- 
hesive soils are seriously affected sample disturbance. Relatively slight 
sample disturbance may reduce the strength obtained from tests less than 
half the actual strength the ground. the other hand, sample dis- 
turbance prior and triaxial compression 
tests and translatory shear tests under normal pressures one more tons 
per square foot usually will result the test being performed specimen 
with void ratio considerably smaller than the natural void ratio the soil 
inthe ground. The result test disturbed specimen density condi- 
tion quite different from any which will attained the soil under field 
loading conditions. Although the effects shearing strength sample dis- 
turbance and decreased void ratio increased density due consolidation 
following sample disturbance are compensating, any approximation between 
test and natural strengths accidental. 

Tests compacted soils used the construction earth dams 
fills frequently have been performed the material compacted near 
the optimum water content. this condition the soil has between 80% and 
90% its voids filled with Newnam, Jr., Assoc. Am. Soc. 
E., has shown that the air voids terms total volume cohesive soil 
compacted the standard Proctor optimum density will fall between and 
Tests cohesive soils this condition are not indicative the strength 
such materials field applications which they may become saturated. 


Available Soils Airport Newnam, Jr., Proceedings, Fifth Texas 
Conference Soil Mechanics and Foundation Eng., Univ. Texas, Austin, Tex., Pt. II, February, 1942. 
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the opinion the writer such specimens should tested, not compacted, 
but nearly possible the condition which they will exist after saturation 
under the field conditions. other words, the specimens should com- 
pacted, should saturated under external stress conditions and water pres- 
sures which are approximately comparable those which will occur the 
prototype, and then tested. The results such tests performed the writer 
show smaller strengths and strength characteristics entirely different from those 
the soil initially compacted. general, most such test specimens de- 
crease volume during saturation. The only exceptions have been specimens 
compacted relatively dry condition densities greater than the standard 
Proctor optimum compaction and saturated under external normal loads less 
than 0.5 ton per ft. The maximum observed expansion the latter group 
specimens was 2.2% total volume. Saturation under prototype condi- 
tions prior testing eliminates the dangers, true imaginary, over-com- 
paction cohesive soils because strengths from tests partly saturated 
specimens not create false sense security the mind the designer. 

Interpreting and Applying Shear Test the has 
stressed the importance interpretation and application the data obtained 
from shear tests, has justifiably refrained from detailed discussion this 
subject, particularly applied cohesive soils. The results shear tests 
must eventually reduced the Coulomb equation simple modification 
it. Otherwise the complications introduced when shearing strength ap- 
plied actual stability problem become intolerable. Therefore interpreta- 
tion shear tests resolves itself into reduction the results, which are 
usually complicated the effects pore water stresses and sample dis- 
turbance, into values which can substituted some form the Coulomb 
equation with reasonable accuracy for the range probable prototype stresses. 

The committee’s report gives the impression that analyses stresses the 
prototype are vital factor the application the results tests for the 
shearing strength soil. the opinion the writer, the probabilities 
even crudely approximate analyses mass stresses soil which the actual 
stresses approach the shearing strength are remote. The solutions theory 
elasticity, which are practically the only tools now available for the approxima- 
tion soil stresses, are limited the range reasonably elastic behavior 
the soil. For most soils this range does not extend beyond 50% the soil 
strength and the deformations caused stresses corresponding the last 
50% the soil strength may easily ten times those for the first 
factory solutions for problems which the soil shearing strength the con- 
trolling factor depend knowledge the shape and position the surface 
along which shearing rupture will take place. both the shape and position 
rupture surface are known from field observations, the problem becomes 
statically determinate. stresses and deformations can then 
confined secondary considerations such the determination the points 
application external reactions the sliding soil mass. The distinction 
here parallel the basic difference between elastic analysis and limit design 
structural engineering. many the available stability solutions only 
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the shape the most probable rupture surface known and repeated statically 
determinate solutions are made until the position for the surface the given 
shape found which results minimum resistance rupture. 

The essential point applications shear test results studies safety 
against failure that analyses mass stresses and field observations mass 
stresses prototypes are little value unless the soil considered, actually 
is, the point failure. Field observations the mode failure, the 
actual shapes the surfaces along which shearing rupture takes place na- 
ture, are the essence such applications. This qualification the committee’s 
report refers exclusively mass stresses, opposed intergranular and pore- 
water stresses. Field observations pore-water stresses prototypes are 
tremendously important but are beyond the scope the points the writer 
wishes make. 

The significant feature the Swedish circle method for the analysis 
stability slopes, which was mentioned the committee, that was based 
on, and has been confirmed by, field observations the shapes actual sliding 
surfaces slopes. This feature completely independent shear test re- 
sults the way which the test results are applied statical analysis. 
Any improvement testing equipment procedure, interpretation 
test results, improvement the solutions obtained this method. The 
method itself can improved only possibly more precise approximations 
actual sliding surfaces. 

The writer full agreement with the conclusions the committee. 
particular, conclusion (e), with the concluding remarks concerning soil stresses 
interpreted mean conditions failure prototypes, deserves careful thought 
every engineer engaged soils studies. 
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DISCUSSIONS 


VISCOSITY AND SURFACE TENSION EFFECTS 
V-NOTCH WEIR COEFFICIENTS 


Discussion 


AND RALPH POWELL 


Soc. indeed refreshing read paper such that presented 
Professor Lenz. With this paper, the explored realm fluid mechanics and 
dimensional analysis has been widened, and the author congratulated 
the capable manner which this has been accomplished. 

Since the publication Reynolds’ outstanding number papers 
have been presented the flow pipes which attempt has been made 
obtain wider correlation and make the results more general application. 
Although further advance has been, and will be, made this connection, the 
paper stands out the pioneering attempt make the results 
general use the engineering profession. The solution the practical pipe 
problems now commonly made means these generalized pipe friction 
coefficients. 

the generalized method became more commonly accepted desirable for 
the solution practical pipe-flow problems, attention turned obtaining 
generalized coefficients for the various measuring devices. Excellent general- 
izations the coefficients for venturi tubes" and orifice meters!? have been 
presented. The current paper extends the same type reasoning the 
coefficients the V-notch weir and covers the widest range conditions con- 


paper Arno Lenz, Assoc. Am. Soc. E., was published March, 1942, 


Proceedings. this has appeared Proceedings, follows: June, 1942, 
Kalinske, Assoc. Am. Soc. 


Prof., Mechanics and Louisiana State Univ., Baton Rouge, La. 
Asst. Prof., Eng. Mechanics, Louisiana State Univ., Baton Rouge, La. 
Received the Secretary June 29, 1942. 


Experimental Investigation the Which Determine Whether the Motion 
Water Shall Direct Sinuous, and the Law Resistance Parallel Osborne Reynolds, 
Philosophical Transactions, Royal Soc. London, Vol. 174, Pt. III, 1883. 


“Fluid Meters—Their Theory and Applications,” Research Publication, E., 4th Ed., 1937. 


Discharge for Viscous Liquids,” Tuve and Sprenkle, Instruments, 
Vol. 1933, pp. 201-206; and Vol. 1935, pp. 225, and 232-234. 
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tained any the papers weirs which have been published date. 
Although certain refinements and additions may follow, the paper must 
considered outstanding assisting the extension the generalized 
reasoning. 

Eq. does not follow obviously from Eqs. and might expected 
from cursory reading the explanatory paragraph. and are 
those the projections the and 
C-W planes line the plane repre- 
sented Eq. 13. Such plane illus- 
particular line this plane. For given 
angle notch, the experimental data plot 
series parallel lines the same 
plane, one line for each fluid. 

Since the exponents, and Eqs. 
and not refer the same things 
those Eq. 13, confusion might have 
been avoided using letters 
tions the projections the lines differ- 
ing from those relating the plane. 

would increase the value the 
paper Professor Lenz would prepare 
three-dimensional figure showing the plane which the line (Eq. 13) appears 
and showing the projections this line that give Eqs. and 12. The values 
and from 11, 12, and should indicated, well the three 
values for these equations. This could done for one angle the notch 
—say 45°. least two experimental curves must drawn. 

would appear that undue emphasis has been placed upon the failure 
Eq. apply for the full range the values and Experience has 
taught engineers that certain factors predominate over certain range, whereas 
others may govern over another. Thus, there are three definite ranges for the 
coefficient drag for the flow past cylinder, for the value for the flow 
through pipes. There are also three ranges for the value for the flow 
through pipe orifice, but there are only two ranges for the value for the 
flow through venturi tube. There are three ranges for the coefficient for 
the V-notch weir—namely, that for which the flow primarily governed the 
surface tension the liquid, that for which primarily governed the vis- 
cosity the liquid, and then that for which independent the surface 
tension and the viscous effects. seems rather unreasonable expect one 
empirical equation (even though derived rationally) cover the various 
ranges without loss accuracy. should suffice simply state the applicable 
range without the need for feeling apology. 


Emery Jun. Am. Soc. the third paragraph the 
“Synopsis” the author makes the statement: “It was found that the weir 


Hydrographer, Utah Copper Co., Salt Lake City, Utah. 
the Secretary July 1942. 
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coefficient increases with both viscosity and surface tension, and that de- 
creases with increase head and angle increase viscosity 
and surface tension does increase the weir coefficient, naturally would follow 
that increase head and angle V-notch would decrease the weir 
cient, for the head and angle V-notch increases the effect viscosity and 
surface tension decreases, and hence the weir coefficient decreases. true 
that the weir coefficient will increase with viscosity and surface 
statement just the opposite from that which the writer would expect, for 
viscosity the reciprocal and not measure fluidity. the author’s 
statement holds true for highly viscous liquids, then, for given head, glucose, 
honey, molasses would flow faster than water. The writer noticed for the 
oils used (with the exception oil that given temperature oil with 
higher viscosity had lower surface tension (see Fig. 19). This might indicate 
that the effect decreasing surface tension was overcome the effect due 
the increase but this contrary the author’s statement 
ysis Data: Development General Equation’’): gen- 
erally supposed that viscosity plays much more important than surface 


tension determining 


Kinematic Viscosity, per Sec 


Fic. 19.—RELATIONSHIP BETWEEN VISCOSITY AND SURFACE TENSION FOR THE GIVEN OILS 


the same oil used and the temperature varies (Fig. 12), the surface 
tension increases with the viscosity. This relationship just the reverse 
that noticed for different oils the same temperature. These two opposite 
relationships, together with the author’s experimental results, bear out the 
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author’s statement that plays much more important réle than 
surface and prove the falsity the writer’s suggestion herein. 

What, then, the answer the apparent reversal natural laws? 
Fig. 13, for values Reynolds’ number greater than certain critical value 
(for example, 845 for 28° V-notch) the plotted points fall the author’s 
straight line; but for values Reynolds’ number less than this critical value 
there rapid drawing away from the straight line, apparently forming 
parabolic curve. Thus, for the lower values Reynolds’ number—or, 
other words, the viscosity increases, becomes more important—the co- 
efficient (C) will decrease instead increase. some value Reynolds’ 
number the tangent the true curve becomes zero approximately 400 
for 28° V-notch). For values Reynolds’ number less than this the slope 
the tangent positive and for values Reynolds’ number greater than this 
the slope negative. Thus, for the highly viscous liquids the coefficient will 
decrease the viscosity increases, which would expected. 

1643 Torricelli his Motu Gravium Naturaliter Accelerato enunciated 
the law that the base all weir 


“The velocity fluid passing through orifice the side 
reservoir the same that which would acquired heavy body 
falling freely through the vertical height measured from the surface the 
fluid the reservoir the center the 


the use Torricelli’s law and calculus, Eqs. and derived. This 
derivation simple and may found most hydraulic textbooks. The 
author arrives the same formula using method developed Edgar Buck- 
1915—namely, the This method dimensional 
analysis, because its value has been proved model experiments hydraulic 
laboratories, now often used. Because this increase importance, 
detailed investigation the development this theorem may order. 

the physical quantities, Gi, Go, Gn, such length, force, 
velocity, density, etc., suffice describe any physical phenomenon com- 
pletely, then they are mutually dependent some known unknown manner. 


This mutual dependence may stated mathematically function 
independent variables: 


other words, the mutual dependence the quantities can expressed 


polynomial independent variables. Any polynomial these vari- 
ables may written: 


Experiments and the Forms Empirical Buckingham, Transactions, 
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or, since one the terms 


equal unity: 


his the Buckingham makes the statement: 


all the terms any physical equation must have the same 
dimensions and the coefficients have dimensions because they are 
merely ratios the dimensionless coefficients 


This merely another way expressing one the first principles physical 
mathematics—namely, dogs, cats, and pigs cannot physically added to- 
gether; they must restated terms some common attribute. other 
words, each term Eq. Eq. must have the same dimensions. Upon 
the validity this statement depends the proof the theorem and hence the 
usefulness dimensional analysis. Since each term Eq. Eq. must 
have the same dimensions, each term Eq. will dimensionless. 

Every physical quantity composed two distinct characteristics— 
magnitude and dimension. Magnitude dimensionless and may repre- 
sented numerical constant whereas the dimension the quantity 
determined certain fundamental units such mass, length, time, 
some power which known can determined soon the quantity itself 
known. mechanics the number these fundamental units three— 
mass, length, problems dealing with thermal and electromagnetic 
well mechanical quantities, the number fundamental units will not 
exceed five. The dimension quantity may represented, therefore, 
and and are symbols for these units. Thus, 


(24) 
Assuming for simplicity that mechanical quantities only are involved, Eq. 
becomes, 
noted previously, Eq. dimensionless; hence, 

which represents the expression which has the same dimension unity. 
25, therefore, becomes, 

and from Kq. 26: 
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The unknowns Eqs. are The remainder will known 
any particular problem. 


Again referring Buckingham’s paper 


“Tf the number fundamental units needed, there always, 
among all the [quantities least one set which are independent 
one another and might themselves used fundamental units 
paid attention the question representing and preserving the 
fundamental units satisfactory primary 


used dimensions and the quantities that are left, 
then from the foregoing 


which the brackets are used symbol “the dimension of.” Sub- 
stituting the dimensions Eq. Eq. 29a, then: 


Letting equal the ratio the terms both sides the equal sign 
29: 


Also, 


which the brackets are again used symbol dimension 
demonstrated Eqs. 28, possible write from Eq. 


As, 


Eqs. 28,-the unknowns are and, since there are many 


equations there are unknowns, Eqs. can solved. Eqs. are 
completely solved and Eq. becomes: 


The only unknowns now are the value the exponents for each 
and each term. Since the Eq. 32a also one the terms Eq. 
divided itself, also may included the summation, and therefore 


(30a) 
(31) 
” 
Vv. 
"7" 
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Eq. 32a can written: 


which the same Eq. the paper. 


fundamental units: 


_ 


and from Eq. 


When the exponents Eqs. are determined: 


= 


= 


When these values are substituted Eq. 32c, Eq. the result, from which 
the author develops Eq. 


get Eq. the author replaces ‘‘V each the products 
sional equivalent, the particular function Eq. were known, 
this substitution would not entirely correct because the true relationship 
between and Eq. 6a; but, since the particular function not known, 
this replacement permissible. least dimensionally correct. more 
direct approach Eq. would choose and the three units 
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and 


@ 


Substituting these values Eq. 32c and solving for demonstrated 
the author, will yield Eq. 10. 

The writer takes exception the author’s statement Data: 
Development General Equation”): functional relationship, be- 
tween and for given weir must obtained from experimental 
This statement should be: functional relationship present must 
obtained from experimental 

There theoretical solution every phenomenon nature. The 
trouble lies with Man. the relationship simple that Man has 
not looked that direction yet, complex that present exceeds 
the capacity the mind solve it. The fact still remains, however, that 
there theoretical solution the problem. Eq. will the form 


Bz, B,, 


letting the value the function and approach zero, 


Bz, B,, 
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for simplicity all terms but one are assumed either zero unimportant: 


(see Eq. 13). the effect Weber’s number (W) in- 


viously, these equations can used only between certain limits. Beyond 
these limits Eq. must used, for the experimental data when plotted 
logarithmic paper seem pull away from the author’s straight line (see Fig. 
13) and form curve resembling parabola. Just what the particular 
tion Eq. probably will revealed approximated further experi- 
mental investigation. 


(compare Eq. 12). noted pre- 


desired, not only lays basis for the measurement viscous liquids 
means notch weirs, but also indicates the way which the flow less 
viscous fluids, like water, may affected changes temperature. For 
many years the writer has felt that, instead trying study the effect 
temperature upon the coefficient notched weirs which measure water, 
great refinement measurement, would more satisfactory find the 
influence viscosity the use more viscous liquids where the effect would 
greatly increased and therefore more readily measured. Mr. Lenz now 
shows that surface tension also has important effect and gives quantitative 
results. 

However, the writer feels constrained mention several things that 
personally wishes had been done differently. One minor point that 


1.5 
the ideal velocity due the head This gives 


times 
much Eq. 9a. Similarly, would follow the definition the Weber 


number given Hunter Assoc. Am. Soc. 


twice the value given Eq. 9b. However, the following comments, the 

author’s definitions have been used. The writer also very much wishes that 

the actual numerical data had been included the paper. Desiring test 

the possibility representing the author’s experimental results formuls 

different form, obtained copy the original from the University 
Associate Prof. Mechanics, Ohio State Univ., Columbus, Ohio. 


Received the Secretary July 1942. 


Mechanics for Hydraulic Engineers,” Hunter Rouse, McGraw-Hill Book Co., 
New York, Y., 1938, 23. 


4“The Effect Viscosity and Surface Tension Upon V-Notch Weir Coefficients,” Arno Lent 
thesis presented the University Wisconsin 1940, partial fulfilment the requirements 
degree Doctor Philosophy. 
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Wisconsin, and was surprised find that even this did not give the original 
measurements. 

1910 showed that the coefficient for water flowing 
through V-notch weir could very accurately represented over wide 
range heads formula the form 


Scaling values from the curves the paper and plotting against indicated 


that all the author’s values for 90° notch above the minimum heads shown 
Fig. fit this formula with maximum discrepancy about 0.5%. This 
using 0.552 for all four liquids and the value for each liquid that gives 
the best fit. 

The writer wished investigate formula the form 


By B, 
(43) 


“— Os 


would require that, Eq. 43, and was found that the best 
fit for formula this type and the author’s data 90° notch weirs was 
given 

0.435 0.201 


C= 0.552 + Ries + wie (44) 


All the points shown Fig. (except those below the minimum head indicated) 
were computed this formula and the discrepancies tabulated. The algebraic 
mean the discrepancies was 0.0003 and the arithmetic mean, 0.0045. This 
about the same mean discrepancy the author’s formula gives, and 
seems the writer that Eq. somewhat simpler than Eq. least 
has the advantage separating the viscosity and surface tension 
interesting note that, according Eq. for water, the effect surface 
tension nearly three and one-half times great that viscosity. 
expected, for the oils the viscosity had more effect than the surface tension 
—nearly twice much the case oil 70° 

should also noted that, whereas the range viscosities studied 
quite large, the variation surface tension relatively small. The property 
that effective here what might called the surface 


where the specific weight. The expression area, say feet 
Squared, and 

kinematic surface tension 
When one consults handbook for possible liquids that would give wider 
range kinematic surface tensions than those used the author, the search 
Engineering (London), Vol. 90, October 28, 1910, 598. 
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not very rewarding. Mercury has surface tension some seven times 
great water, but weighs 13.6 times much, that its kinematic surface 


tension less—about the same that the oils used the author. Water 
seems stand the top the list (except that some water solutions sub- 
stances like CaCl, are somewhat higher). Carbon tetrachloride has 
matic surface tension about one fourth that water and would probably 
the best liquid compare with water get more accurate measure 
the effect surface tension. Another liquid that might used water 
which small amount saponin and still smaller amount oleic acid have 
been This will not give nearly wide range carbon tetrachloride 
(in fact little less than the oil used the author), but would have the 
advantage having the same viscosity water alone, and therefore varying 
only one factor time. 

The writer has not taken time develop similar Eq. for 
other notch angles, and study the variation the coefficients with 
would only after such study that the relative merits formula the 
form Eq. could compared with the author’s Eqs. and 14, and course 
there may some other form that better than any these. 

The author has given little attention the interesting question the 
variation with and when either both the latter are below what 
calls the minimum values. Water the best liquid with which study 
low Weber numbers, and Barr’s experiments ran down Weber number 
about 54. These results, put the units this paper, and with the tempera- les 
ture assumed 46.4° (the mean the range given), are plotted Fig. 20. fit 
(Some these data were used the elsewhere; but there the coefficient 

the 

(44) 
Minimum 
ind 
the 
400 2000 4000 6000 
Weber Number, 

bot 
hea 
unfortunately was plotted against the Reynolds number, procedure which tha 
this paper shows have been quite beside the point.) Curve Fig. nes 
the most striking, and interesting that this semilogarithmic plotting 


Study Foam Formation Aqueous Virgil Leland Hansley, 
doctor’s thesis, Ohio State Univ., Columbus, Ohio, 1928, 
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perfectly straight from 166. The similarity between this 
curve and that orifice plotted against also worthy note. 
The fact that curves and are different may partly due the fact 
that the writer’s assumption the same temperature for the three cases 
introduces appreciable error, but more probable that the slight difference 
the three weirs was responsible. so, before the V-notch weir can ever 
become instrument for absolute measurement (dependable without calibra- 
tion), the size plate, thickness edge, and angle bevel must standard- 
ized more definitely. original data included record the tem- 
perature, shown the sample that was printed. possible that these 
records are still existence, that and could even yet found for each 
his runs?) 

For water the temperature assumed, and were computed for 
number different heads, and from these was computed and 14, 
and 44, and the results plotted Fig. the author’s edge was in., 
would expected that his results this temperature would lie between 
curves and but nearer curve (Curve the one from which Barr’s 
values, plotted Fig. 15, are taken.) make the discrepancy more than 
would require changing the minimum from 300 about 600. 
matter fact, any formula the kind proposed the author (or the writer) 
could not expected fit beyond the point inflection the curve, and for 
least low the lowest the curve—that is, 340. This the curve 


fitted Mr. and changing his formula into one similar Eq. 
gives 


0.435 0.102 


This fits curve with maximum discrepancy about 0.04%. There are 
other pairs numerical coefficients, course, which would fit equally well— 
any that would reduce 


0.0157 
(47) 


Undoubtedly similar formulas could deduced fit curves and 
the point inflection with not much greater discrepancies. 

That the author found, for very large values and that becomes 
independent either not surprising view the analogous facts regarding 
the coefficients for orifices, nozzles, and pipes (and, the writer believes, for 
channels well); but somewhat surprising that the limit reached for 
both and precisely the same value Also, the fact that high 
heads all Barr’s weirs gave values less than 0.5886 would seem indicate 
that this limiting value depends the weir thickness (as well its sharp- 
ness and the smoothness and size the upstream face). 

Finally, physical chemists have discovered that certain solutions, such 
500 ppm peptone water 100 ppm saponin water, have what 
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This quite different from ordinary bulk vis- 
and seems due the formation plastic surface film. not 
likely that the liquids used the author had this property, but would seem 
probable that, liquid had it, its flow over notched weir would affected 
thereby. Therefore, experiments were made saponin and oleic 
acid solution suggested herein, this matter should investigated. 


Study Foam Formation Aqueous Virgil Leland Hansley, unpublished 
doctor’s thesis, Ohio State Univ., Columbus, Ohio, 1928, pp. 


lat 
has 
gra 
ter 
tha 
exe 
bac 
nat 
larg 
and 
floo. 
bott 
abo 
Dana 
ton, 
i 


AMERICAN SOCIETY CIVIL ENGINEERS 


Founded November 1852 


DISCUSSIONS 


UNUSUAL EVENTS AND THEIR RELATION 
FEDERAL WATER POLICIES 


Discussion 


oversimplification, the author has iso- 
lated the germs revolution that taking place the policy the United 
States with respect water resources. The combination 
physical, one economic, and one altruistic—is producing (in fact 
has already produced) complete reorientation the water-utilization pro- 
gram. one can challenge the fact that the decade was charac- 
terized unusually large and destructive floods. The author’s statement 
that similar climatic and hydrologic conditions have “probably” not been 

_exceeded during the past century accepted. 

However, for proper evaluation this condition one the combination 
shaping national policy, the reader requires fuller historical 
background. The impact unusual climatic and hydrologic conditions upon 
national life depends upon the character national life the time their 
occurrence. illustrate, science reveals that there sound—there are 
only sound waves—if there ear intercept the sound waves. 
large degree with the occurrence unusual floods. Population density, 
economic wealth, and encroachment structures upon the natural flood plains 
and channels rivers are determining factors appraising flood effects and 
therefore the wisdom unwisdom any given flood-control policy. 

When levee construction along the Mississippi River began New Orleans, 
La., 1717, small earthen dikes from high were all that 
sary afford protection from floods. that time this appeared wise 
flood policy, although future years proved that these same levees created 
bottleneck that accentuated flood dangers; and, other sections along the river 
above New Orleans built dikes, robbing the river its flood plain, the race for 


paper Hoyt, Am. Soc. E., was published February, 1942, Proceedings. 
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Dana Wood, James Sweet, and Campbell. 
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higher and higher levees became flood-creating flood-preventing device, 
robbing Peter pay Paul. 

The passage the Swamp-Land Act 1850 established another ‘‘water- 
regulation policy” and that, from the backsight today, had many unwise 
features. The author’s paper deals with land and water, and may also 
pointed out that numerous recent land-policy discussions, appearing print, 
have cast much unjust criticism the early land-settlement policies this 
country, particularly the Homestead Act. Fairness requires that the early 
policy makers the United States, dealing with land and water, credited 
with enacting laws and putting into effect policies deemed wisest them 
the light the problems their time. Even today there are lawsuits Cali- 
fornia because efforts prevent river floods that overflow and replenish 
the ground water supply the lower plains area. 

This subject national water policy calls for more extensive plotting 
points and circumstances before anyone can indicate trend line with 
confidence. old farmer once said never had any trouble getting two 
fence posts line, but when had three ran into difficulty. 

Land and water policies should founded upon scientific and engineering 
inquiry. Unfortunately, impact economic crisis,” which the 
author refers (see has been more persuasive with the legislative 
policy makers. questionable whether federal water policy shaped 
the ebb and flow unemployment can judged engineeringly sound. 
also questionable whether the use soil and water conservation Consti- 
tutional peg for economic adjustments agriculture should admitted 
the engineering profession factor shaping national water policies. 
This not say that the engineer should avoid all consideration matters 
beyond those that are technical and immediately The danger 
that any water policy may literally pulled pieces for other purpose than 
make legal and Constitutional pegs which hang variety new 
social and political philosophies. 

The current water policy enmeshed the that has 
been built up, awaiting the sign when jobs are needed after the war ends. 
Many these projects are for the purpose conserving and using water 
resources. experience gained the past decade has shown, 
soundness does not necessarily control the order development when the 
objectives sought lie outside the bounds integrated water policy. 

Emergency relief through make-work projects was made readily acceptable 
the generous use the federal treasury. Appropriations through the 
various emergency relief acts the period exceed $15,000,000,000. 
Never before has the United States accepted the lure debt and borrowing 
witnessed the decade reviewed the author’s paper 
policy. 

This decade witnessed generous use federal loans and grants states and 
municipalities for the development water resources. The long-established 
policy the Army Engineer Corps that required local contributions toward the 
cost river-improvement structures was seriously threatened the 
Flood Control Act, and was practically eliminated the 1938 Flood Control 
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Act. The complete elimination local contributions tends make flood 
control proposition. Heretofore, the willingness share the 
cost was the best index local need flood control. The principle, rather 
than the amount, local contribution important. may be, the author 
states, that the great losses during the decade from floods and 
droughts deeper impression policy makers, than ever (see 
but public acceptance was achieved largely because the 
federal Santa Claus. Here one obvious reason the rapid transition 
wide popular movement belief into active action which the 
author speaks (see Here kind inverted “dollar 
Dollar policy making carried with certain objectives and 
created certain obligations. This has long been the situation growing out 
various forms federal aid—to wit, the beneficiary surrenders something, 
usually measure local control. 

surprise, therefore, that many new federal agencies were created— 
some Congress, some Executive Order—to exercise the rights and powers 
yielded states and communities. Widespread federal planning came into 
vogue. now embraces every possible land and water use, and extends from 
the smallest mountain stream sea level. The Nation now period 
intensive federal water-policy making accompanied rapid decline state 
and local control. 

The extension federal water policy does not stop planning. being 
projected into the fields operation and management, prompted, times, 
political rather than engineering considerations. Again the 
comes into play. The licensing policy the Federal Water Power Act 
established 1920 has been abandoned for direct federal construction water- 
power plants and federal marketing the energy generated. This latter, 
marketing phase, still much muddled state. The political desire for 
public power has been equal greater influence shaping present water- 
power policies than any circumstance the flood 
which the author refers. 

“Multiple-use projects” has become slogan which has political well 
engineering significance. This multiplicity uses has greatly extended the 
field federal participation and has served divest states much their 
control water resources. use” presumes ability the part 
planners make fine economic adjustments water-use values that, from 
engineering approach, can little more than wishful thinking. Take, for 
example, the question allocating cost between flood-control and water-power 
benefits. Army engineers used measure allowable flood-control expenditures 
capitalizing flood losses. The best that can said for this practice that 
served indicate the relative importance the regions damaged. The need 
for flood protection, and therefore the amount expenditures, determined 
the measure and extent the fear the endangered people feel, due the 
recurrence floods. this fear that generates the pressure the members 
Congress, state legislators, and the fathers” appropriate the money 
for flood control. The methods control are for the engineers determine. 
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With water power, the former practice was measure its economic value 
comparison with the cost steam power. The current policy seems 
allocate power only such share the capital investment will make 
government power comparison with prevailing prices private 
utility companies. This might acceptable policy its object were for 
other purpose than assure the disposal the government’s power by- 
product multiple-purpose projects. The trend the policy, however, 
undertake federal dictation retail rates for electric service, long recog- 
nized state function. 

The policy has not developed into true water-conservation 
policy such was envisioned years ago such men George Maxwell and 
the late Frederick Newell, Am. Soc. The nearest approach 
water-conservation policy recent years found the program for upstream 
development furthered the Wheeler-Case Act and the Taylor Grazing Act. 
This departure water policy, well presented the author, stands 
marked contrast the earlier flood-control and navigation-improvement 
policy which was rigidly confined the lower reaches major streams. 

The stimulus fact finding and hydrologic research, which has come about 
not only through the demand for more and better data but also through the 
generous allocation appropriation funds these ends, valuable part 
the new policy from the engineering standpoint. hoped that the 
politicians and planners eventually will guided more and more this 
research than politics and social philosophy. 

The writer concurs with the author’s concluding expectation that only 
future experience can show the soundness unsoundness the prevailing 
water policy and, therefore, the need for modifications. Only gearing the 
planning practical, and, therefore, immediate future can the needs today 
and the probable needs tomorrow balanced with the greatest measure 
justice and practicability. Such planning, maximum value, requires 
constant adjustment the light social and economic changes and, one may 
well add, the light what the author describes abnormalities climatic 
and hydrologic conditions. its very best sense politics must matter 
issues; and national water policies, regard both practice and planning, 
should based upon economic and engineering facts and should free 
such issues. 


developments and interrelations federal water policies during the decade be- 
tween 1930 and 1940 presented the author. The paper particularly 
noteworthy because presents vast amount material condensed form 
readable length. 

Perhaps somewhat more significance should placed the great flood 
1927 the Mississippi Valley influencing change federal water policies, 
and the revolution federal water policies might well dated from that time 
rather than 1930. was the great flood 1927 that really demonstrated that 


Head Engr., Hydr. Data Div., TVA, Knoxville, Tenn. 
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flood control national problem far too big for local interests even good 
times. The federal legislation that resulted 1928, following this flood, 
recognized flood control for the first time national responsibility. Although 
levees had been built for many years previously along the Mississippi River, this 
work had been done under the guise navigation rather than flood control. 
The studies the Army engineers, embodied the Jadwin plan, and later de- 
velopments showed that the ultimate flood protection the Mississippi Valley 
required something more than levees, and comprehensive reservoir plan was 
developed for the Ohio River Basin. This plan was stepping stone the flood 
control acts 1936 and 1938. 

The creation the TVA 1933 marked not only the inception mul- 
tiple-use project the federal government, the author states, but also 
created for the first time complete watershed development project with all 
phases that development under one federal agency. the author has indi- 
cated, other federal legislation has provided for cooperation between various 
governmental agencies, such the Corps Engineers and the Department 
Agriculture, but only the TVA has the complete water control large 
drainage basin, including water control both the land and the streams, 
been authorized and established under one agency. The success this method 
handling the water problems region has resulted legislation being 
introduced into Congress create other similar authorities. 

This completely unified type organization has facilitated the planning 
water-control projects, utilizing the contributions various sciences such 
hydrology, forestry, and agriculture, and appears have distinct advantages 
the development complete water-control plan for region over the more 
cumbersome procedure planning two more separate governmental 
agencies, each jealous its own particular field the culmination the final 
plan. Integral with the TVA type organization recognition the basic 
hydrologic information needed for complete water-control planning from all 
standpoints and maximum utilization the data the several sciences in- 
volved. Illustrative this may mentioned the pioneer work the TVA 
the extension the network rainfall stations throughout the Tennessee 
Valley, particularly those higher altitudes for which accurate information 
rainfall amounts and intensities previously existed. 

most beneficial result recent federal water policies has been the develop- 
ment closer cooperation between various government agencies concerned 
with water projects and the recognition the interrelationships the interests 
these several agencies specific problems. The TVA, for example, although 
organized handle completely the water problems its region has, from its 
beginning, recognized the value, and pursued policy, close cooperation with 
existing governmental agencies avoid duplication effort. Stream-gaging 
work for the TVA done the Geological Survey with which there has 
been close and most satisfactory cooperation throughout the life the 
Authority. This relationship has resulted increase the scope and im- 
portance the work the Survey the Tennessee Valley region. Where 
formerly stream-gaging stations were operated the Tennessee Valley, and the 
results published years later Water-Supply Papers, immediate problems 
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planning and construction TVA projects have resulted stream-flow 
information being obtained from increased number stations and kept 
date for the day-to-day needs the Authority. The making stream-flow 
measurements below large dam the Tennessee River the morning 
used connection with operations the afternoon quite different from 
collecting records worked and published years later. 

Another example cooperative relationship that between the Weather 
Bureau and the TVA with respect weather forecasting, which utilized 
the operation the many reservoirs the Tennessee Valley. This unique 
arrangement provides the TVA with weather forecasts twice daily oftener 
and the TVA turn forecasts river stages throughout the Tennessee Valley, the 
information both weather and rivers being given the public the form 
cooperative daily bulletin. was recognized Weather Bureau officials 
that, although the Bureau required law forecast river stages throughout 
the United States, since the TVA necessarily had much more detailed 
forecasting connection with its reservoir operations, would inefficient 
for the Weather Bureau duplicate this work. 
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AND RALPH STEWART 


Assoc. Am. Soc. writer wishes call 
attention two mistakes made the derivation the theory creep 
presented the paper Professor Marin. 

Creep Law.—For the basis his theory the author (see Eq. 7a) selects 
the following relationship between the unit stress and the unit strain 
material and may assume any integral fractional positive value above unity. 

order that stress-strain relationship could used advantage the 
derivation theory flexure should valid both for the positive and 
the negative values and s—that is, both for the state tension and the 
state compression, requirement which, course, satisfied the familiar 
formula that expresses Hooke’s. This condition, however, true with regard 
Eq. only when the exponent signifies odd integer. the other 
hand, when assumes the value even integer fraction, this equation 
under the conditions negative stress becomes inconsistent, since the left- 
hand side then negative, and the right-hand side positive when 
even integer and imaginary when fraction. 

Eq. would correct, course, from mathematical point view 
and were used the sense the numerical values stress and strain. This, 
however, was not the intention the author, and with this interpretation Eq. 
would scarcely suitable for the purpose intended. 

The error discussed herein leads inconsistencies further derivations. 
Thus, the basic differential equation (Eq. 18) relating the value the bending 
moment the beam with the value the curvature the elastic line 


incorrect when even integer fraction, since according this equation 


paper Joseph Marin, Assoc. Am. Soc. E., was published May, 1942, Pro- 
ings. 
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negative bending moment bends the beam the same sense positive 
bending moment when etc., produces imaginary curvature when 
fraction. 

The same error naturally carried into the solutions the statically inde- 
terminate structures, discussed later the paper. Take, for instance, Example 
depicted Fig.6. this figure denote the numerical 
values the end moments, since the arrows representing these quantities are 
shown their true directions. The relationship between these moments 
expressed Eq. inspection which shows that, should taken 
positive even integer, positive values and may found satisfy 
this equation—that is, satisfy the conditions the problem. 

The author has been unlucky selecting for the illustration his theory 
Table the numerical values and 7—that is, those specific 
values for which his theory still holds. Had used any other values the 
fallacy his theory would have manifested itself the absurdity the 
numerical results. 

Moment-Distribution peculiarity the stress-strain law used 
the author makes superposition deflections beams caused different 
loadings inadmissible. illustration this principle, imagine simply 
supported beam rectangular cross section loaded with concentrated 
couple the end and free from any other applied loads. The couple 
causes rotation the end through angle Imagine now new couple 
the same magnitude applied the same end the same direction. 
The second couple doubles the original moments caused the first couple, 
follows directly from statics, and doubles the original stresses anywhere 
the beam, may concluded from inspection The strains, how- 
ever, are more than doubled view the nature the stress-strain law 
assumed, because the beam fibers are working now the flatter parts the 
stress-strain curve than they were before. Larger strains result larger 
deformations, and the additional angle rotation the end produced the 
second couple alone, larger, therefore, the original angle that is, 
the presence the first state stress the beam has made the beam more 
flexible the additional loading. 

This behavior the beam has important bearing the validity the 
moment-distribution procedure, and makes evident the fact that the beam 
just considered the distribution factor the second loading smaller than 
the first loading, provided that the other members attached the same joint 
were still unstressed before the second loading. This dependence the dis- 
tribution factor (and the carry-over factor well) the state stress 
already existing the beam may demonstrated also for the fixed-ended 
beam going through the necessary algebraic derivation. The author did 
not detect this dependence because his beam had been unstressed previous 
the application the end rotation. Dependence the value distribution 
factor the state stress the beam makes the version the distribution 
procedure presented the author erroneous. 

modified moment-distribution procedure allowing for the state stress 
existing the beam conceivable theoretically but scarcely possible 


q 
Ac 


September, 1942 TURNER MECHANICS CREEP 1207 


practice. would too laborious, since each step the course distribu- 
tion would require new set values the factors each span the 
structure. 

The conclusion reached concerning the applicability the method 
moment distribution should not surprising one realizes that this method 
based the principle superposition. 

Summarizing the discussion may stated that the theory creep 
presented Professor Marin valid, from the mathematical point view, 
only the following two cases: 


When the exponent odd integer; and 
When the signs the bending moments are known beforehand all 
parts the structure. The value under these conditions may arbitrary, 
but some modification the author’s formulas will necessary the sign 
the moment not the same everywhere the structure. 


The moment-distribution procedure inadmissible even cases when the 
theory holds. 

must also kept mind that the theory presupposes constant loads and 
seems inapplicable simultaneous action loads different duration. 


lated for bringing into question the validity treating Hooke’s hypothesis 
law nature. may doubted whether Hooke was first present the 
suggestion that strain approximately proportional stress, since appears 
the Mariotte-Leibnitz theory flexure following the publication 
problem. 

that early date, caloric heat, light, electricity, and magnetism were 
regarded imponderable chemical elements. This viewpoint (that com- 
plete ignorance intermolecular forces which oppose elastic change form) 
warranted the adoption the infinitesimal point hypothesis Boscovich, 
which embodied complete disregard the temperature change incident 
the application the 

Now this temperature change which takes place upon the application 
the load furnishes measure the magnitude the load, described the 
writer first practical recognition appears the work testing 
the rolling mill where, varying the speed the machine testing, the 
operator could satisfy engineers who wanted soft steel others who wanted 
medium steel from the same metal with different speed testing. Thus, 
from the same bar, was possible develop, for soft steel, 
yield point with ultimate strength 58,500 60,000 per in., 
thus satisfying the college-trained engineer who wanted soft steel; and again, 
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increasing the speed pulling the bar, yield point 44,000 45,000 
per in., and ultimate strength 60,000 62,000 per in., could 
developed, thus satisfying the engineer who wanted medium steel! 

The energy this temperature change upon application the load depends 
the final state and independent intermediate states opposed the 
external work, which measured kinetically the mean force times the 
deformation (by the initial equilibrium the external and internal 
forces meant). 

The shift internal energy resisting the external force kinetically 
double the external work, and from this relation the physical properties the 
metal may the elastic strain thus involves not only 
elongation but temperature change well, apparent that the kinetic 
instantaneous balance equilibrium external force and internal resistance 
will shift the static balance the thermo-elastic temperature change 
dissipated radiation. further appears that the extent the specimen 
homogeneous the temperature change will proportional the applied 
force and, when this thermal limit elasticity exceeded, range imperfect 
elasticity encountered before the development the yield point the metal. 
this range imperfect elasticity developed, the cooling effect tension 
reduced until radical heating effect developed, with flow the metal. 

The shift from kinetic static equilibrium was discovered Kelvin, who 
attributed the effect viscosity instead the radiation temperature change 
incident the application the load.. Kinetically, the internal energy shift 
double the external work, but this kinetic external work but five sixths 
the static, may demonstrated many ways different formulas. 
Hooke’s hypothesis sense law nature; inexact unless the 
time element taken into consideration. 

his work the manufacture steel, found that 
applying the load the structural specimen the rate 0.1 in. per min, 
yield point 39,000 per in. was developed, whereas with the fast speed 
the machine yield point 45,000 per in. was developed. 

cover the phenomena viscosity within the thermal limit elasticity, 
one may state that, the load remains place constant intensity until 
the temperature change its application has been dissipated radiation, 
the static elongation proportional the load applied and there creep 
increased elongation. The situation altered residual strain from cold 
bending other harsh treatment. 

Residual strain presents kind elasticity its own, best illustrated 
cold bending some square rods half circle with 3-in. radius; then 
placing the horseshoe type specimen vise and sawing kerfs in. deep and 
in. apart, the distance between the ends the bent piece will decreased 
when taken from the vise. Reversing the process and sawing the kerfs 
the outside the bend, the ends will again spring will present sharper 
curvature; then sawing the kerfs both sides, the horseshoe will open out 


“Elasticity, Structure, Strength and Chemical Action Materials,” Turner, Pt. 


7 
7 
7 
7 
1 


September, 1942 TURNER MECHANICS CREEP 1209 


and decrease its curvature. This change form caused upsetting the 
internal equilibrium between the internal residual and elastic areas 

When extensometer measurements stresses concrete floors were fad 
with student investigators, often happened that they measured, according 
Hooke’s law, residual changes form which they misinterpreted elastic 
strains one and one-half one and three-quarters times great the yield 
point the metal they were investigating. (Measurements across cold bends 
kinks small rods give rise erratic results.) This kind residual creep 
has never received the attention deserves. 

Prior experiments repeated stress, the elastic limit iron 
and steel was considered little less than the yield point the metal. With 
range stress from tension compression exceeding the primitive yield 
point, failure occurred under sufficiently large number repetitions. 

Instead recognizing range imperfect elasticity beyond the natural 
elastic limit, the metaphysical range stress formulas Launhart and Wey- 
rauch was developed based upon the assumption fatigue within this natural 
limit where none was developed experiment. 

With the passage time these ranges stress formulas have been dis- 
carded. 

The writer’s experiments thermo-elastic measurement stress place 
the thermal limit elasticity approximately five eighths the yield point, 
harmony with the work Wohler. Creep increased deformation will 
occur with time under stresses beyond the thermal limit elasticity because, 
when this limit exceeded, work has been done displacing the molecular 
arrangement. 

Recalescence points melting cooling metals are recognized points 
where the molecular formation undergoes modification. For example, soft 
ingot metal the cubic system presents molecular kernels eight atomic 
spheres. Upon rolling, this structure broken down, and body-centered, nine- 
atom kernels appear. Again, drawing into fine wire, the body-centered forma- 
tion changes face-centered formation, and the increase strength under 
the law attraction the square the number atoms the kernel. 

High-Temperature higher temperatures, shift molecular 
structure may take place manner similar that occurring the recalescent 
points the metal, developing the so-called phenomena creep high 
temperature. 

natural phenomenon has been much misrepresented metaphysical 
hypotheses that elastic reaction. Volume, rather than elongation, 
the correct measure the effect axial force. Tension increases volume 
and compression reduces volume, and, although the same force causes equal 
volumetric changes whether from tension compression, the temperature 
changes differ less magnitude for tension and more for compression. 
accordance with the law volume, increase volume endothermic and 
decrease volume exothermic. This fundamental physical law destroys 
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the basis Lorentz and Einstein’s the interconvertibility 
energy and mass and proves the irrationality that characteristic energy 
matter developed mathematically from literal expression Hooke’s 
hypothesis. conservative system, the distribution energy depends 
the relation masses and irrespective the individual unit, demonstrated 
the philosophical treatment The 

The fantastic idea that one-quarter gram matter has intra-atomic 

energy equivalent that developed the engineer burning 5,000 tons 
coal metaphysical Not only has the elastic action solids been 
confused theory, but, elastic resistance gases compression, physicists 
‘have invented perfectly imaginary ideal gas which the molecules have 
volume and occupy space. Because all gases have been liquefied and 
liquids are almost incompressible, one needs only deduct the volume the 
liquefied gas from the gross volume find the net volume capable com- 
pression accordance with the laws Charles, Boyle, and Gay-Lussac. 

Professor Marin’s creep theory for beams embodies the assumption that 
plane section before bending remains plane after bending. Among the infinite 
number sections the length the beam, one only simple beam re- 
mains plane after bending. That section divides the moment area equal 
parts. applied the restrained and continuous beam, notwithstanding 
the creep correction Professor Marin, the maximum tensile stresses mid- 
span and those over the support will calculated anywhere from 15% 30% 
greater than actual and the maximum compressive stresses will approxi- 
mately correct. The deflections will much less than those computed 
the restrained and continuous beams and greater than those computed the 
ordinary formula for the simple beam. 

The most curious inconsistency the ‘‘Ut tensio sic (extension pro- 
portional force) lies the fact that, under longitudinal compression, lateral 
extension occurs without any force component the applied force account 
for it. Some have endeavored explain this misconception shearing 
stress; but the shears, beam action, can developed only principal 
stresses opposite signs, and one these principal stresses has external 
force its direction action account for its development. 

The elastic temperature creep (that due radiation temperature change 
application the load) discloses the irrationality the dynamic stress 
theory because the kinetic external work but five sixths the static; doubling 
would one and two thirds. Because the yield point before radiation 
six fifths that after radiation complete, the general effect the suddenly 
applied load developing the yield-point value but one and one-third times 

the gradually applied static load, instead twice the load, taught hand- 
books and textbooks today. 


Webster’s New International Dictionary (unabridged), Ed., 1934, 847—m 


Modern Lorentz, Ginn and Co., Boston, Mass., and New York, Y., 1927; 
Atoms,” Haas and Uhler, Van Nostrand, New York, pp. 112-113; and 

The Encyclopedia Britannica, 9th Ed., 207; also Structure, Strength, and Chemical 
Action Materials,” Turner, Pt. III, pp. 49, 51, and 141. 


Comprehensive Treatise and Theoretical Chemistry,” Joseph Mellor (Long- 
mans, Green), London, England, 1922 


n 

0 


September, 1942 TURNER MECHANICS CREEP 1211 


Initial tests newly cured reinforced-concrete beams disclose that steel 
stress times its lever under moderate loads customarily 20% 30% less than 
the applied moment. This inequality creeps out existence under repeated 
loading. Textbook authors point evidence the operation tensile 
resistance concrete; but, turning the concrete beam bottom side up, 
evidence its insignificant value concrete tension flexure appears, 
requiring more logical explanation. 

the concrete matrix crystallizes, expands, inducing tensile stresses 
the steel magnitude 4,000 6,000 per in. The reaction these 
steel stresses induces state residual compression the concrete about 
the metal and, because the reaction residual strain antipodal its elastic 
counterpart shown the cold-bent horseshoe bars, the resisting moment 
that steel stress times its lever plus residual compression the concrete 
times the same arm, with the joint effect equilibrating the external moment. 

Some deformed reinforcing bars indicate their top and bottom surface under 
their final pass through the rolls. If, 18-in. specimen, saw kerfs are 
made the vise intervals, the bar will curve slightly because the 
residual strain from rolling has been thrown out balance the reduction 
sectional area. the example the initially cured concrete beam, 
such residual strains are dissipated repeated loading. After years 
painstaking investigation, concluded that sound metal, rather 
than being fatigued brought nearer rupture alternating stresses from 
tension compression within the natural elastic limits, was improved 
elasticity, ductility, and tenacity. 

interesting example apparent creep downward motion the neutral 
plane reinforced-concrete beams was developed the National Bureau 
Standards?’ under the direction the late Richard Humphrey, Am. Soc. 
E., connection with the joint committee the Society. locating the 
position the neutral plane, Mr. Humphrey and his committee assumed that 
this locus would found where the longitudinal compression passed through 
zero, homogeneous beams, but that this locus crept downward the 
load increased. the reinforced-concrete beam, the shears between the 
tension the steel and compression the concrete are equilibrated the steel 
instead the locus the centers rotation. 

the stress mechanism homogeneous beam, the lines greatest com- 
pression above the neutral plane pair with the lines greatest tension below 
and intersect the neutral plane 45°. the concrete beam, because the 
shears are equilibrated the bottom about the steel, the states right 
shearing stress centers rotation are twisted 22.5° from their position 
the homogeneous beams. This causes the lines greatest compression 
pair with the lines greatest horizontal shear, intersecting the neutral plane 
22.5° instead 45° that the neutral plane compressed longitudinally 
proportion the applied load, causing the locus zero compression 
creep downward the load increased. Elastic deflection concrete 
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beam, therefore, increases more rapidly than the load because this com- 
pression and shortening the neutral plane. 

Failure grasp this relation has led the unscientific conception 
parabolic stress distribution the compression zone and general confusion 
and mystery regarding the co-action metal and concrete, limiting the con- 
ception the field investigation that concrete ribs functioning steel 
beams would do. 

progressive achievement resulted from the deliberations the Humphrey 
committee because they failed grasp the idea monolithic action economic 
construction and the nature and distribution not only shearing stresses 
beams but twisting stresses plates. The stress mechanism both 


and reinforced-concrete beams and plates has been developed 


the 

The theory elastic reaction treated under Hooke’s hypothesis incon- 
sistent with the statics abstract forces because lateral extension occurs 
under longitudinal force with component that force account for its 
development. thermo-elastic theory this inconsistency disappears and 
Young’s modulus may coordinated with specific heat; coefficient expan- 
sion; unit weight; atomic weight and number; strength; and melting and 
boiling temperature through (the mechanical equivalent), (the gravita- 
tional constant), the gas-volume coefficient (constant the gas- 
pressure coefficient (constant volume). 

Thermo-elastic theory, treating actual rather than imaginary substances, 
embodies nothing more complicated than common school arithmetic, which 
should render popular with the poor student who has labored with the 
tedious differential equations mathematical theory, unrelated the physical 
properties actual substances. 

The college textbook theory the nature tension and compression 
embodies the spheres exclusion developed 1650, and the foundation 
mathematical elasticity the infinitesimal point hypothesis Boscovich 
that same modernistic date and generation. 


permanent deformation highly stressed structure, which the author calls 
has been given little attention civil engineering literature. 
reading this paper the writer was impressed with the fact that contained 
valuable information the subject. 

the writer’s early practice noticed that the 20-yr-old wooden joists 
certain heavily traveled bridges had acquired pronounced sag due creep. 
For long beams continuous over number supports, the effect creep 
obviously would tend equalize the positive and negative moments after the 
manner the plastic flow overstressed structure. fact, the distinction 
between creep and plastic flow narrow and seems depend the time 
element. 

Structure, Strength and Chemical Action Materials,” Turner, Pt. III, 
plates VIII and IX. 

Engr. Bridge and Structural Design, City Los Angeles, Los Angeles, Calif. 
Received the Secretary July 13, 1942. 
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designing machinery there evidently justification for including the 
effect creep the stress analysis, but designing continuous bridges and 
building frames the neglect shear deformations may produce greater inaccu- 
racy than the neglect creep. 

One difficulty including the effect creep complicated structure 
that some points the stresses would sufficient cause creep but that 
other locations (except rare occasions) the stresses might too light 
cause creep. This would make the stresses due creep discontinuous, whereas 
the elastic stresses would continuous. 

two places the author indicates conceptions and conclusions that the 
writer considers erroneous—namely, the reference end-moment distri- 
bution character and nature, and the statement 
that should used for computing moments which include the effect creep. 


Basic Constants 
Flexure 


The basic constants beam flexure for any beam are illustrated Fig. 
which shows the elastic curves, their tangents, and their long chords for the 
three different conditions that yield the basic constants flexure. The 
illustrations can considered depicting railway curves any geometrical 
curves, and the properties may considered geometrical rather than me- 
chanical. basic constants flexure read from these figures are shown 
Fig. are the best and the minimum number constants necessary 
determine the geometry the triangles which circumscribe the elastic 
curves. 

The end-moment distribution constants are not basic but are con- 
stants. For tapering members other special conditions they become complex. 

The writer would define solution basic character one which uses the 
basic constants flexure and which does not use any form converging ap- 
Such solution shown Figs. 15. involves fewer 


l- 
n- 
its 
ta- 
as- 
| 
the 
ical 
ion 
ned 
eep. 
reep 
the 
tion 
time 


STEWART MECHANICS CREEP Discussions 


Load Kip per 
7.4=F.E.M. 


0.6 


6.85 


Maxwell’s Theorem Check 13.44 


INCLUDING THE EFFECT CREEP 
Fixed-End Moment) 


8.333 


Theorem Check 11.=11. 


1.89 3.7878 6.0606 3.03y 


(F.E.M. Fixed-End Moment) 


1214 
(a) 
2.565 
13.44 
1 
~ 
8.60 10.06 
11.0 
4 2 1 
11.0 
2 6 
a 3. 5. 6 
5 


September, 1942 STEWART MECHANICS CREEP 1215 


arithmetical operations than the author’s solution. The arithmetic simpler, 
that there are long columns figures add, and its result more accurate 
than the author’s. will compute moments any predetermined degree 
accuracy one operation. does not involve the use simultaneous 
equations. Except for the last operation, has automatic check 
theorem reciprocal deflections. The writer therefore contends 
that the author’s statements regarding methods computation are incorrect 
and that they tend perpetuate now common practice failing select the 
most effective method for analyzing large class structures. 
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promise effecting the next major advance the design and construction 
underground structures. The past three four decades have witnessed great 


improvements the manufacture cement and steel, and lesser extent 
standardization the choice aggregates, that uniformity the ordinary 
construction materials may now assured. The technique concrete pro- 
portioning has become fully recognized that, with proper engineering or- 
ganization, the failure structures due careless ignorant workmanship 
eliminated. The analysis stresses, even though many underground struc- 
tures are indeterminate, has been the subject extensive research and has 
become rationalized such extent that full dependence may placed 
the structures withstand their assumed loads. Design assumptions based 
upon anticipated earth pressures, however, still are irrational the extreme 
and contain about many unknowns must have confronted the builders 
the Egyptian pyramids. 

The results soil studies the Chicago subway, outlined Mr. Peck’s 
paper, are valuable contribution the engineering profession field that 
still relatively unexplored. They indicate that, moderate expenditure, soil 
samples carefully taken and properly analyzed the laboratory can safely 
relied upon disclose the physical characteristics the soil. The field mea- 
surements have demonstrated that upon these test results reasonably sure 
assumptions the value the lateral pressures against open cuts can 
predicated. The results also indicate that, tunneling through wide variety 
Chicago clays, with water contents ranging from 15% 60%, the structures 


Symposium was published June, 1942, Proceedings. 
Deputy Coordinator, Protective Service, Office Civilian Defense, Chicago, 
Received the Secretary June 22, 1942. 
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tend deflect the shape that will produce horizontal pressures approaching 
the full vertical loads. circular tubes that are sufficiently flexible, for ex- 
ample, follows that the bending moments commonly assumed exist are 
not found the finished structure. The results further indicate that, 
open cut, the proper spacing wales and bracing and pre-stressing 
the struts prevent deflections under load, the adjacent ground movements 
settlement can minimized. 

These results further than theory. They are practical value. They 
mean dollars saved the elimination wasteful design. They mean dollars 
saved the prevention damage from faulty construction. 

The design assumptions adopted for the Chicago subway were accordance 
with widely accepted practise. Although was felt that they were ultra- 
conservative, they were fully justified precedent. fact, the absence 
exact data the time the work was started, any departure from the current 
practise would have been 

The method most commonly used assume that the soil perfect 
fluid having arbitrary density. This density usually the result ma- 
jority opinion and has rational background. the credit the engineer- 
ing profession, usually assumed great that underground structures 
not fail when designed. further factor safety, instead adopting 
only one value for density, stresses are analyzed the basis horizontal 


structures each instance are made adequate meet the most severe these 
conditions. apparent that these assumptions disregard cohesion and other 
soil properties and are independent the depth the work. The experi- 
mental results the Chicago subway indicate that both the soil properties 
and the depth are factors that affect pressures. general, the fluid theory 
may said irrational and most instances will lead uneconomical 
design. 

The chief inherent fault the fluid hypothesis that the total lateral 
pressure thus assumed increase the square the depth, and have 
value modified somewhat whatever angle internal friction assumed. 
fact, however, the test results show that total pressures moderate depths 
are less than those computed the fluid hypothesis, although extreme depths 
the reverse may true. Further, the Chicago clays are shown act con- 
nection with excavation work they have angle internal friction 0°. 
follows that any assumptions that neglect the significant soil properties are 
irrational and, although generally the accepted assumptions are the safe 
side, occasions may arise wherein such assumptions would lead unsafe de- 
signs. Mr. Peck’s suggestions that the unconfined compression strength 
the soil and the depth the cut used factors determining total lateral 
pressure, and that pressure diagrams developed from this data, seem 
offer logical and relatively simple approach the fixing this very essential 
element design. 

Many the other results soil tests the Chicago subway have proved 
great value establishing safe construction procedures. The rates move- 
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ment the clay tunnel headings were observed spearheads driven into 
the clay advance the headings, and from these measurements practical 
cycles operations were established. Similarly, lateral movements were de- 
termined from so-called tests. close observation the extent 
the movements and the time intervals involved, practical limits for spacing 
lateral braces were fixed and measures taken for pre-stressing the tunnel 
linings avoid undue unsafe ground settlement. 

open-cut excavations the Chicago subway, extended observations 
were made outlined detail Mr. Peck’s paper. This supplied the neces- 
sary data for the determination the safe vertical spacing the wales and 
bracing, the depth penetration the sheeting, and the magnitude pre- 
stressing the struts required prevent untoward movements the adjacent 
ground and keep settlements within practical limits. 

Although much the experimental data obtained from soil tests the 
Chicago subway were collected during the progress the work, and were thus 
too late affect the design, this branch the engineering service was great 
benefit construction procedures. further felt that the results will 
inestimable value for future designs both the city departments new public 
works and also the engineering profession whole. These tests not only 
furnish more exact knowledge the physical properties Chicago clays than 
has heretofore been available, but they also point the way more rational 
basis design and safer construction procedures. 


described Professor Housel, was designed convey 60-in. steel water 
main beneath the classification yards and main-line tracks the Michigan 
Central Railroad Company and two grade separations. There were twenty- 
two tracks grade passed under and seven tracks carried the separa- 
tions, which supported the ‘‘hump” tracks the Michigan Central classifica- 
tion yards. The volume traffic the tracks, the character the soil, and 
probable future grade separation led the decision carrying this important 
transmission main beneath these obstacles the tunnel. 

The tunnel constructed was some 1,200 long between center lines 
shafts. The inside diameter was in. with average wall thickness 
about in. Reinforcing steel was used which consisted 1-in. square bars 
the inside top and bottom quadrants the lining and round bars 
the outside the spring-line quadrants. Spacing averaged about 12-in. 
centers. Longitudinal steel consisted round bars 18-in. centers. 

The soil encountered tunnel depths was typical soft Detroit blue clay 
which known “squeezing” clay. Therefore the specifications required 
the work done with compressed air. With the pressures used there was 
generally little subsidence the ground above the tunnel, and slight sub- 
sidence grade-separation abutments was recorded. particular construc- 
tion difficulties were encountered other than the usual difficulties attendant 


Mgr. and Chf. Engr., Dept. Water Supply, Detroit, Mich. 
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compressed air tunneling soft blue clays where unremitting vigilance 
the price successful construction. 

Mining operations were conducted during the day with usual progress 
per shift. Crown planking was not used except beneath the rail- 
road tracks where some eight ten planks were carried. The concrete lining 
was placed, the night shift, within about the tunnel face that 
time, except over week ends, would more than tunnel 
for period hr. Over week ends the tunnel face was well breasted. 

During the progress the tunneling, grout was forced over the concrete 
bulkheads every other day through grout holes spaced maximum 10-ft 
centers the crown with additional grout piping the upper quarter points 
where they were believed necessary. few grout holes were also placed 
the spring line and the invert checks; but grout could placed except 
through the crown grout holes and occasionally the upper quarter points. 
When the tunnel lining operations were completed grout holes were cleaned 
out and the tunnel grouted refusal. 

the writer’s opinion that the grouting tended toward uniformity 
soil pressures during the early months. Grouting pressures did not find their 
ultimate release until some time after the tunnel was constructed. This was 
noticeable holes drilled through the tunnel lining after construction showed 
water and air under pressure. This could only have originated from the grout- 
ing there was water encountered the tunneling operations. Upon 
first making the tap outside the lining the water was plainly under pressure 
but was quickly relieved and these holes soon dried up. 

The characteristic elliptical pressure observed the later measurements 
was also observed pressure cells installed other tunnels during construc- 
tion the Springwells Unfortunately access these cells 
longer available due the use the tunnels. 

The pressure cell measurements again emphasize the necessity for speed 
soft ground tunneling operations. Further, the Central Avenue Tunnel 
relatively rigid structure and, such, essentially different from thin-walled 
structures, which, due large deflections, utilize the passive resistance the 
soil maintain uniform pressures around their circumference. 

Professor Housel has presented ingenious analysis elastic structure 
plastic clay soil which cannot help but challenge attention. should 
stated that his analysis based entirely the concept soft clay truly 
cohesive material. The fact that closes his analyses through the various 
conditions loading merits careful consideration. 


practicing engineer contrasted with the soil-research engineer, 
refreshing read paper soil mechanics that not filled with theoretical 
mathematical formulas and that can followed intelligently the engineer 
(1936), Filtration Plant,’’ Eugene Hardin, Transactions, Am. Soc. E., Vol. 101 
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who has the practical problem determining the probable load that will 
supported undergound structure such tunnel. Professor Housel 
has applied theory 10-yr experiment and has arrived workable solu- 
tion age-old problem which the structure built will plastic 
soil. This valuable contribution the profession, any one will appre- 
ciate who has had the problem determining pressures conduits installed 
the tunnel method. 

True, Anson Marston, Past-President and Hon. Am. Soc. E., and 
his associates the Iowa Engineering Experiment Station have conducted 
research over long period years pressures conduits, and the result 
that monumental has made the design certain types conduits 
science instead guesswork; but the Iowa research was concerned with con- 
duits placed open trenches top the ground and then covered with 
earth, whereas Professor Housel’s work deals with structures which the 
load comes not from earth placed against the structure man but 
readjustment the plastic earth itself when part its support has been 
removed. 

Engineers who encounter the problem tunnel pressures for the first time 
are surprised learn that tunnel structures are designed more rule 
thumb than actual design data. The reason, course, that until quite 
recently very little was known the magnitude the pressure sup- 
ported, and still less the distribution that pressure around thé periphery 
the tunnel. 

There are two distinct parts the problem pressures tunnels (as 
revealed Professor Housel’s paper): One the initial pressure, which governs 
the design the temporary lining the tunnel; and the second the maximum 
final pressure that will govern the design the permanent lining. 

Some years ago when the organization with which the writer associated 
entered the metal tunnel liner field, the writer was amazed discover that 
there was rational basis which design liner-plates and way but 
costly experience for the contractor determine how much strength his plates 
would require. The practice was (and still is) place the “burden proof” 
the contractor, who was expected use his experience other jobs 
determine what liner-plates should buy for this particular contract. 
may buy heavier plates than necessary, which case one will ever know how 
much money might have been saved science had replaced guesswork. the 
contractor optimistic, the chances are that the plates will too light and 
then faced with all kinds trouble such the danger physical injury 
workmen, putting compressed-air apparatus, getting structural reinforce- 
ment members, trying get the plate manufacturer take back few 
carloads plates and supply heavier ones overnight. 

There are, course, certain limitations the application the method 
proposed Professor Housel. the first place, applies only soils 
which the shearing resistance function true cohesion and independent 
internal friction. the language the practicing engineer, this means what 
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generally classified clay, and except for tunnels bedrock probably in- 
cludes large percentage the cases encountered practice. Then, 
stated Professor Housel, the pressure equations not permit direct solu- 
tion without the introduction simplifying conditions. air used, 
this one the simplifying conditions through the introduction equal 
pressures the top, side, and bottom the tunnel. making use measured 
shear values, necessary select what the author terms 
Fortunately, enough experience already available permit the engineer 
establish this ratio with reasonable success. 

The investigation has brought light one phenomenon which may puzzle 
many, and that the possibility decreasing load the structure for the 
first three months its life. However, that not important the engineer 
contractor responsible for the temporary lining the fact that the initial 
load the temporary lining does not increase for this period, which should 
give sufficient time for placing the permanent tunnel structure. 

The results time required reach final equilibrium check quite closely 
with the results Dean Marston’s experiments, previously 
although naturally the final loads could not expected similar. 

Given the shearing values the soil (which, course, will have 
determined from soil samples obtained the site for that purpose) and suffi- 
cient other information permit assigning reasonably accurate ‘‘overload 
the engineer now has new tool for use designing both temporary 
lining for tunnel and the permanent structure well. 

Other applications this same principle doubtless will discovered 
connection with other types structures subjected loads produced earth 
masses which the condition equilibrium has been disturbed. 

Frequently, one hears the criticism that the studies soil mechanics re- 
search men are too theoretical and abstruse. one excellent example 
practical application soil mechanics, application that can remove 
lot guesswork and save contractors and owners sizable sums money. 


CHAMBERLAIN,” Assoc. Am. Soc. the construc- 
tion the Chicago subway, the writer one the consultants followed with 
interest and profit the soil investigations outlined this Symposium. 
state that construction methods during the building these works were in- 
fluenced, great degree, the results these investigations would 
understatement. The farsightedness Messrs. Harrington and Burke 
instituting these investigations contributed greatly the safety construc- 
tion and economy. Careful examination the results will make for greater 
safety and considerable saving future construction. Design future 
construction, this class ground, may approached with much clearer 
understanding the problems involved. 

the main, noted, the material found north the Chicago River 
stiff blue clay the depths involved the construction, and the capacity 


this material arch and resist shear that the successful use the so-called 
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method may attributed. The upper stratum 
hard yellow clay must have added little except weight the overburden 
since, for years, has been broken along the full length the work during the 
laying public utilities and other structures. 

matter common knowledge that the stiff blue clays found this 
region have very strong structure, evidenced the many thousands 
feet small circular tunnels that have been built without temporary support. 
Some these tunnels have had outside diameter many 
these works the permanent lining extended almost the face the excava- 
tion all times, thus reducing minimum the amount ground open 
onetime. The time that any ground without support about the same 
that required place temporary lining. The settlement these projects 
usually confined area about twice the width the structure, and the 
amount settlement increased the softer clays. 

combat this settlement method construction known the liner- 
plate method was developed, which light steel plates were placed contact 
with, and jacked hard against, the ground soon was exposed. The 
plates were supported first the bench and then needle beam placed 
near the center the tunnel. After the predetermined length tunnel 
had been excavated and the temporary lining placed, the jacks were removed 
the concreting advanced. Grouting the permanent lining was 
daily the end the finished concrete lining. The outstanding feature 
this method was that part the overburden was without support any 
time prior the completion the permanent lining. 

large number circular liner-plate tunnels have been constructed 
clay and other materials the twenty years without any consider- 
able amount settlement, including the trunk sewer (16.5 outside di- 
ameter) the City Detroit. four the Detroit tunnels there was only 
the slight uniform settlement the pavement caused the initial movement 
the clay that had been found unavoidable Chicago. interesting 
note that, Chicago also, the settlement was confined area less 
than twice the width the structure. 

The use the liner-plates became very popular method construction 
until the desire for long traveler forms for concrete, greater daily progress, and 
the application plates large-diameter tunnels caused the introduction 
internal supporting ribs. The introduction ribs, and the elimination 
the complete 24-hr cycle operation from the start excavation the 
completion the permanent lining, changed the liner-plate method back the 
old heading and bench method. 

For number years the use steel ribs with steel-plate lagging has been 
standard Chicago. Many sewer tunnel sections the Sanitary District 
have been built with complete success. Since the ground conditions the 
‘northerly sections the subway work varied but little from those which 
prevailed greater part the Sanitary District, sufficient reason could 
advanced for not using similar type structure. The fact that the cross 
section the subway tunnels was greater than the cross section those pre- 
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viously built did present many problems, the works were built 
thickly populated part the city, and settlement had reduced 
minimum. 

was recognized that the use any method construction involving 
open heading and bench excavation depends great measure the capacity 
the ground arch and support itself for least period time necessary 
the placement the temporary lining. records actual measure- 
ments the magnitude direction the loads handled the temporary 
lining were available, was necessary depend the results observed 
work previously constructed. 

Work was started contract S-5 (Fig. what was considered good 
clay. Steel ribs with liner-plates bolted them, but set away from the clay 
provide space for packing and grouting, were used. became apparent 
the early stages this work that the ground was not being held during the 
several operations essential the completion the full rib. This fact con- 
tributed excessive settlement and caused the development weight that 
several occasions resulted severe distress the ribs. 

The aforementioned situation offered fertile field for investigation, and 
actual measurements the loads being developed were made. Careful study 
these data made possible improve methods immediately, and con- 
tract S-6, clay the same stiffness, was possible reduce settlement 
materially. The principal changes that were made contract were the 
development plate that could blocked wedged contact with the 
clay when placed, and the use wall plates the spring line. all the im- 
provements, the most important was the use the wall plates, was pos- 
sible means jacks place the full arch bearing and hold until the 
posts could placed and wedged the foot blocks. These changes became 
standard for the remaining work. 

Many the difficulties that now seem inherent work this character 
can eliminated, greatly reduced, careful study data such those 
contained this Symposium. hoped that similar investigations 
may made other classes ground. 

South the Chicago River, after careful examination available data 
having with the subsurface conditions, and after evaluation the pos- 
sible difficulties due heavy loading the large structures along the line 
the work, was decided that the clay was too soft and plastic permit the 
use the method construction that had been started the northerly 
contracts. Therefore, the shield method was adopted, and proved very 
satisfactory. 

The clay both the sections driven with shields was stiff enough permit 
line and grade maintained with the jacks and judicious use taper 
rings. This fact made possible control the grade the street the 
amount and speed the flow the plastic clay through the pockets the 
shield. The surface the street was under continuous observation level 
parties and, when the rise fall the street exceeded predetermined amounts, 
the heading was notified and corrections were made the openings the 
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shield. points was found that the clay was soft that the surface 
would not respond pressure, indicating that the use any but the shield 
method would have allowed complete collapse. 

interesting note that, for almost the full length the shield-driven 
tunnel, the volume pea gravel used fill the void caused the tail the 
shield exceeded the theoretical void. This would indicate that the action 
the shield caused rupture the clay even these soft areas. 

There one conclusion drawn Professor Terzaghi from the investiga- 
tions this work with which the writer not entirely accord. the 
matter future design tunnel structures the magnitude the ones 
shown Fig. 15, shows thin-walled structure, the design which 
states based, part, the strength the clay through which the tunnel 
driven. were possible state that the ground through which any con- 
siderable length tunnel driven will uniform, and that the future 
there will disturbance the ground adjacent the structure, the 
writer would willing explore the possibility such change with more 
enthusiasm. Consideration must also given the fact that concrete and 
other lining material may not measure laboratory standards. The 
personal equation which makes for perfection workmanship must also 
kept mind all times. 
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scholarly paper great practical interest, which functions complex 
variable are used advantage. Since the analytical methods adopted will 
valuable for general reference, appropriate indicate that much 
the qualifying material preceding Eqs. 107, 108, and 109 follows directly from 
the specification (from pure mathematics) conditions that are necessary and 
this specification are the familiar Cauchy-Riemann equations which are neces- 
sary verify that Eqs. 107, 108, and 109 satisfy basic Eqs. 104 and 105 
when 

connection with Eqs. and 49, the writer was pleasantly surprised 
find that the 1938 edition the excellent tables Jahnke and Emde cur- 
rently However, believed that tables the so-called 
and type Bessel functions are more commonly used the United 
States than those giving Hankel functions the order zero for complex 
values the argument. Consequently, may useful reference 


give the requisite transformation relation which developed the standard 
treatises 26, 


paper Westergaard, Am. Soc. E., was published April, 1942, Pro- 
ceedings. 


Chf. Design Engr., TVA, Knoxville, Tenn. 
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“ker” and “kei” functions, Eqs. and become: 


(1) 


Modern tables these with convenient asymptotic expan- 
for the larger values the argument, are readily available. 
thoroughgoing coordinate treatment the general theory Bessel functions 
best obtained through the medium complex variable theory and preferably 
introduced way the Laplacian equation.?’ 

The conversion the Laplacian operator cylindrical coordinates, Eq. 
can obtained, course, elementary methods; but, since the engineer 
assigned analytical work has frequent occasion effect transformation 
many similar types differential invariants from Cartesian orthogonal 
coordinates, may stated passing that such conversions are 
obtained with striking ease and directness means the methods tensor 
analysis. For example, the case Eq. The line element is—- 


and the non-zero components the metric tensor are— 


Using notation, the tensor expression for the Laplacian 
operator, valid any system coordinates, space any number dimen- 
sions, 


The only non-zero value the Christoffel three-index symbol the fore- 
going case {22,1} and the conversion immediate: 


The author states that the coefficient defined Eqs. and isa 
logarithmic potential, and mentions electrical potential analogies that exist 
connection with the derivation Eq. 65. familiar and perhaps closer 
analogy occurs the theory elasticity which both direct and logarithmic 
potentials are used calculating the displacements produced any point 


Ed., pp. 55-64. 
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solid body the application pressure distributed over area the 
bounding 


Dana Am. Soc. determination the stresses 
plate the vicinity transverse concentrated load one the most 
troublesome problems plate analysis. The ordinary theory bending 
plates not valid for stresses point load, and necessary correct the 
results accordance with more exact theory take into consideration the 
fact that most practical applications the load actually distributed over 
small area. variety special methods have been used the past obtain 
solutions for limited number cases. The author has succeeded develop- 
ing original and powerful treatment the problem which possesses great 
generality. 

The sets place coefficients and area coefficients given the paper provide 
solutions for many useful cases. Data for simply supported and clamped 
rectangular plates finite length were not given the paper, and may 
interest derive the place coefficients for these cases. 

For rectangular plate with simply supported edges and 
and with load uniformly distributed over small circular area 
the load given are: 


and 


which 


cosh? 


Substituting Eq. 153 Eqs. 152 and rearranging: 


m=1 


“The Mathematical Theory Elasticity,” Love, Cambridge Univ. Press, 1934, 4th Ed., 


p. 1 
Applied Mechanics, The Univ. Texas, Austin, Tex. 
Received the Secretary June 29, 1942. 


“Uber die Biegung der allseitig unterstiitzten rechteckigen Platte unter Wirkung einer 
shenko, McGraw-Hill Book Co., Inc., New York, Y., 1940, pp. 156-167. 
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From Eq. 52, with 
comparing Eq. 154 with Eq. 13: 


mai 


m= 


TABLE 
FoR SUPPORTED AND 
CLAMPED RECTANGULAR 
PLATES WITH CENTRAL 
CONCENTRATED LOAD 


CLAMPED 


—0.6174 


—0.4591 
—0.4516 


Notes: =ratio length plate 


y-direction length Values 


are taken from author’s paper. 
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Discussions 


Then, 


cosh? 


cos The series the fore- 


going expressions for and converge 
rapidly. Numerical values for and 
for the case central concentrated 


load 


given Table 

For similarly loaded rectangular 
plate, with all four edges clamped, the 
moments any point are found super- 
posing upon the moments for simply 
supported plate the moments due the 
distributed edge couples the clamped 
plate. The solution for 
rectangular plate with central concen- 


and several ratios are 


trated load has been and from this solution may shown that the 

corrections added the place coefficients Eqs. 155 obtain the place 

coefficients for clamped rectangular plate with central concentrated load are: 
m—1 


cosh 


(correction) 


cosh 


with values given the reference cited.* 


and and B,, are numerical coefficients 


Applied Mechanics, Vol. No. September, 1939, A-114; also Vol. No. March, 1940. 


and 

and 

ca 
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Numerical values for the place coefficients for the clamped plate are also 
given Table 


stresses plates loaded under transversal forces which are distributed over 
comparatively small areas deserves studied most carefully engineers. 
Those interested the practical applications the mathematical theory 
bending elastic plates must welcome this concise study the effects 
concentrated loads. established fact that the foundations for the 
handling the problems the theory flat plates have long been developed. 
For example, mathematicians have developed means for determining distribu- 
tions certain potential functions space over areas, the theory which 
can furnish the basis for solving many plate problems. Despite these early 
developments (and also valuable recent advances made engineers) for solving 
such problems, the practical applicability the bending theory elastic plates 
still somewhat limited certain facts. the intention this discussion 
point out least few these facts, some which were revealed the 
paper. 

The methods available for computing stresses flat plates loaded con- 
centrated forces belong this group questions. must noted that 
certain difficulties that are encountered when the maximum bending stresses 
such plates must evaluated are not based some inherent weaknesses 
the plate theory. Rather, they are due the fact that engineers dealing with 
these problems perhaps have not yet applied the special methods which are 
needed for solving such cases. The contrary the truth concerning the afore- 
mentioned difficulties—namely, that the theory able furnish the required 
and entirely satisfactory means for finding the peak values the stresses 
plates loaded concentrated forces. 

well known that the section which single transverse load 
applied beam the bending moment has discontinuity. the bending 
theory beams this discontinuity does not introduce any difficulties, and the 
bending stresses the vicinity the area which “concentrated load” 
applied are computed from the same formula the stresses more distant 
points. The assumption (never mentioned) that the practical applications 
such load necessarily should distributed over sufficiently large contact 
area, with the adjoining member the construction transmitting the load, 
that further around the pressure area need considered 
most cases. 

When the corresponding case concentrated load occurs plate, 
engineers can, and should, take somewhat the same attitude toward the plate 
theory. The extreme cases very strong load concentration small area 
can excluded from consideration since they also introduce practical difficulties 
construction. Dean Westergaard states clearly, one should distinguish 
whether the concentrated load acts the interior the surface the 
near along one the edge sections. These two cases need somewhat 


Mechanical Engr., Westinghouse Research Laboratories, East Pittsburgh, Pa. 
Received the Secretary July 15, 1942. 
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different methods attack. The present discussion has been confined only 
the first case. 

Few instances are known the writer which authors who have dealt 
with similar cases have gone the trouble establishing the workable expres- 
sions needed for expressing the distribution the bending stresses the vi- 
cinity and within the area through which the concentrated load acts. Although 
the form solution Lagrange’s equation (Eq. 23) for load which 
concentrated mathematical point known the vicinity this point, 
other important properties the solution rarely have been visualized dis- 
cussed those who make use these developments practical applications. 

The static feeling engineer tells him that should not, and does not, 
matter much how load actually distributed over area, which sufficiently 
small (not ‘‘too when compared its dimensions with the thickness 
the plate and not large’’ when compared with the smallest lateral di- 
mension the plate, measured parallel its plane), one wants determine 
the stresses only points the plate, which are farther away from the loaded 
area. Practically very little use has been made this good instinctive 
feeling the mathematical applications the plate theory. other words, 
useful and helps the discussion the case the load replaced first 
single force concentrated mathematical point. Clearly this 
mathematical abstraction, but very useful one, among the many other ones 
used, provided one confines its proper limits. 

Engineers probably have not been conscious the fact that finding the 
proper integral Lagrange’s equation the plate theory under the given 
boundary conditions which corresponds the case the single force concen- 
trated point the first necessary step toward establishing the expressions 
required for computation the peak stresses also near the finite loaded area. 
Experience has shown that the deflections and the form the elastic surface 
can computed any point, including the area loaded, using this solution. 

the mathematical theory the distribution potentials around given 
bodies areas, use made similar ideas when potential functions are intro- 
duced which (1) satisfy the differential equation the potential, (2) become 
finite single point required manner, and (3) satisfy the prescribed 
boundary conditions along the surface the body along curves the case 
two-dimensional problems. These functions are there known 
which solve special cases boundary conditions. The bending 
theory flat plates closely related the theory logarithmic potentials. 
Evidently, finding the solution for given contour line plate under given 
boundary conditions and carrying single force arbitrary point identical 
with this problem finding the Green function Lagrange’s plate equation 
the fourth order, satisfying given boundary conditions. The computation 
the stresses near within the area which concentrated force acts 
plate thus must consist first determination the Green function corre- 
sponding the load concentrated the center the loaded area. 

After the solution for single concentrated force the center the loaded 
area has been found constructed the methods which are available for 
such cases the three and introduced the author, 
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can computed readily comparing with Eq. known, however, 
that logarithmic tube’’ which expresses the bending stresses plate 
near the point application the load increases infinity near this point, 
and expressions based the load concentrated point fail and near the 
point application. The author commended for showing clearly that 
the most general case three additional functions, and must 
constructed defined Eqs. and 12. These functions are needed the 
bending stresses are computed near within the loaded area. They need 
computed only once for given pressure distribution and shape the 
loaded area and are not dependent the shape the contour line the plate. 
After these three functions have been constructed, the peak stresses, which are 
usually found within the loaded area the plate, can computed. 
interesting note that few cases the plate literature evidently these 
integrals have been constructed implicitly when the complete solutions were 
found the ordinary methods integration. 

noted that the author has made frequent use expressions containing 
functions the complex variable The theories plane stress 
and the bending plates are both based the biharmonic equation and 
closely related the theory functions. plane hydrodynamics when 
certain patterns potential flow are discussed the same methods are used. 

this connection perhaps the attempts his doctor’s 
thesis 1928, may deserve mention. used the methods conformal 
mapping for the numerical construction the bending moments parallel 
strips and rectangular plates loaded single force and satisfying the Navier 
boundary conditions. efforts were directed toward the numerical de- 
termination Green’s function for the moment sum. The results his 
elaborate computations, which were presented curves plotted for the bending 
moments and and for the twisting moment should become even 
more valuable today the three place coefficients, and introduced 


Dean Westergaard, were computed for the various cases which were discussed 
this thesis. 


Correction for Transactions: April, 1942, Proceedings, page 516, Eqs. 
” ” 
and 39, change three places. 


mit freiaufliegenden rechteckigen Platten unter Berg- 
Forschungsarbeiten, Heft 302, 1928. 
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ROBERT FOLLANSBEE, AND DAVENPORT 


Am. Soc. the early contri- 
butions Mississippi River hydrology were made chiefly the Army en- 
gineers, the Navy made its contribution even before the classic work Hum- 
phreys and Abbot. Acting under instructions from the Secretary the Navy, 
Lt. Robert Marr, master the Memphis Navy Yard, investigated the 
discharge the Mississippi River Memphis from April July, 1849, and 
from March, 1850, February, 1851. made measurements means 
floats, and deduced the mean velocity from surface observations means 
floats sunk various depths. The reduction factor was slightly more than 
10%. also obtained daily gage heights, and from these gage heights and 
the results the discharge measurements computed the daily discharge for 
the period covered. description this work was presented the 
ington Astronomical and Meteorological Observatory.” 

The author refers the earnest approaches made the early contributors 
Mississippi River hydrology toward the more rigorous standards developed 
during recent years. that connection may noted that Humphreys 
and Abbot’s study the vertical-velocity curves obtained during their investi- 
gations indicated that the velocity different depths varied the abscissas 
the parabola which had its z-axis parallel and below the water surface. 
However, they did not pursue further the study the properties the para- 
bolic curve. That was reserved for British engineer, Maj. Allen Cunningham, 
who, during the 1870’s, made 565 sets vertical-velocity curves the Ganges 
Canal India. concluded that, considering the vertical-velocity curve 
common parabola, its properties indicated that the mean velocity equals 
(velocity 0.211-depth velocity These are practically 


paper Jarvis, Am. Soc. E., was published March, 1942, Proceedings. 
Dist. Engr., Geological Survey, Denver, Colo. 
Received the Secretary June 26, 1942. 
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the 0.2-depth and the 0.8-depth. was not until 1906, however, when 
Barrows, Am. Soc. E., then member the Geological Survey, 
made study vertical-velocity curves and found that their form was that 
parabola, the mean ordinate which was closely the average the ordi- 
nates the 0.2-depth and the 0.8-depth, that the procedure measuring the 
velocity these depths and taking the average the mean velocity was 
adopted the Geological Survey its standard. 


early records the Mississippi River constitutes valuable service furnishing 
the flow this dominant waterway for period con- 
siderably more than century. 

The long-time compilations runoff presented the paper are subject 
use variety ways, not only for fixing the mean runoff over long period 
years but for studies trends, and variation magnitudes, for possible 
correlation with other factors. The author has included comparisons devi- 
ations runoff with those rainfall. also has mind that the data may 
form foundation which assemble, compare, and interpret land-use 
information. 

particularly desirable such uses runoff data that the relative re- 
finement the data known that one may consider the extent, any, 
which evidence relationships, apparent sought for, may influenced 
fortuitous systematic errors. not uncommon for such errors obscure 
invalidate the significance the evidence. desirable, the author 
indicates, suggest modifications that will make the results more nearly what 
required for such uses. also very essential know the inherent limi- 
tations the results that they may not applied unwisely. the data are 
satisfactorily free such errors, the further question should considered 
whether the significance the factor investigated, such land use, may 
masked other factors, particularly those climatic variation. 

The writer has been greatly impressed with the tremendous quantity data 
concerning the history and behavior the Mississippi River and feels quite 
humble discussing subject which his knowledge encompasses compara- 
tively little. The points raised herein are not intended the way 
criticisms, but more the way questions which hoped the author may 
throw light. 

The longest record developed that the Mississippi River Natchez 
from 1817 1938, period 122 years, the results being shown Tables 
and The technique used the determination the record the earlier 
years may learned perusal Professional Paper No. which presents 
results for out the years from 1818 1860. These records apparently 
were based stages prior 1848 solely Natchez, and subsequently mostly 
Carrollton and Donaldsonville, La. There seems very little information 
for appraising the accuracy the basic stage records. 


Hydr. Engr. Chg., Div. Water Utilization, Geological Survey, Washington, 
Received the Secretary July 29, 1942. 
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The stage-discharge relation for the Natchez station was determined from 
discharge measurements January and February, 1858, Natchez, and for 
the remainder the year was determined Vicksburg and transferred 
Natchez. developing the relation mean curve was apparently drawn 
through the loops described the individual measurements rising and falling 
stages. The inaccuracy involved the determination the discharge for 
given day would tend compensative computing the discharge for the 
year. The mean annual discharge was determined essentially application 
the twelve mean monthly stages the year the rating table, adding the 
results, and dividing twelve. 

The discharge determinations consisted velocity measurements made 
means timing double floats, the most common form being surface float 
cork which was attached paint keg, the depth submergence which was 
fixed cord attached the cork. The kegs were passed all depths from 
the bottom the top, but uniform depth seems have been adopted 
for general use. The mean velocities the floats were determined for divisions 
200 width, and the products the mean surface velocity and the areas 
the cross sections for such divisions were computed and added for the entire 
cross section give the Soundings were made the 
vicinity Natchez January and the vicinity Vicksburg February, 
September, and December during 1858, which apparently gave the basis for the 
areas cross sections used the discharge computations. 

The was corrected ratio designed adjust 
the float velocity the mean velocity the section and also take account 
the prevalent condition the wind. Similarly, the discharges Carrollton 
and Donaldsonville were computed from measurements made 1851 and 1852. 

This early river measurement technique, briefly summarized herein but 
described fully Professional Paper No. was unquestionably the best 
afforded the knowledge and facilities the time and was conducted with 
supreme care. may compared with that used the present time the 
Geological Survey for the gaging station 


Measurements.—Discharge measurements are made from 
the downstream side the highway section the bridge. The highway 
section the upstream side the bridge, and the railway section 
the downstream side. The metering equipment lowered through 
opening between the floors the highway and railway. The equipment 
operated means reel and boom mounted low-wheeled truck 
frame. The truck frame also carries four-cylinder automobile engine 
which supplies power for operating the metering equipment and moving the 
truck from point point along the bridge. Depth water and position 
meter are shown depth indicator attached the reel. During 
medium and high stages the depth water and velocity are such that the 
metering equipment carried downstream and the indicated depth too 
large. Conversion true depth effected use two correction tables, 
one which gives the correction slant length supporting line above 
water surface, and the other gives the correction curved length line 
below water surface. The use each table requires measurement the 
vertical angle the supporting line above the water surface, and this 


Water-Supply Paper No. 747, Geological Survey, Washington, C., 1933, 16. 
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measurement accomplished use protractor applied the line 
just below the bridge floor. The corrections indicated depth are 
lessened the use large sounding weights. The maximum sounding 
weight used the Vicksburg station 300 pounds. 


control method. Changes stage-discharge relation determined 
frequent measurements.”’ 


From the establishment the Vicksburg gaging station 1930 the end 
the water year 1941, 1,524 discharge measurements have been made this 
station, generally two three days apart, record the vagaries the flow 
the river. Fig. illustrates the plotting discharge measurements 1938 
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excess 700,000 per sec. The writer unable state how typical this 
period may have been. However, the plotted points suggest the complications 
involved the determination the stage-discharge relation for application 
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the stage records. Soundings are made for each discharge measurement and 
the cross sections thus determined indicate considerable variation during short 
periods. Comparisons areas cross sections for gage heights 
from soundings from 1933 1939 indicate increase area about 10%. 
This increase may have been associated with the effects cutoffs below 
Vicksburg. However, illustrates another changing aspect the behavior 
the Mississippi River. 

The application the stage-discharge relation determined for Natchez 
1858 involves questions reliability stage records, the accuracy the 
compensative adjustment use mean curve drawn through rising and 
falling branches the plotted measurements, uncertainty cross-sectional 
areas applicable the time measurement, and other matters pertinent 
pioneer technique. Much more serious than these questions, however, are 
those pertinent the application the 1858 stage-discharge relation preceding 
the period which the determining observations were made far back 
1818. Illustrations Improvement the Lower Mississippi River for 
Flood Control and delineate pronounced progress channel 
changes the Mississippi River below Vicksburg and Natchez between 1820 
and 1861 and also, apparently, the construction levees. 

Fig. which shows rating curves the Mississippi River Vicksburg from 
1858 1937, indicates that that period the mean rating curves that have 
been defined have varied within limits approximating +5%. apparently 
not demonstrated, however, that changes influential the Vicksburg-Natchez 
section the river may not have been substantially more operative the 
period 1817 1858 than the later period. 

The critical feature respect the use the earlier records the cloud 
uncertainty resting over their reliability. The author recognizes this situation 
his supposition that modifications may suggested that will make the 
results more nearly what required for the analysis land-use information. 
notable that the Mississippi River Commission has not included such 
modifications its compilations past records—for what reason the writer 
not informed. their present status, important caution earnestly that 
any applications the records their inherent nature respects accuracy 
kept mind. 

Improvement the Mississippi River for Flood Control and Navigation,” prepared 


Elliott, under the direction Jackson, Waterways Experiment Station, Vicksburg, 
Vol. III, May 1932, Plates VI(f), VI(h), and XXXII. 
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paper for evaluating the economic benefits from proposed projects intended 
wholly partly for flood prevention. The applicability the procedure ob- 
viously depends upon the purpose that the investigation serve. 
evaluation flood-prevention benefits conceivably might made for large 
landholder, group landowners, determine whether the construction 
works for the protection the property would profitable private 
enterprise; might undertaken for municipality flood control dis- 
trict determine whether the local benefits would justify the local taxpayers 
undertaking flood protection works; may undertaken determine 
the economic justification from the viewpoint the adoption federal 
flood control multiple purpose project. The procedure presented the 
paper evidently has view this third purpose. 

prudent, even imprudent, private owner would moved principally 
the anticipated enhancement the value his property deciding upon 
the advisability investing his own funds flood control works. This en- 
hancement might credited with the flood losses prevented, which are listed 
Item 2(a) the paper, well with the more profitable utilization the 
land after protection, listed Item 2(b), but the latter usually the com- 
pelling consideration. any case, would scarcely undertake any consider- 
able expenditure unless least the paper profits showed wide margin over 
the estimated cost. The close figuring shown the author’s example ordi- 
narily would not have much financial appeal. 

Flood control works, whether the nature storm drains more exten- 
sive enterprises, ordinarily are undertaken municipalities local flood 
control districts broad grounds public necessity convenience, rather 
than because detailed evaluation financial benefits. Generally, such 
works are adopted either rectify inconvenient intolerable situation 


paper Otis Page, Am. Soc. E., was published March, 1942, Proceedings. 
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that has grown because municipal development, prevent the recur- 
rence recently experienced flood disaster. Such works are built satisfy 
insistent public demand, and close figuring financial benefits not 
often required. exceptional cases, sure, legal requirements may entail 
the assessment benefits against the properties affected, and the considera- 
tions presented the author might applicable the preparation the 
assessment roll. 

the investigation proposed federal flood control project, evalua- 
tion the benefits required the Flood Control Act June 22, 1936, 
which Congress laid down the policy: 


that the Federal Government should improve participate 
the improvement navigable waterways their tributaries, including 
watersheds thereof, for flood control purposes, the benefits whomso- 
ever they may accrue are excess the estimated costs and the lives 
and social security people are otherwise adversely 


Sound policy undoubtedly requires that public works undertaken only 
when the benefits exceed the cost; but these benefits are not neces- 
sarily, nor even generally, material ones which can evaluated en- 
gineering computation. The more important works contemplated the 
Federal Flood Control Act are those for the prevention mitigation 
recurrence destructive floods which have caused loss life well 
property, have produced widespread inconvenience and discomfort because 


the interruption rail and highway traffic and other public services, and have 
inflicted suffering and misery upon the considerable number persons who were 
driven from their homes. The writer does not believe the intent 
Congress that the justification for such projects should rest nice evaluation 
the material benefits, but rather that the rule reason applied de- 
termine, general grounds, whether the value benefits which cannot 
definitely computed money will exceed the estimated cost. 

the other hand, the original Flood Control Act included the authoriza- 
tion number projects intended remove continuing flood hazard 
agricultural lands which the local inhabitants had exhausted their resources 
attempted protection, were unable raise funds for the necessary works. 
The procedure described the author would appear particularly applicable 
the investigation proposed projects this character. 
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Discussion 


FREDERICK MERRITT, AND RALPH STEWART 


continuous structures moment balance has two main features reeommend 
it. First, enabling designers estimate final moments, considerable work 
may Second, its self-checking character may also save time, that 
errors may either immediately located neutralized further balancing. 
Accounting for sidesway moments much handicap other methods 

The method moment balance may applied somewhat different 
form with the advantage ease generalization and saving computation. 
Let: six times the rotation one end simply supported beam unit 


times the rotation one end this beam due unit moment applied the 
other end; the ratio least moment inertia span length; actual 
end moment minus fixed-end moment; and the rotation the end any 
beam. Clockwise moments will considered positive and counterclockwise 
moments negative. Then, referring Fig. general, neglecting deflection 
beam ends: 


Kon 


6 E Kop go = Mop Bop Mpo 

and 


paper Cornish, Esq., was published May, 1942, Proceedings. Discussion 
this paper has appeared Proceedings, follows: May, 1942, Messrs. Matthews, and William 
Larsen; and June, 1942, Messrs. Eremin, Grinter, and Aleck. 
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TABLE SHEET 


Mem- Mem- 


(a) 


+0. +0.6 +2.5 
Final moment 3.4 +2.6 


that the carry-over factor toward for the beam, considering fixed- 
ended. Summing yields: 


which the change the fixed-end moments the far ends. Since 


Substituting Eq. Eq. 14a: 


Using Eq. 17, the method moment balance may applied continuous 
frames containing beams variable moment inertia. The final moments 
are obtained adding the results the last balance the fixed-end moments. 
For beams constant moment inertia, Eq. becomes: 


which 


Using Eq. 18, the frame Fig. was analyzed for comparison results 
(see Table 3). The same changes fixed-end moments assumed the 
author were also taken preliminary estimate. Joint was balanced first. 
One half the change the far ends was recorded for and and the sum 
placed under the sigma. This sum was then added the fixed-end moment 
(in this case zero), and the result, entered Table was distributed, with 
sign changed, and proportion their K-values. The change 
the fixed-end moment was obtained adding the two values entered for 


No. Quantity Mem- Mem- Mem- Mem- Mem- 
ber ber ber ber ber ber ber 
ree K 3 5 8 5 4 5 14 5 3 8 
‘ 
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TABLE 3.—(Continued) 


ber ber ber ber ber ber ber ber ber ber 


each member. the final moment were desired this stage, the addition 
would include the fixed-end moment, and the change moment would not 
entered. will noticed that the end the first balance results agreeing 
with those the author were obtained. 

The method moment balance has the same disadvantage practical 
application other methods analysis converging increments—that is, 


locating loads for maximum moment. Because the unknown end condi- 


tions interior spans, difficult apply these methods problems 
moving loads. 


states that the trial and revision procedure which describes rapid and 
that the arithmetical work almost entirely self-checking. These claims are 
attractive and led the writer experiment some with the author’s demon- 
stration problem shown Fig. 

Table shows the results three independent solutions this problem, 
well the final moments computed the author. For the practical 
design building frame this type necessary consider various com- 
binations loading, including wind earthquake loads. This requires that 
the frame solved for the effect fixed-end moment each its joints. 
The first solution Table gives the effect each individual fixed-end 
moment and carried accuracy 0.01 ft-ton. This first solution was 
performed without the use simultaneous equations and without the use 
any kind successive approximations. 

The solution Springer (see Fig. was performed conventional 
end-moment distribution min, including the time required make the 
sketch and write the stiffness factors and fixed-end moments. The solu- 
tion Blume (see Fig. 10), which the same order accuracy the 
author’s two-place solution, required only min after the sketch showing 
stiffnesses and fixed-end moments was prepared. 


Engr. Bridge and Structural Design, City Los Angeles, Los Angeles, Calif. 
Received the Secretary July 1942. 
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TABLE Four INDEPENDENT (see Fig. 


Description 


Effect fixed-end moment, AB. 
Effect fixed-end moment, BA. 
Effect fixed-end moment, FE. 


Combined moment............ 


Joint A Joint B | Joint C 
+0.01 +0.38 —0.77 +0.38 +0.01 
{ 
3.0 0.8 2.2 0.1 0.1 
3.40 3.09 0.88 2.21 0.01 0.01 


+184 


0.357 0.357 
8 
-038 
+058 +031 
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TABLE 4.—(Continued) 


—1.29 +0.57 +-0.72 +-0.16 +0.27 —0.22 —0.21 —0.02 —0.05 +0.07 
—0.07 +0.06 +0.01 —0.11 +0.36 —0.13 —0.11 +0.01 —0.07 +0.06 
+0.34 +0.58 —0.92 —3.16 —2.18 —2.99 +8.33 —0.92 +0.34 +0.58 
—0.10 —0.21 —0.33 +0.89 +1.15 +1.70 
—0.99 +1.28 —1.76 +8.90 +1.37 +2.41 

( 1.3 0.3 3.4 1.8 3.6 8.8 3.8 1.4 2.4 
0.99 1.26 0.29 3.45 1.76 3.68 8.89 3.79 1.35 2.42 

1.0 1.3 0.3 3.4 ef 3.6 8.9 3.7 1.4 2.3 


The author deserves credit for having made scholarly investigation the 
analysis rigidly jointed frames. However, between the superior speed 
the conventional end-moment distribution solution and the superior quality 
the first solution Table appears the writer that the author’s method 


0.375 0.625 0.36 0.28 0.36 0.625 0.375 

0.31 0.23 0.46 0.285 0.19 0.285 0.24 0.46 0.23 0.31 


also may mentioned that for one-story structures the type 
solution will considerably faster than either end-moment distribution 
the method. 


Correction for Transactions: May, 1942, Proceedings, page 684, Eq. 
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DISCUSSIONS 


PROFILE CURVES FOR OPEN-CHANNEL FLOW 


Discussion 


Assoc. Am. Soc. reading this paper, one 
wavers between two conclusions—that Professor Gunder wrong, that the 
backwater curves shown standard texts are seriously error. Fortunately, 
neither conclusion correct. The principal difference lies the shape the 
type curves channels mild, zero, adverse, critical slope. According 
Fig. the and curves start perpendicular the bottom the 
channel and are concave downward for some distance before becoming concave 


upward. Previously published sketches these curves have shown them 
dicular the bottom and concave downward all the way the critical depth 
line. Heretofore has been shown simple horizontal line. 

Actually, the discrepancies are not serious one might infer from exami- 
nation Fig.1. The author himself states that the point inflection lies near 
the bottom the channel—very close the Bazin Ganguillet-Kutter 
expressions are used for the Chézy coefficient. Even the Manning relation- 
ship used the discrepancies are not important. drawing small-scale 
illustrations such Fig. the true shape the curves often must modified. 
That Fig. has been somewhat distorted shown the fact that the curves 
start from the bottom right angles it, whereas, according the text 
the paper, they should start vertically. This peculiar feature the curves, 
which the most surprising result Professor Gunder’s studies, and which 
makes the reversal curvature necessary, can easily shown consistent 
with the assumptions ordinarily made computing backwater curves. 

The rate loss head through friction, with respect distance along the 
channel, assumed given Chézy’s formula 


paper Dwight Gunder, Esq., was published April, 1942, Proceedings. Discussion 
this paper has appeared Proceedings, follows: June, 1942, Carl Rohwer, Assoc. Am. Soc. 


Associate Engr., Iowa Inst. Hydr. Research; Associate Prof., Hydraulics and Structural Eng., State 
Univ. Iowa, Iowa City, Iowa. 


Received the Secretary August 1942. 
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the range infinitesimal depths close the bottom the channel the 
velocity head very great, and changes rapid rate, whereas the elevation 
head scarcely changing all. rectangular channel the rate change 
velocity head, with respect distance along the channel, given 


(21) 


Since the change elevation head comparatively insignificant, the friction 
loss derived almost entirely from the change velocity head. Equating 
the values given Eqs. and 21, and changing signs allow for the implicit 
negative sign one finds that for depths close zero 


From this expression evident that any formula used which ap- 
proaches zero approaches zero, the corresponding water-surface slope will 
infinite. The limit entirely theoretical, course, representing flow con- 


‘dition impossible attainment. 


The evaluation friction the region near this limit, where the velocity 
flow many times the Belanger critical, very uncertain. customary 
assume that the rate loss head friction, variable flow, the same 
uniform flow the same depth and velocity. The rate loss uniform flow 
then computed one the standard open-channel formulas. doubtful, 
however, whether any open-channel formula has proven range extending 
velocities greater than twice the critical. Observations extremely high 
velocities are difficult, and flow slopes steep enough sustain such velocities 
almost certain become pulsating flow rather than uniform 
Although uniform flow the high velocities existing the lower portions the 
type curves nearly impossible, variable flow this range can exist. Ex- 
amples computed profiles for common case are shown Fig. dis- 
charge 507 per sec flows out from under gate into channel with 
horizontal bottom. Flow from left right the figure, and two possible 
conditions flow starting from different gate heights areshown. The Manning 
formula used, with taken 0.016. the bottom the gate above 
the horizontal bottom the channel, the water backed depth 
behind the gate. the bottom the gate above the bottom the 
channel, the depth upstream becomes either case the water-surface 
profile computed Eq. follows the same curve, crossing the critical depth 
vertically point 3,370 downstream from the gate when the opening 
ft, 6,250 downstream when the opening ft. 

Water-surface profiles computed using the constant value which corre- 
sponds the median value are also shown The curve downstream 
from the 10-ft gate reaches the critical depth distance 3,500 ft, 
excess the value given Eq.19. The curve starting the 4-ft gate reaches 


Transactions, Am. Geophysical Union, 1940, Pt. 594. 
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the critical depth 6,760 downstream, distance greater than that given 
Eq. 19. Although evident that the discrepancy would greater for 
smaller gate openings, the 4-ft gate opening already beyond the range 
values likely encountered practice. The discrepancies could de- 


Ft) 


2000 4000 
Along Channel, Feet 


Constant 


rather than the median value. For example, with the gate ft, value 
corresponding depth gives result differing from that computed 
less than 1%. 

With regard simplicity, not much easier use constant value 
computing the curves than use Eq. 19, but the use constant value 
computing and curves permits the use Bresse’s with 
great saving labor over that required using Eqs. and 18. 

From practical viewpoint, there seems little choice between results 
obtained using variable Chézy coefficient and those obtained using constant 
coefficient. The effect channel shape will ordinarily greater importance. 
Until comparison has been made with data from variable high-velocity flow 
wide rectangular channels, neither method can considered have superior 
accuracy. The author’s prediction reversal curvature the type 
curves, phenomenon which may observable channels any shape, 
should help stimulate interest obtaining actual measurements. 

Tables Bresse’s function for adverse well mild and steep slopes are given 


Steady Flow Open Sherman Woodward and Chesley Posey, John Wiley Sons, New 
York, Y., 1941. 
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Founded November 1852 


DISCUSSIONS 


NUMERICAL PROCEDURE FOR COMPUTING 
DEFLECTIONS, MOMENTS, AND 
BUCKLING LOADS 


Discussion 


JOHN WILBUR, RALPH STEWART, 
AND STEFAN FRAENKEL 


Joun Assoc. Am. Soc. Newmark’s 
paper has been stimulating the writer, particularly since places emphasis 
mathematically simple approaches problems elastic stability. 
though Professor Newmark shows that simultaneous equations may used 
insure that the assumed deflection curve compression strut essentially 
proportional its shape computed elastically the basis the assumed 
deflection curve, his detailed treatment devoted the substitution 
method successive approximations for the solution simultaneous equations. 

The writer solved the problem illustrated Fig. the basis five 
simultaneous equations, and obtained substantially the same results. This 
solution was relatively straightforward, although the form the simultaneous 
equations was such that was convenient adopt method successive 
approximations their solution. Having successfully completed this solu- 
tion, was then decided investigate the accuracy simpler direct solutions, 
based two simultaneous equations only. Symmetrical cases only were 
considered; otherwise three equations would have been necessary. The re- 
sults these relatively simple solutions were quite satisfactory, illustrated 
the following discussion. 

Consider first the case Fig. 13. Assuming that the deflections the 
strut the center and quarter points are and respectively, the equivalent 


concentrated the center and quarter points are computed from 


the relations for parabolic loading given Fig. 5(b). These loads, together 
with the resultant net reaction, are shown Fig. the moment area 


paper Newmark, Assoc. Am. Soc. E., was published May, 1942, 
Proceedings. Discussion this has appeared Proceedings, follows: May, 1942, Bruce 
Johnston, Assoc. Am. Soc. E., and June, 1942, Ketchum, Jr., Assoc. Am. Soc. 
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method: 


and 


and 


About Center Line 


Symmetrical 


PLs 


Before applying this the case illustrated Professor New- 
mark Fig. 15, where not constant throughout the length the strut, 
was convenient develop equations for equivalent concentrated loads for 
parabolic loading curves for the case which the spans the two adjacent 
segments are not equal. Referring Fig. 25, these equations are follows: 


and 


A 
| 
= 
t 
q 
| 
~ 
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Assuming that the deflections the strut the center and the 0.2-points 
are and respectively, the equivalent concentrated -load the 0.2- 


point computed from the foregoing relations for parabolic loadings, whereas 


the equivalent concentrated -load the center computed the basis 


parabolic loadings equal adjacent sections. These loads, together with the 


Symmetrical About 
Center Line 


resultant net reaction, are shown Fig. 26. the moment-area method, 


and 


P 
(ED 
cr | | 
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and 


Since, for strut composed sections constant the deflection curve 
must compounded sine curves, interesting exact solution the fore- 


going problem may made follows, with reference Fig. 27. 


Point 
Common 
Tangency 


For the case under consideration, the true deflection curve, shown the 
full line, composed segments the curves 


41 


and 
the point common tangency, where 


=> 10 (5 Ly — 3 Ls) 


and 


ad 
Xo \ 
a . 
the following relations hold: 


6b) 
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These relations reduce 


(5 3 L3) aw Ls 


and 28b become 
6; sin — (5 Li — 3 Ls) = 


(29a) 
and 


Dividing Eq. 29a Eq. 


(5 L; — 3 Ls) 


1 


This gives the length 
the strut under consideration. The critical load then given 


4.50 


the opinion the writer, Professor Newmark has performed definite 
service since has dissociated one type problem elastic stability from 
the shroud the formal solution differential equations. hoped 
that relatively simple methods solution will eventually developed for 
other types problems elastic stability. 


determining the elastic curves beams are presented compact form the 
author. The paper has value engineer’s reference library since the 
analyses demonstrated have heretofore been scattered different treatises. 
The use progressive load, shear, and moment increments establish deflec- 
tions and the alinement elastic curves familiar the analysis 
excellent demonstration the use the device termed correction 
(which, analogy, the same “linear correction deflec- 


Engr. Bridge and Structural Design, City Los Angeles, Los Angeles, Calif. 
22a Received the Secretary July 1942. 
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appears the writer that part the paper dealing with beams 
variable section can improved analytical procedure and also that rather 
important deficiency the author’s illustration the use his computations 
justifies rewriting this part the paper. 

clarify and justify this opinion, statement regarding the basic constants 
beam flexure necessary. basic constant defined being either 
simple constant quantity constant ratio between two variable quantities. 
constant that more complex than basic constant will referred 
derived constant because derived the use two more basic constants. 

The basic constants beam flexure are well illustrated the dimensions 
governing railway curves. circular railway curve completely determined 
the length tangent and the angle intersection the tangents. The 
transition spiral often used the end circular curve fully determined 
the length its tangents and their angle intersection. the over-all 
length the spiral curve known (as the span beam known), then one 
tangent and the angle intersection are sufficient. The basic constants 
ing the flexure beam that has more than two supports are similar and 
are follows, the distance between supports being known: 


(1) The basic stiffness the left end the beam. This the ratio 
the moment the angle intersection between the tangents the elastic 
curve when the beam hinged the right end and moment applied 
the left end. 

(2) The ratio the length one tangent the elastic curve the length 
the beam for the same condition flexure (1). 

(3) and (4) These are the same (1) and (2), except that the left end 
the beam taken hinged and the moment applied the right end. 

(5) The angle intersection between the tangents the elastic curve 
when the beam acts simply supported beam subjected its loads. 

(6) The ratio the length one tangent the elastic curve the length 
the beam for the simply supported condition. 


” 


The first four constants are which are independent the 
loading. The last two are “load which depend the loads. 

the six basic constants flexure, only five are independent, since the 
principle Maxwell’s theorem reciprocal deflections (angular) will enable 
any missing beam constant computed from the others. understood 
beam flexure that the tangents may taken equal their projected 
length the unsprung beam. 

end slope not basic flexure constant. This can seen from the 
fact that cantilever beam the slope the fixed end known zero; 


this knowledge, combined with the length tangent the elastic curve, 


insufficient determine the curve. If, however, the length tangent 
and the intersection angle between tangents are known, the curve determined. 

The Hardy Cross stiffness factor not basic constant because 
function end slope and also two basic stiffnesses, one each end the 
beam. 


Fj 
i 
7 
; 
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Much valuable time has been lost and still being lost engineering 
designers through injudicious selection flexure constants analyzing the 
various types structures. The author’s treatment the beam Figs. 
and would tend encourage rather than correct this tendency. revised 


M=180 
Concentrated Angle Changes 696 808 414 231 123 


Loads About Left End 


! 
Distance from Left End See Diagram Below. 


' 


flection=b-a= x 


Spread =b= 3129 2608 2086 1565 1043 5216 


' 
‘ 


J 
Deflection 


Basic Flexure Diagram. 
Draw One for Each End Beam. 


Basic Stiffness =0.772 (Relative) 


Carry-Over Factor =0.288 Left Right 


0.575 (Any Unit Relative Unit Length May Used) 
Cross Stiffness Factor Stiffness 1.267 
(Carry-Over Factor Right Left 2.548 =0.738) 


computation the properties this beam therefore offered Fig. 28, and 
the justification for this revised computation illustrated Fig. 29. The 
columns Fig. are the same members whose properties are determined 
Fig. 28. The deck constant section with relative stiffness shown. The 
sidesway moments and the appurtenant lateral force shown Fig. are 
computed the following consecutive steps: 

Write deck traverse angles and Multiply stiffness, 2.97, 
angle value, get moment 8.91; divide moment 8.91 stiffness 1.485 


j 
1.342 
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get column traverse angle Add angles and get bottom column traverse 
angle multiply angle stiffness 0.772 get moment 5.404. Add all 
column moments and divide frame height get lateral force. The effect 
any other lateral force will direct proportion. The lateral deflection 
the structure follows direct by-product. 


1.485 


can seen that the foregoing computation sidesway moments, using 
basic constants, can done less than one minute, whereas computing these 


moments, with the derived constants, which are the only ones mentioned 
the author, would take much longer, would give less accurate result, and 
would fail give the deflection the structure direct by-product. The 
same situation would prevail with vertical loads. 

The usefulness the computation deflections and elastic properties 
individual beam obviously increased bringing the results into form 
that can used for the quick computation the deflections and moments 
structure which the beam becomes member. 


Jun. Am. Soc. the writer was 
called upon investigate the stresses 125-ft derrick boom. Since such 
long and slender members are subjected considerable deformation which, 
turn, influences the stresses, was necessary determine the deflection this 
boom. The methods outlined Professor Newmark’s able paper lent them- 
selves well this investigation. 

The boom (see Fig. 30) consisted four angles—two 6-in. 4-in. }-in. 
angles, which formed the top chord, and two 4-in. 4-in. angles, 
which constituted the bottom chord. preliminary investigation indicated 
that the strengthening the middle section (which was long) with 3-in. 
bars the top and 3-in. bars the bottom was advisable. 
That the section shown Fig. 30, which was also used the investigation 
related herein. The section the boom varied considerably outside the 
middle part, shown the variations the values and 
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general view the boom shown Fig. 31. The boom will analyzed 
horizontal position. Because the fact that its cross section unsym- 
metrical, the neutral axis does not coincide with the geometrical center line, 


48" 


N 


Geometrical Center Line 


24" 
26.955" 


Wx 


Boom 


Neutral 


Line Action Axial Load 


Geometrical Center 
Line (Horizontal) 


31.—GENERAL 125-Fr Derrick Boom 


the distance between them any point being denoted The line 
action the axial loads determined the right end the boom the 
location the pin, and the left end the intersection the boom line and 


—8.39 1.13 2.27 2.95 2.95 2.95 2.95 2.95 2.27 1.76 0.00 
42,057 5,363 7,098 |11,843 |11,843 |11,843 |11,843 |11,843 7,098 3,387 5,071 


the vertical line action the load (It should understood that 
includes both live load and the left reaction the boom due its own weight.) 
This intersection (which, course, imaginary) occurs for the case the 
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flat boom 2.06 in. below the geometrical center line. The distance any 
point between the geometrical center line and the lire action the axial 
load will denoted summary the values and for various 


1600 
1800 


2770 


190.2 


4' 
59! 27' 


125° 


Fig. shows the loading the boom due its own weight. The dead- 
load reaction the left end causes axial load 51,250 lb, and the maximum 
axial load due live load 64,950 computing the deflections the 
following conventions were adopted: 


(1) Downward deflection positive; 

(2) Inch and pound units are employed the computations; 

and 

(4) The boom divided into ten parts, that ft. 


TABLE DOWNWARD DEFLECTION 


Moment inertia 42,057 7,098 11,834 


Moment due weight boom 1.888 3.313 4.313 
469 


Description 


Distributed angle change 

Concentrated angle change 

Assumed average slope 

Trial 

Linear correction deflection 
Resultant deflection 1,747.7 
Resultant deflection inches 1.31 


The total deflection, was considered made three constituent parts: 


(a) deflection caused bending moment due the weight the 
boom; 

(b) deflection due axial dead load; and 

(c) deflection due axial live load. 


| 
1,287.8 855.8 467.0 
2,887.8 3,743.6 
295.4 443.1 
3,183.2 4,186.7 
3.14 
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The following relations exist between these deflections: 


bw 5a = Ap (32a) 
and 


Determination necessary computations are recorded Table 
and further explanations are required. course, would not have been 
necessary compute the equivalent concentrated angle changes, since there 
were discontinuities the and the simplified method men- 
tioned Part the paper could have been used. 

Determination and component the deflection the 
upward movement caused the moment that the axial dead load has about 
the neutral axis. That this movement upward can ascertained from the 
fact that, except the left end, the line action the axial load below 
the deflection This can justified follows: Assume the boom 
fully supported throughout its length that there deflection. Then 
remove the supports except for the reactions, and the boom will deflect under 
the simultaneous action bending moment due weight and axial load due 
weight. The moment arm the axial load will the distance between 
the line action the axial load and the deflected neutral axis—namely, 
the time the intermediate supports are removed, and 
therefore will not depend upon the deflection The moment arm given 
minus sign inasmuch the moment will cause upward, negative, 
deflection. All computations for and the values are given Table 


Pornts: 


Common factor 


11,834 11,834 11,834 11,834 7,098 3,387 5,071 
4.777 4.925 4.741 4.237 3.188 1.760 
404 416 400 358 449 520 1/E 
399.5 411.3 395.7 380.2 445.7 421.5 
4,278.1 3,934.3 3,194.8 2,075.1 509.7 
590.8 738.5 886.2 1,181.6 1,329.3 +1,477.2 
4,801.4 5,016.6 4,820.5 1,839.0 N/E 
3.60 3.7 3.62 3.17 1.38 


Determination and A.—The axial live load also will cause upward 
deflection, since its line action lies below the deflected neutral axis. The 
moment arm the axial live load, however, will much smaller since the 
boom deflects the amount before the live load applied. The arm 
Ap), the minus sign being given because the resulting 
upward deflection. not possible use Eq. estimating the 


313 
467.0 
the 
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TABLE 4.—DETERMINATION UPWARD DEFLECTION DUE 


Description 


Moment inertia 42,057 7,098 11,834 


Assumed deflection due axial dead load —0.06 
Moment (51,250 times moment arm) 

Distributed angle change +29.10 

Assumed average slope .00 —50.90 

Trial (1) —303.24 
Trial deflection (2) +2.42 +2.42 
Linear correction deflection +223.88 +201.51 +179.12 
Resultant deflection —172.60 
Resultant deflection 0.075 0.130 


DEFLECTIONS 
—0.130 —0.164 
Loap 


2.26 2.98 


TABLE UPWARD DEFLECTION 


Description 


Moment inertia 42,057 5,363 11,834 


Assumed due axial live load —0.08 —0.078 —0.067 
Moment (64,950 times moment arm) 411.0 —116.2 
Distributed angle change +21.65 +15.90 
Assumed average slope —50.00 —28.35 
Trial (1) +45.56 —32.79 —45.24 
Trial deflection (2) —0.81 +1.81 +1.33 +0.68 
Linear correction deflection —40.27 —35.80 —31.32 
Resultant deflection —42.90 —67.26 —75.88 
Resultant deflection —0.032 —0.050 —0.057 


1258 
—4.55 
—233.4 
+19.68 
+1.64 
+156.73 
217.57 
0.164 
Line 
2 = 
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Pornts: 


Common factor 


11,834 11,834 11,834 7,098 3,387 5,071 
+17.99 +17.08 +16.09 +30.72 
+16.75 34.74 1.82 +67.91 
—378.06 —361.31 —326.31 —274.75 —206.84 —118.78 
+1.50 +1.42 +1.34 +1.68 +2.56 
242.14 247.85 235.33 206.24 160.38 93.83 
0.182 0.186 0.177 0.155 0.120 0.071 


LEADING THE SAME DEFLECTIONS 


—0.182 —0.186 —0.177 —0.155 


3.42 3.58 3.44 3.02 


Common factor 
11,834 11,834 11,834 11,834 7,098 3,387 5,071 
—0.036 —0.019 —0.026 —0.017 —0.031 
—0.806 —0.419 —0.356 —0.576 —0.377 —0.691 
+4.43 +2.30 +1.95 +3.17 +3.45 +13. 
1.30 0.13 2.43 +4.38 +7.55 11.00 
—49.54 —49.41 —47.18 —42.80 
+0.37 +0.19 +0.16 +0.26 +0.29 +1.11 
—71.59 —64.92 —55.96 —27.61 
—0.057 —0.054 —0.049 —0.042 —0.033 
—0.0574 —0.0534 


| 
: 
— 
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upward deflection, the moment arm not proportional the deflection 
caused the lateral loads. The writer believes that fair approximation 
can made such cases mentally comparing the ratio the moments 
due axial load and lateral load and applying this ratio the deflection 
caused the lateral load. believed that even quick comparison will 
yield estimate which will give the correct answer more than two sets 
computations. Table shows the computations for and the values the 
total deflection 
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torsion pendulum determine the amount resonance expected when 
structure known frequency subjected earthquake. This ingenious 
procedure makes use accelerometer record obtained during the quake. 
torsion pendulum operation based the use displacement curve, de- 
rived double integration accelerometer record, has been used the 

parison results obtained the TABLE 5.—Comparison 


two methods shown Table ERATION RESULTS 
which based the East-West 
ACCELERATION 
component the Helena quake Structure 
The slight differences are 
bly due the choice values for 0.25 0.69 
computing the displacement curve. 0:26 
Since the accelerometer instruments 2.00 0.109 0.06 


are generally arranged put 


Scaled from author’s spectrum curve. Based 
motion starter, which responds integration made the Coast and 


ground movement, small part Survey- 

the record near the beginning 

lost, and there is, therefore, some room for exercise judgment when the 
values these constants must selected. The results obtained the two 
methods appear, however, essential agreement. 


paper Biot, Esq., was published January, 1942, Proceedings. Discussion 
this paper has appeared Proceedings, follows: April, 1942, Messrs. George Rich, Hoff, 
and Merit White; May, 1942, Messrs. Heck, Frank Neumann, and Jacob Feld; and June, 1942, 
Messrs. George Housner, and Homer Hadley. 

Senior Engr., Bureau Reclamation, Dept. the Interior, Denver, Colo. 
Received the Secretary June 1942. 
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gives one the impression that large caliber mathematical artillery has been 
brought bear problem that should have succumbed simpler methods. 
default simpler method suggest, perhaps such criticism ungracious; 
and this method analysis seems quite worth studying its 
own account, entirely apart from the problem which here applied. 
the writer understands it, Mr. Tapley’s method (omitting some the details) 
equivalent the following. has computed the kinematic viscosity 
imaginary, less viscous liquid (which would have produced Froudian 
similarity his model) and the Reynolds number resulting from using the 
actual size and velocity his model, but this imaginary viscosity instead 
the actual viscosity. running tests several different slopes, obtained 
enough variation his mean velocities give appreciable variation the 
Reynolds and Froude numbers and derive formula for the coefficient which 
was only from 11,810 29,970, but very clever mathematical device 
(Eqs. 29), the validity which the writer not position discuss 
present, the range extended the entire field from zero infinity. 
Therefore, the value the coefficient for the aforementioned fictitious Reynolds 
number obtained and the problem solved. The writer wishes make 
just few observations. 

would seem that Eq. only one the forms that Eq. might take. 
For example, might not the equation expressed 


course, the experimental results seem answer this particular question 


the negative, but general must not assumed that (F, always can 
put the form (F) (R). 


paper Tapley, Assoc. Am. Soc. E., was published March, 1942, 
ceedings. 


Associate Prof. Mechanics, Ohio State Univ., Columbus, Ohio. 
Received the Secretary July 20, 1942. 
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the writer, the most interesting point the paper the use 
“momentum velocity instead the usual “discharge 
computing and R—that is, instead defining average velocity discharge 
divided area, the author has defined momentum divided mass. 
This gives about 20% higher velocities the examples the paper. Since the 
ratio for the different columns Table quite constant, would 
seem that this innovation would hardly affect the final result this particular 
problem. However, raises the very interesting question whether, 
general, would improvement over the ordinary usage. would seem 
that could hardly expected ever replace, entirely, the ordinary defini- 
tions and because knowledge the velocity distribution would 
required before either could computed. also raises the question 
whether one should ever use square root twice the 
kinetic energy divided the mass. 


The ratio has been variously denoted and various authors. 


and Members, Am. Soc. E., used Woodward, 
Am. Soc. E., Posey," Assoc. Am. Soc. E., used and 
Hunter and Victor Assoc. Members, Am. Soc. E., 
used Km. The value 1.20, which indicated Table quite high, 
but the true value still higher because, the author admits, has neglected 
the variation velocity the vertical direction. Since this would con- 
siderable, the true value was probably between 1.25 and 1.30—almost 
much for laminar flow. This must due the unsymmetrical cross 
section for, the work Professor shows, any reasonable assump- 
tion the frictional resistance would not account for such high value 
the section were symmetrical. 

should also noted that Eq. and the other equations derived from 
it, special attention required any the velocities are negative (that is, 
the upstream direction). This has been pointed out Professor 

minor criticism the paper that the meaning the definitions 
and not clear the casual reader. Another that, model this 
small, would seem that, among the factors listed Eq. 3a, surface tension 
would have included. fact, the writer’s opinion (without any 
rational reason back up) that the model should have been least twice 
large—that is, scale instead scale. 


Correction for Hydraulic Flow,” O’Brien and Johnson, Engineering 
News-Record, August 16, 1934, 215. 


Fluid Mechanics,” O’Brien and Hickox, McGraw-Hill Book Co., Inc., 
New York, Y., 1937, 271. 


Hydraulics Steady Flow Open Channels,” Woodward and Posey, John 
Wiley and Sons, Inc., New York, Y., 1941, 47. 


Mechanics for Hydraulic Hunter Rouse, McGraw-Hill Book Co., Inc., New 
York, Y., 1938, 54. 


Kinetic Energy and Momentum Correction Factors for Pipes and for Open Channels Great 


Victor Streeter, Civil Engineering, April, 1942, pp. 212-213. 
“Variations Correction Factors for Pipes and Open Channels,” Posey, ibid., July, 1942, 
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BENDING MOMENTS THE WALLS 
RECTANGULAR TANKS 


Discussion 
DANA YOUNG, Am. Soc. 


Dana Am. Soc. comparative data presented 
Mr. Silverman Table are very interesting, and satisfying see that 
the results are reasonably close agreement. drawing the curves 5(a) 
Figs. and attempt was made provide design data that would 
the safe side. The significant fact that the actual moments are considerably 
less than for vertical horizontal beam strip. 

The calculation plates under any except the simpler boundary conditions 
best tedious task, and not expected that the detailed calculations 
should made every designer. Fortunately, not necessary so, 
any more than necessary for designer integrate the equation the 
elastic curve every time wishes know the moments clamped beam. 
table moments for all useful cases what needed; the results given 
represent start this direction. 

paper Dana Young, Assoc. Am. Soc. E., was published November, 1941, 
Proceedings. Discussion this paper has appeared Proceedings, follows: February, 1942, 
Silverman, Assoc. Am. Soc. 


Prof. Applied Mechanics, The Univ. Texas, Austin, Tex. 
15a Received the Secretary August 17, 1942. 
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